Changeset 654 for trunk/LMDZ.MARS/libf


Ignore:
Timestamp:
May 10, 2012, 10:02:41 PM (13 years ago)
Author:
tnavarro
Message:

bug correction if all dust disappears in cloud formation

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r645 r654  
    277277        DO l=1,nlay
    278278         DO ig=1,ngrid
    279           IF (pq(ig,l,igcm_ccn_number) +
     279          IF ((pq(ig,l,igcm_dust_number) +
     280     &      ptimestep* (pdq(ig,l,igcm_dust_number) +
     281     &        pdqcloud(ig,l,igcm_dust_number)) .le. 1.)
     282     &   .or. (pq(ig,l,igcm_dust_mass) +
     283     &      ptimestep* (pdq(ig,l,igcm_dust_mass) +
     284     &        pdqcloud(ig,l,igcm_dust_mass)) .le. 1.e-20)) THEN
     285         pdqcloud(ig,l,igcm_dust_number) =
     286     &     - pq(ig,l,igcm_dust_number)/ptimestep
     287     &     - pdq(ig,l,igcm_dust_number) + 1.
     288         pdqcloud(ig,l,igcm_ccn_number) = 
     289     &     -pdqcloud(ig,l,igcm_dust_number)
     290         pdqcloud(ig,l,igcm_dust_mass) =
     291     &     - pq(ig,l,igcm_dust_mass)/ptimestep
     292     &     - pdq(ig,l,igcm_dust_mass) + 1.e-20
     293         pdqcloud(ig,l,igcm_ccn_mass) =
     294     &     -pdqcloud(ig,l,igcm_dust_mass)
     295          ENDIF
     296          IF ((pq(ig,l,igcm_ccn_number) +
    280297     &      ptimestep* (pdq(ig,l,igcm_ccn_number) +
    281      &        pdqcloud(ig,l,igcm_ccn_number)) .le. 0.) THEN
     298     &        pdqcloud(ig,l,igcm_ccn_number)) .le. 1.)
     299     &   .or. (pq(ig,l,igcm_ccn_mass) +
     300     &      ptimestep* (pdq(ig,l,igcm_ccn_mass) +
     301     &        pdqcloud(ig,l,igcm_ccn_mass)) .le. 1.e-20)) THEN
    282302         pdqcloud(ig,l,igcm_ccn_number) =
    283303     &     - pq(ig,l,igcm_ccn_number)/ptimestep
    284      &     - pdq(ig,l,igcm_ccn_number)
     304     &     - pdq(ig,l,igcm_ccn_number) + 1.
    285305         pdqcloud(ig,l,igcm_dust_number) = 
    286306     &     -pdqcloud(ig,l,igcm_ccn_number)
    287307         pdqcloud(ig,l,igcm_ccn_mass) =
    288308     &     - pq(ig,l,igcm_ccn_mass)/ptimestep
    289      &     - pdq(ig,l,igcm_ccn_mass)
     309     &     - pdq(ig,l,igcm_ccn_mass) + 1.e-20
    290310         pdqcloud(ig,l,igcm_dust_mass) =
    291311     &     -pdqcloud(ig,l,igcm_ccn_mass)
Note: See TracChangeset for help on using the changeset viewer.