[RegCNET] Bug in the interactive lake model

[bixq]

Dear RegCM-lakers:

We have just finished the draft version of RegCM-4.0_lake. We'll commit
it later after Nellie do some benchmarks.

Michael and Azar, if you (or other RegCM-lakers) have time, you're welcome
join us, have a test and raise your comments. The package is now at:
http://www.ictp.it/~bixq/RegCM-4.0_lake.tar.gz

I have made the following modifications in the above package:

1. at Terrain step, about 180 lakes are added in subroutine lakeadj
    (mod_fudge.f90), the elevation of lakes have been adjustes, and the
    mean depth of lakes have been written into DOMAIN file;

2. get rid of fort.40, instead, the code reads lake depth from #14
    record of DOMAIN file; There is a limitation here, if the user
    want to consider the variation of lake depth for certain big lake,
    he/she have to write a simple code add the information into DOMAIN
    file;

3. get rid of the swap files (fort.41 and fort.42), instead, I add several
    2D (or 3D) fields and use them as swap, this make things easier;

4. I have just checked the serial run, the code can write LAK output
   (Nellie will do further check). The issue Michael raised, "the
    interactive lake module doesn't work" probably due to Zeng's air-sea
    flux schemes used, and there exists inconsistency among them. We'll fix it.
    You know, by default, Zeng is called after vecbats, but lake is in
    vecbats.

5. The parallal mode requires further work.

To Azar Zarrin: please send me the whole version which you used, so we 
could have a full comparison.  I received those codes you sent me 10 days 
ago, when I compiled them with rest RegCM3, it failed.

Thanks and best regards,
Xunqiang Bi and Nellie Elguindi

On Mon, 20 Sep 2010, bixq wrote:

>
> Dear all:
>
> Thank you for finding the bug !
>
> I am working on the lake model now and hopefully we could find the solution 
> and make it work by next Monday (Sep.27) in RegCM-4.0.
>
> Regards,
> Xunqiang Bi
>
> On Sun, 19 Sep 2010, Michael Notaro wrote:
>
>> To the RegCM3 developers,
>> 
>> It seems that the interactive lake module doesn't work
>> correctly in either serial or parallel mode.  In both
>> cases, the precipitation values become zero in the lake subroutines
>> and there is no snow on ice given the lack of precipitation.
>> Could you developers of the model please comment on the problems
>> of the lake module and how to resolve these issues?  To my
>> understanding, these issues continue not just in version 3 but
>> were not resolved in version 4 either.
>> 
>> Thanks,
>> Michael
>> 
>> 
>> 
>> 
>> On 09/19/10, Taha Al Hosary  <TAlHosary at ncms.ae> wrote:
>> 
>>> Dear all,
>>> When I ran the RegCM model I found that the results of precipitation 
>>> amount changes in whole by beginning of period. What is the most 
>>> appropriate start for good predict for the rainfall?
>>> Best Regards,
>>> Taha
>>> 
>>> 
>>> -----Original Message-----
>>> From: regcnet-bounces at lists.ictp.it 
>>> [mailto:regcnet-bounces at lists.ictp.it](javascript:main.compose() On Behalf 
>>> Of Azar Zarrin
>>> Sent: Tuesday, September 07, 2010 7:32 PM
>>> To: Regcnet at lists.ictp.it
>>> Subject: [RegCNET] Bug in the interactive lake model
>>> 
>>> Dear RegCM developers,
>>> 
>>> I'm running RegCM model coupled with Hostetler lake model for a cold 
>>> time period that ice and snow would be present over the lakes. I'm 
>>> getting the reasonable values for ice, temperature and evaporation 
>>> over the lakes, but the snow values are always zero.
>>> I've checked the lake.F and ice.F routines. The precip values which I 
>>> have printed from Ice. F are zero, while there are big amount of 
>>> precip and hence snow values outputted from the atmospheric model in 
>>> the SRF field. It seems that the precip isn't being passed into the 
>>> ice routine from the BATS and atmospheric model.
>>> 
>>> Would you please verify this bug and let me know how to solve this 
>>> error.
>>> 
>>> Many thanks,
>>> Azar
>>> 
>>> 
>>> ----
>>> Azar Zarrin
>>> Postdoctoral Research Associate
>>> Center for Climatic Research
>>> University of Wisconsin-Madison
>>> 1225 West Dayton Street, Room 1152
>>> Madison, Wisconsin 53706
>>> _______________________________________________
>>> RegCNET mailing list
>>> RegCNET at lists.ictp.it
>>> https://lists.ictp.it/mailman/listinfo/regcnet
>>> 
>>> _______________________________________________
>>> RegCNET mailing list
>>> RegCNET at lists.ictp.it
>>> https://lists.ictp.it/mailman/listinfo/regcnet
>> _______________________________________________
>> RegCNET mailing list
>> RegCNET at lists.ictp.it
>> https://lists.ictp.it/mailman/listinfo/regcnet
>> 
>
> ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
>  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
> ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   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
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~