[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