Changeset 2443 for trunk/LMDZ.MARS/util/simu_MCS.F90

Timestamp:

Jan 15, 2021, 5:41:46 PM (4 years ago)

Author:

abierjon

Message:

Mars GCM:
Changes in aeroptical.F90 :

• add the flag nf90_64bit_offset when creating the output file. Without it, the program couldn't deal with large files (concat) with multiple aerosols
• modify the code structure to make it a bit less memory-intensive
• add topdust in computed aerosol opacities
• add the possibility to compute absorption opacities, instead of extinction : the user must specify 'ext' or 'abs' in aeroptical.def
• update the description of aeroptical in util/README

+ adapt simu_MCS to take into account aeroptical output variables when choosing

the binning reference

+ add outputs for topdust in physiq_mod

AB

File:

• trunk/LMDZ.MARS/util/simu_MCS.F90 (modified) (1 diff)

trunk/LMDZ.MARS/util/simu_MCS.F90

@@ -1206 +1206 @@
 !===============================================================================
     ! Observer variables to extract :
-    IF ((trim(GCMvarname).eq."dsodust").or.(trim(GCMvarname).eq."dso") &
-      .or.(trim(GCMvarname).eq."dustq").or.(trim(GCMvarname).eq."dust") &
-      .or.(trim(GCMvarname).eq."reffdust")) THEN
+    IF ((index(GCMvarname,"dust").ne.0).or.(index(GCMvarname,"dso").ne.0)) THEN
+      ! if the variable name contains "dust" or "dso". Especially for the targeted variables :
+      ! dustq,dustN,dsodust,reffdust,opadust, and their equivalents for stormdust & topdust
       OBSvarname  = "dust"
       numbin_name = "numbindust"
     ELSE IF ((trim(GCMvarname).eq."h2o_ice").or.(trim(GCMvarname).eq."wice") &
-           .or.(trim(GCMvarname).eq."reffice")) THEN
+           .or.(trim(GCMvarname).eq."reffice").or.(trim(GCMvarname).eq."opawice")) THEN
       OBSvarname  = "wice"
       numbin_name = "numbinwice"