Changeset 3125 for trunk/LMDZ.MARS/libf

Timestamp:

Nov 14, 2023, 3:25:37 PM (16 months ago)

Author:

abierjon

Message:

Mars GCM:
Add a sanity check on tau_pref_gcm in aeropacity when running without rescaling (dustscaling_mode=/=1), hence resolving ticket #156
+ add some comments and a sanity check in conf_phys to ensure freedust flag is not used with dustscaling_mode=1.

AB

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• aeropacity_mod.F (modified) (1 diff)
• conf_phys.F (modified) (3 diffs)

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

@@ -911 +911 @@
         ENDDO
         tau_pref_gcm(:) = tau_pref_gcm(:) * odpref / pplev(:,1)
+        
+        ! Sanity check : if tau_pref_gcm > 50 with no rescaling, there is probably an error somewhere
+        if (maxval(tau_pref_gcm(:)).gt.50) then
+          write(*,*) "Error in aeropacity: ",
+     &     "tau_pref_gcm>50 while there is no rescaling to ",
+     &     "the scenario (dustscaling_mode =/= 1). ",
+     &     "Better stop here."
+          write(*,*) "If you want to run with non-rescaled dust (GCM6)",
+     &    " be sure to have good dust fields in start.nc using newstart"
+          call abort_physic("aeropacity",
+     &     "tau_pref_gcm>50 while no rescaling",1)
+        endif
+        
       ELSE
         ! dust opacity strictly follows what is imposed by the dust scenario

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

@@ -448 +448 @@
 
 ! dust rescaling mode (if any)
+!     =0 --> freedust, tauscaling=1
+!     =1 --> GCM5.3-like, tauscaling computed wrt tau_pref_scenario
+!     =2 --> tauscaling=1 but dust_rad_adjust (rescaling for the radiative transfer only)
          if (freedust) then
            dustscaling_mode=0
@@ -455 +458 @@
          call getin_p("dustscaling_mode",dustscaling_mode)
          write(*,*) "dustscaling_mode=",dustscaling_mode
+         if ((dustscaling_mode.eq.1).and.freedust) then
+           print*,'freedust and dustscaling_mode=1 are incompatible !'
+           call abort_physic(modname,
+     &          "freedust and dustscaling_mode=1 are incompatible",1)
+         endif         
 
 #ifndef MESOSCALE
@@ -494 +502 @@
      &            should be used !'
            call abort_physic(modname,
-     &          "rdstorm or topflows requires activ",1)
+     &          "rdstorm or topflows requires active",1)
          endif
          if (rdstorm.and..not.lifting) then