[Ncep.list.emc.monsoondesk] Bug which caused CFSv2 giving different results when layout or/and no procs for ocean model was changed.
gibies george
gibies at tropmet.res.in
Sat Feb 15 05:50:09 UTC 2014
Prajeesh,
I think you are right, and that is a very serious bug in the ICE model for
CFSv2. But my suggestion is a slightly different solution. It can be done
by properly declaring the input variable in subroutine
"ice_top_to_ice_bottom".
subroutine ice_top_to_ice_bottom (lcplocn, gfs_flx, coldstart_daily_flx,
Ice, part_size, part_size_uv)
type (ice_data_type), intent(inout) :: Ice
real, dimension (isc:iec,jsc:jec,km), intent(in) :: part_size,
part_size_uv
..................................
.............................
................
_____________________________________________________
I think, it is possible because, the domain boundaries are
globally declared in the module as follows.
use ice_grid_mod, only: Domain, isc, iec, jsc, jec, isd, ied, jsd,
jed, im, jm, km
*Gibies George, *
*CSIR-Senior Research Fellow,*
*Climate and Global Modelling Division,*
*Indian Institute of Tropical Meteorology,*
*Dr. Homi Bhabha Road,*
*NCL (P. O.), Pashan,*
*Pune 411008, India.*
*Office Teleph : +912025904423*
*Office mail : gibies at tropmet.res.in <gibies at tropmet.res.in>*
*________________*
*Vettikuzhiyil House*
*Ezhakkaranadu South*
*Puthencruz-682308*
*Kochi*
*Mobile number : **+919028523910,* *+919420407392,** +919881707796*
*Personal mail : gibiesgeorge at gmail.com <gibiesgeorge at gmail.com>*
*http://sites.google.com/site/gibiesge/
<http://sites.google.com/site/gibiesge/>_______________________________________________________*
*"Please do not print this e-mail unless you really need. Because e**very
3000 pages cost a tree. Save Trees, Save mother Earth." *
On 14 February 2014 17:02, Prajeesh A.G CAT <prajeesh.cat at tropmet.res.in>wrote:
>
> In the subroutine *ice_top_to_ice_bottom* in the source file
> *ice_model.f90* of mom4, for coupling with gfs the following codes where
> added:
>
>
>
>
>
>
>
>
>
>
>
>
>
>
> *if (gfs_flx .AND. (.NOT. coldstart_daily_flx) ) then do j=jsc,jec
> do i=isc,iec if (part_size (i,j,1) .GT. 0.99) then Ice %
> flux_u (i,j) = Ice % flux_au (i,j) Ice % flux_v (i,j) = Ice %
> flux_av (i,j) Ice % flux_t (i,j) = Ice % flux_at (i,j)
> Ice % flux_q (i,j) = Ice % flux_aq (i,j) Ice % flux_sw(i,j) = Ice
> % flux_asw(i,j) Ice % flux_lw(i,j) = Ice % flux_alw(i,j)
> endif enddo enddoendif*
>
> The if statement in the above line is for checking if it is Ice free or
> not, (part_size > 0.99 means Ice free).
>
> But the lower bound (lbound) and upper bound (ubound) of part_size(:,:,1)
> and Ice% Variables (Ice%flux_au(:,:), Ice%flux_au(:,:), ......) *are
> different.*
>
> The following is an values of lbound and ubound of part_size and Ice%flux_
> au from out.poe:
>
> 8: lbound(part_size) 1 1
> 8: ubound(part_size) 12 410
> 8: lbound(Ice%flux_au) 85 1
> 8: ubound(Ice%flux_au) 96 410
> 9: lbound(part_size) 1 1
> 9: ubound(part_size) 12 410
> 9: lbound(Ice%flux_au) 97 1
> 9: ubound(Ice%flux_au) 108 410
> 10: lbound(part_size) 1 1
> 10: ubound(part_size) 12 410
> 10: lbound(Ice%flux_au) 109 1
> 10: ubound(Ice%flux_au) 120 410
>
> If then, the whole above code will go wrong.
>
> And this was the reason for getting different results when changing the
> no: of procs or layout of ocean model in CFSv2.
>
> A simple solution for this could be by modifying the code to:
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
> *if (gfs_flx .AND. (.NOT. coldstart_daily_flx) ) then jj=0 do
> j=jsc,jec jj=jj+1 ii=0 do i=isc,iec ii=ii+1 if
> (part_size (ii,jj,1) .GT. 0.99) then Ice % flux_u (i,j) = Ice %
> flux_au (i,j) Ice % flux_v (i,j) = Ice % flux_av (i,j)
> Ice % flux_t (i,j) = Ice % flux_at (i,j) Ice % flux_q (i,j) = Ice
> % flux_aq (i,j) Ice % flux_sw(i,j) = Ice % flux_asw(i,j)
> Ice % flux_lw(i,j) = Ice % flux_alw(i,j) endif enddo enddoendif*Any
> comments on this are appreciated.
>
> Prajeesh A G
> Scientist 'B'
> CCCR, IITM
> Pune, India
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: https://lstsrv.ncep.noaa.gov/pipermail/ncep.list.emc.monsoondesk/attachments/20140215/c717b885/attachment-0001.html
More information about the Ncep.list.emc.monsoondesk
mailing list