Ignore:
Timestamp:
Jun 21, 2022, 11:05:46 AM (3 years ago)
Author:
aslmd
Message:

rainthreshold constant & delete useless lines

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.GENERIC/libf/phystd/rain_generic.F90

    r2722 r2723  
    8484
    8585   ! for simple scheme : precip_scheme_generic=1
    86    real,parameter :: t_crit=218.0 ! need to be modified to be generic for all species
    8786   real lconvert
    8887
     
    129128
    130129      if(call_ice_vap_generic) then ! to call only one time the ice/vap pair of a tracer
    131 
    132          ! Let's get the index of our tracers (we look for igcm _generic_vap and igcm_generic_ice)
    133          igcm_generic_vap=iq
    134          igcm_generic_ice = -1
    135 
    136          ! We look for the corresponding ice traceur (before or after in the list of traceurs, maybe could be generalised to the whole list)
    137          if (iq .ne. nq) then
    138             if ((noms(iq)(1:len(trim(noms(iq)))-4) .eq. noms(iq+1)(1:len(trim(noms(iq+1)))-4)) .and. (index(noms(iq+1),"ice") .ne. 0)) then
    139                igcm_generic_ice = iq+1
    140             end if
    141          end if
    142          if ((iq .gt. 1)) then
    143             if ((noms(iq)(1:len(trim(noms(iq)))-4) .eq. noms(iq-1)(1:len(trim(noms(iq-1)))-4)) .and. (index(noms(iq-1),"ice") .ne. 0)) then
    144                   igcm_generic_ice = iq-1
    145             end if
    146          end if
    147          if (igcm_generic_ice .eq. -1) then
    148             write(*,*) "ERROR : You set a vap traceur but you forgot to set the corresponding ice traceur, &
    149             or the pair vap/ice is not written one after another in traceur.def"
    150          endif
    151130     
    152131         m=constants_mass(iq)
     
    314293
    315294               do i = 1, ngrid
    316                   t_tmp = zt(i,k)
    317                   if (t_tmp .ge. T_h2O_ice_liq) then
    318                      lconvert=rainthreshold_generic
    319                   elseif (t_tmp .gt. t_crit) then
    320                      lconvert=rainthreshold_generic*(1.- t_crit/t_tmp)
    321                      lconvert=MAX(0.0,lconvert)             
    322                   else
    323                      lconvert=0.
    324                   endif
    325 
    326                   ! lconvert is between 0. and rainthreshold_generic
    327 
     295                 
     296                  lconvert=rainthreshold_generic
    328297
    329298                  if (ql(i,k).gt.1.e-9) then
Note: See TracChangeset for help on using the changeset viewer.