[RegCNET] Another Bug confirmed and fixed in RegCM4.2 with CLM option

[guangshan chen]

Dear Phan Van Tan,

This fixing is not to fix the problem of missing values or zero in SRF and SST.

Some variables in SRF, you can find them in clmoutput.clm2.h0*.nc, such as TSA (2 meters temperature).

I think it will take time to find the reason.

Sorry, I do not have clue now.

Guangshan

On May 12, 2012, at 11:15 AM, Phan Van Tan wrote:

> Dear Guangshan Chen,
> 
> After fixing, the model output in the SRF and STS files are still missing values or zero !!!???
> 
> 
> --- On Fri, 5/11/12, Guangshan Chen <gchen9 at gmail.com> wrote:
> 
> From: Guangshan Chen <gchen9 at gmail.com>
> Subject: Another Bug confirmed and fixed in RegCM4.2 with CLM option
> To: regcnet at lists.ictp.it
> Cc: "Phan Van Tan" <tanpv2000 at yahoo.com>, "chao sun" <1988sunchao at gmail.com>
> Date: Friday, May 11, 2012, 8:48 AM
> 
> Dear all,
> 
> There is another bug confirmed and fixed in RegCM4.2 with CLM option.
> 
> First, thanks for Chao Sun who found it and Phan Van Tan who provided part of fixing.
> 
> This bug is in ./Main/clmlib/clm_atmlnd.F90, line 1370 and line 1374.
> 
> I have confirmed that line 1370 and line 1374 do cause a segment fault due to reallocation of c2r_allout by each processor.
> (if you use one processor, it is no problem. Also you may can not find it by just run the model for couple of days. I find it after
> about one month. At first, I just tested the code for couple of days and did not find it.) After fixing it, the model passed the crashed point.
> It is still running now (two months already).
> 
> The lines 1370 and 1374 should be fixed as follow:
> 
>      if(AERTYP.ne.'AER00D0') then  !Aerosol scheme on
>        nout = 22
>        allocate(c2r_all(nt*nout))
>        if (.NOT. allocated(c2r_allout(numg*nout))) allocate(c2r_allout(numg*nout))
>      else
>        nout = 20
>        allocate(c2r_all(nt*nout))
>        if (.NOT. allocated(c2r_allout(numg*nout))) allocate(c2r_allout(numg*nout))
>      end if
> 
> I made a new patch file based on the old one in ./Main/clmlib/patchset. You can use this new patch to replace the old one.
> 
> 
> Next step, comment deallocate(c2r_allout) (delete this line) in ./Main/clmlib/mod_mtrxclm.F90.
> I had made a patch file for mod_mtrxclm.F90. The attached is a new one.
> 
> Guangshan Chen
> University of Wisconsin-Madison, USA
> 
> 
> 
> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.ictp.it/pipermail/regcnet/attachments/20120512/70912148/attachment.html>