Changeset 4009

Timestamp:

Jan 14, 2026, 8:30:25 PM (10 days ago)

Author:

aslmd

Message:

MESOSCALE: setting off subgrid slopes for mesoscale. still has to add a dimension with index 1 to many surface fields.

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• phys_state_var_init_mod.F90 (modified) (2 diffs)
• physiq_mod.F (modified) (4 diffs)

trunk/LMDZ.MARS/libf/phymars/phys_state_var_init_mod.F90

@@ -88 +88 @@
 #ifndef MESOSCALE
       if(is_master) then
-#endif
         status = nf90_open(filename, nf90_nowrite, ncid)
         if (status /= nf90_noerr) then
@@ -105 +104 @@
           endif
         endif
-#ifndef MESOSCALE
       endif
       call bcast(nslope)
+#else
+      print *, "subgrid slopes temporarily switched off in mesoscale, nslope=1"
+      nslope=1
 #endif
 

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

@@ -3172 +3172 @@
       comm_HR_LW(1:ngrid,1:nlayer) = zdtlw(1:ngrid,1:nlayer)
       !state  real  SWDOWNZ    ij   misc  1  -  h  "SWDOWNZ"   "DOWNWARD SW FLUX AT SURFACE"     "W m-2"
-      comm_SWDOWNZ(1:ngrid) = fluxsurf_dn_sw_tot(1:ngrid)
+      comm_SWDOWNZ(1:ngrid) = fluxsurf_dn_sw_tot(1:ngrid,1)
       !state  real  TAU_DUST   ij   misc  1  -  h  "TAU_DUST"  "REFERENCE VISIBLE DUST OPACITY"  ""
       comm_TAU_DUST(1:ngrid) = tau_pref_gcm(1:ngrid)
@@ -3179 +3179 @@
       !state  real  QSURFDUST  ij   misc  1  -  h  "QSURFDUST" "DUST MASS AT SURFACE"            "kg m-2"
       IF (igcm_dust_mass .ne. 0) THEN
-        comm_QSURFDUST(1:ngrid) = qsurf(1:ngrid,igcm_dust_mass)
+        comm_QSURFDUST(1:ngrid) = qsurf(1:ngrid,igcm_dust_mass,1)
       ELSE
         comm_QSURFDUST(1:ngrid) = 0.
@@ -3208 +3208 @@
       !! -- vdifc must have been called
       if (.not.callrichsl) then
-        sensibFlux(1:ngrid) = zflubid(1:ngrid)
-     .         - capcal(1:ngrid)*zdtsdif(1:ngrid)
+        sensibFlux(1:ngrid) = zflubid(1:ngrid,1)
+     .         - capcal(1:ngrid,1)*zdtsdif(1:ngrid,1)
       else
         sensibFlux(1:ngrid) =
@@ -3215 +3215 @@
      &   *sqrt(pu(1:ngrid,1)*pu(1:ngrid,1)+pv(1:ngrid,1)*pv(1:ngrid,1)
      &         +(log(1.+0.7*wstar(1:ngrid) + 2.3*wstar(1:ngrid)**2))**2)
-     &   *zcdh(1:ngrid)*(tsurf(1:ngrid)-zh(1:ngrid,1))
+     &   *zcdh(1:ngrid,1)*(tsurf(1:ngrid,1)-zh(1:ngrid,1))
       endif