Changeset 2908 for trunk/LMDZ.MARS

Timestamp:

Feb 28, 2023, 2:48:15 PM (21 months ago)

Author:

romain.vande

Message:

Mars PCM:
Further implementation of subslope parametrisation.
Carefull ! This is a devolpment revision and it still need improvements and tests.
Only nslope=1 is possible for now!

!! This commit can produce small differences in restartfi.nc due to numerical roundup even for the case nslope=1 !!

Changes for tsurf:
tsurf_mesh_avg is used for callradite, thermosphere and in a specific case with Richardson based surface layer.

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/physiq_mod.F (modified) (4 diffs)

trunk/LMDZ.MARS/README

@@ -3923 +3923 @@
 Other code changes for nslope.gt.1 (sometimes the grid mesh average is used instead of the value of the subslope)
 
-
+== 28/02/2023 == RV
+This commit can produce small differences in restartfi.nc due to numerical roundup even for the case nslope=1.
+tsurf_mesh_avg is used for callradite, thermosphere and in a specific case with Richardson based surface layer.

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

@@ -1033 +1033 @@
            CALL callradite(icount,ngrid,nlayer,nq,zday,zls,pq,
      &     albedo_meshavg,emis_meshavg,
-     &     mu0,zplev,zplay,pt,tsurf(:,1),fract,dist_sol,igout,
+     &     mu0,zplev,zplay,pt,tsurf_meshavg,fract,dist_sol,igout,
      &     zdtlw,zdtsw,fluxsurf_lw(:,iflat),fluxsurf_dn_sw(:,:,iflat),
      &     fluxsurf_up_sw,
@@ -1057 +1057 @@
                CALL callradite(icount,ngrid,nlayer,nq,zday,zls,pq,
      &              albedo_meshavg,emis_meshavg,mu0,zplev,zplay,pt,
-     &              tsurf(:,1),fract,
+     &              tsurf_meshavg,fract,
      &              dist_sol,igout,zdtlwclf,zdtswclf,
      &              fluxsurf_lwclf,fluxsurf_dn_swclf,fluxsurf_up_swclf,
@@ -1467 +1467 @@
 
           DO ig=1, ngrid
-             IF (zh(ig,1) .lt. tsurf(ig,1)) THEN
+             IF (zh(ig,1) .lt. tsurf_meshavg(ig)) THEN
                wstar(ig)=1.
                hfmax_th(ig)=0.2
@@ -2146 +2146 @@
       if (callthermos) then
         call thermosphere(ngrid,nlayer,nq,zplev,zplay,dist_sol,
-     $     mu0,ptimestep,ptime,zday,tsurf(:,iflat),zzlev,zzlay,
+     $     mu0,ptimestep,ptime,zday,tsurf_meshavg,zzlev,zzlay,
      &     pt,pq,pu,pv,pdt,pdq,
      $     zdteuv,zdtconduc,zdumolvis,zdvmolvis,zdqmoldiff,