[RegCNET] Urgent !!!

[bixq]

Dear all:

If you downloaded the latest RegCM3 version (regcm, regcm_beta and 
regcm SVN) in the last 7 days, please redownload it, there is one mistake 
in it. Sorry for my mistake.

Last week, we got the report from Dr. Zhang Dongfeng about the potential
restart problem related to the tuning parameters (qck1oce and qck1land, 
gulland and guloce) of the large-scale precipitation scheme. We decided
to move the do-loop from param.F into init.F, but unfortunately I put 
into a wrong place (in the restart branch of the if-else-endif), I have 
just found the fixed the error this morning.

If you use an older version (before Nov. 20, 2009), and your parameters
qck1oce=qck1land and gulland=guloce, then I think you can stick to use
your version till you have finished your projects.

Sorry again for my mistake, I promise to be more careful later.

Best Regards,
Xunqiang Bi

On Sat, 21 Nov 2009, bixq wrote:

>
> Hi, DongFeng:
>
> Great, yes, there is a problem there, but not for the default parameters.
> Actually, I just checked the reatart on IBM SP6 computer last week, because 
> the qck1oce=qck1land and gulland=guloce in my run, I got the
> perfect restart.
>
> I'll move the do-loop from param.F to init.F to fix this problem.
>
> Thanks and best regards,
> Bi Xunqiang
>
> On Sat, 21 Nov 2009, zhangdongfeng wrote:
>
>>  Hi Dr. Bi,
>>
>>  I focus on the work of coupling the IBIS into RegCM3. When I set up the
>>  restart, I found the original code has a bug.
>>
>>  In initial run (not ifrest), the qck1(i,j), cugl(i,j), rh0(i,j) are
>>  defined in param.f for both ocean and land according to the values of
>>  satbrt(i,j).
>>  In restart run (ifrest), becaues the satbrt(i,j) (in param.f) is used
>>  before it is read (in init.f), so is not defined and all  equal zero. This
>>  leads to the qck1(i,j), cugl(i,j), rh0(i,j) are all defined by land values
>>  in restart run.
>>
>>  If  qck1oce=qck1land and guloce=gulland, just as the default values in the
>>  model, the restart process is Ok. When different values are used, the
>>  restart run results have problems.
>>
>>  For my run, gulland=0.65 and guloce=0.3 (MIT parameter values), the
>>  differences are found between single run and restart run. So the
>>  satbrt(i,j) should be defined in param.f  for both initial and restart
>>  run.
>>
>>         do j=1,jendx
>>         do i=1,ilx
>>            if (satbrt(i,j).gt.13.9 .and. satbrt(i,j).lt.15.1) then
>>             qck1(i,j) = qck1oce    ! OCEAN
>>             cgul(i,j) = guloce     ! OCEAN
>>             rh0(i,j)  = rh0oce     ! OCEAN
>>            else
>>             qck1(i,j) = qck1land   ! LAND
>>             cgul(i,j) = gulland    ! LAND
>>             rh0(i,j)  = rh0land    ! LAND
>>           end if
>>         end do
>>         end do
>>
>>  Best Regards,
>>  Dongfeng
>>

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   Dr. Xunqiang Bi         email:bixq at ictp.it
   Earth System Physics Group
   The Abdus Salam ICTP
   Strada Costiera, 11
   P.O. BOX 586, 34100 Trieste, ITALY
   Tel: +39-040-2240302  Fax: +39-040-2240449
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~