Changeset 726 for trunk/LMDZ.GENERIC/libf/phystd/largescale.F

Timestamp:

Jul 17, 2012, 5:04:27 PM (12 years ago)

Author:

jleconte

Message:

17/07/2012 == JL for LK

• Generalization of aerosol scheme:
  • any number of aerosols can be used and id numbers are determined consistently by the code. Aerosol order not important anymore.
  • addition of a module with the id numbers for aerosols (aerosol_mod.F90).
  • initialization of aerosols id numbers in iniaerosol.F90
  • compile with -s x where x *must* be equal to the number of aerosols turned on in callphys.def (either by a flag or by dusttau>0 for dust). => may have to erase object files when compiling with s option for the first time.
• For no aerosols, run with aeroco2=.true. and aerofixco2=.true (the default distribution for fixed co2

aerosols is 1.e-9; can be changed in aeropacity).

• If starting from an old start file, recreate start file with the q=0 option in newstart.e.
• update callphys.def with aeroXXX and aerofixXXX options (only XXX=co2,h2o supported for

now). Dust is activated by setting dusttau>0. See the early mars case in deftank.

• To add other aerosols, see Laura Kerber.

File:

• trunk/LMDZ.GENERIC/libf/phystd/largescale.F (modified) (1 diff)

trunk/LMDZ.GENERIC/libf/phystd/largescale.F

@@ -111 +111 @@
 
 !     for varying particle size in rad tran and (possibly) sedimentation
-         if(.not.aerofixed)then
+         if(aeroh2o.and.(.not.aerofixh2o)) then
+
             reffH2O(i,k) = CBRT( 3*zcond(i)/( 4*Nmix_h2o*pi*rho_ice))
             reffH2O(i,k) = max(reffH2O(i,k),1.e-8)