Changeset 706 for trunk/LMDZ.MARS/libf


Ignore:
Timestamp:
Jun 15, 2012, 10:55:19 AM (13 years ago)
Author:
emillour
Message:

Mars GCM:

  • Minor fix in "nuclea.F", max() and min() functions must have arguments of identical types.
  • Bug correction in "physiq.F": only water (and possibly dust and cnn if using microphysics) tracer tendencies should be updated after call to watercloud.

EM

Location:
trunk/LMDZ.MARS/libf/phymars
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.MARS/libf/phymars/nuclea.F

    r695 r706  
    115115          deltaf = (2.*desorp-surfdif-fistar)/
    116116     &             (kbz*temp)
    117           deltaf = min( max(deltaf, -100.), 100.)
     117          deltaf = min( max(deltaf, -100.d0), 100.d0)
    118118
    119119          if (deltaf.eq.-100.) then

  • trunk/LMDZ.MARS/libf/phymars/physiq.F

    r705 r706  
    10041004           
    10051005! increment water vapour and ice atmospheric tracers tendencies
    1006            pdq(1:ngrid,1:nlayer,1:nq) =
    1007      &       pdq(1:ngrid,1:nlayer,1:nq) +
    1008      &       zdqcloud(1:ngrid,1:nlayer,1:nq)
    1009            
     1006           pdq(1:ngrid,1:nlayer,igcm_h2o_vap) =
     1007     &       pdq(1:ngrid,1:nlayer,igcm_h2o_vap) +
     1008     &       zdqcloud(1:ngrid,1:nlayer,igcm_h2o_vap)
     1009           pdq(1:ngrid,1:nlayer,igcm_h2o_ice) =
     1010     &       pdq(1:ngrid,1:nlayer,igcm_h2o_ice) +
     1011     &       zdqcloud(1:ngrid,1:nlayer,igcm_h2o_ice)
     1012
     1013! increment dust and ccn masses and numbers
     1014           if (microphys) then
     1015             pdq(1:ngrid,1:nlayer,igcm_dust_mass) =
     1016     &         pdq(1:ngrid,1:nlayer,igcm_dust_mass) +
     1017     &         zdqcloud(1:ngrid,1:nlayer,igcm_dust_mass)
     1018             pdq(1:ngrid,1:nlayer,igcm_dust_number) =
     1019     &         pdq(1:ngrid,1:nlayer,igcm_dust_number) +
     1020     &         zdqcloud(1:ngrid,1:nlayer,igcm_dust_number)
     1021             pdq(1:ngrid,1:nlayer,igcm_ccn_mass) =
     1022     &         pdq(1:ngrid,1:nlayer,igcm_ccn_mass) +
     1023     &         zdqcloud(1:ngrid,1:nlayer,igcm_ccn_mass)
     1024             pdq(1:ngrid,1:nlayer,igcm_ccn_number) =
     1025     &         pdq(1:ngrid,1:nlayer,igcm_ccn_number) +
     1026     &         zdqcloud(1:ngrid,1:nlayer,igcm_ccn_number)
     1027           endif
     1028
    10101029! We need to check that we have Nccn & Ndust > 0
    10111030! This is due to single precision rounding problems
Note: See TracChangeset for help on using the changeset viewer.