Changeset 1604 for LMDZ5


Ignore:
Timestamp:
Dec 22, 2011, 12:26:02 PM (13 years ago)
Author:
lguez
Message:

Removed two "fi" with no corresponding "if" in "makdim".

In procedure "coefkz", when calculating cloud fraction "zfr", "zq" can
be zero at high altitude. Avoid division by zero, set "zfr" to 0 when
"zq" is 0.

In "physiq", allocatable arrays "tabijgcm", "longcm", "latgcm",
"igcm", "jgcm" cannot be arguments of "phys_output_open" if they have
not been allocated. Allocate them with zero size when
'npCFMIP_param.data' cannot be opened.

Location:
LMDZ5/trunk/libf
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/libf/grid/dimension/makdim

    r1146 r1604  
    1 for i in $* ; do
    2    list=$list.$i
     1for i in $*
     2  do
     3  list=$list.$i
    34done
    45fichdim=dimensions${list}
    56
    6 if [ ! -f $fichdim ] ; then
     7if [ ! -f $fichdim ]
     8    then
    79# si le fichier de dimensions n'existe pas, on le cree
    810
    9   if [ $# -ge 3 ] ; then
    10      im=$1
    11      jm=$2
    12      lm=$3
    13      n2=$1
    14      ndm=1
     11    if [ $# -ge 3 ]
     12        then
     13        im=$1
     14        jm=$2
     15        lm=$3
     16        n2=$1
     17        ndm=1
    1518
    1619# Le test suivant est commente car il est inutile avec le nouveau
    … …  
    2831
    2932
    30   else if [ $# -ge 2 ] ; then
    31       im=1
    32       jm=$1
    33       lm=$2
    34       ndm=1
    35   else if [ $# -ge 1 ] ; then
    36       im=1
    37       jm=1
    38       lm=$1
    39       ndm=1
    40   else
    41          echo il faut au moins une dimension
    42          exit
    43   fi
    44 fi
     33    else if [ $# -ge 2 ]
     34        then
     35        im=1
     36        jm=$1
     37        lm=$2
     38        ndm=1
     39    else if [ $# -ge 1 ]
     40        then
     41        im=1
     42        jm=1
     43        lm=$1
     44        ndm=1
     45    else
     46        echo il faut au moins une dimension
     47        exit
     48    fi
    4549fi
    4650
    … …  
    6064EOF
    6165
    62 fi
    63 
    6466\rm ../dimensions.h
    6567tar cf - $fichdim | ( cd .. ; tar xf - ; mv $fichdim dimensions.h )

  • LMDZ5/trunk/libf/phylmd/coef_diff_turb_mod.F90

    r1067 r1604  
    389389!           calculer la fraction nuageuse (processus humide):
    390390!
    391           zfr = (zq+ratqs*zq-zqs) / (2.0*ratqs*zq)
     391          if (zq /= 0.) then
     392             zfr = (zq+ratqs*zq-zqs) / (2.0*ratqs*zq)
     393          else
     394             zfr = 0.
     395          end if
    392396          zfr = MAX(0.0,MIN(1.0,zfr))
    393397          IF (.NOT.richum) zfr = 0.0

  • LMDZ5/trunk/libf/phylmd/physiq.F

    r1565 r1604  
    12201220      INTEGER :: nbtr_tmp ! Number of tracer inside concvl
    12211221      REAL, dimension(klon,klev) :: sh_in ! Specific humidity entering in phytrac
     1222      integer iostat
    12221223
    12231224cIM for NMC files
    … …  
    15081509      nCFMIP=npCFMIP
    15091510      OPEN(98,file='npCFMIP_param.data',status='old',
    1510      $          form='formatted',err=999)
     1511     $          form='formatted',iostat=iostat)
     1512            if (iostat == 0) then
    15111513      READ(98,*,end=998) nCFMIP
    15121514998   CONTINUE
    … …  
    15401542     $tabijGCM, lonGCM, latGCM, iGCM, jGCM)
    15411543c
    1542 999      CONTINUE
     1544            else
     1545               ALLOCATE(tabijGCM(0))
     1546               ALLOCATE(lonGCM(0), latGCM(0))
     1547               ALLOCATE(iGCM(0), jGCM(0))
     1548            end if
    15431549         ENDIF !debut
    15441550 
Note: See TracChangeset for help on using the changeset viewer.