Changeset 2443 for trunk/LMDZ.MARS/libf/phymars

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/libf/phymars/physiq_mod.F (modified) (2 diffs)

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -3167 +3167 @@
 
            if (rdstorm) then  ! writediagfi tendencies stormdust tracers
-             call WRITEDIAGFI(ngrid,'reffdust','reffdust',
-     &                        'm',3,rdust*ref_r0)
              call WRITEDIAGFI(ngrid,'reffstormdust','reffstormdust',
      &                        'm',3,rstormdust*ref_r0)
@@ -3226 +3224 @@
 
            if (slpwind) then
+             call WRITEDIAGFI(ngrid,'refftopdust','refftopdust',
+     &                        'm',3,rtopdust*ref_r0)
+             call WRITEDIAGFI(ngrid,'topdustq','top Dust mass mr',
+     &                        'kg/kg',3,pq(:,:,igcm_topdust_mass))
+             call WRITEDIAGFI(ngrid,'topdustN','top Dust number',
+     &                        'part/kg',3,pq(:,:,igcm_topdust_number))
              select case (trim(dustiropacity))
               case ("tes")