Changeset 3910 for trunk/LMDZ.PLUTO/libf/phypluto

Timestamp:

Sep 3, 2025, 11:19:36 PM (6 months ago)

Author:

tbertrand

Message:

Pluto PCM:

• adding option to use N2 ice fractional maps (n2frac) read in the startfi.nc
• adding option in newstart to correct (tsurf and tsoil) for too warm or too cold N2-free (correct_t_non2) or N2-rich (correct_t_n2) patches

TB

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• condense_n2.F90 (modified) (3 diffs)
• phyetat0_mod.F90 (modified) (3 diffs)
• phyredem.F90 (modified) (3 diffs)
• phys_state_var_mod.F90 (modified) (3 diffs)
• physiq_mod.F90 (modified) (2 diffs)
• surfdat_h.F90 (modified) (1 diff)

trunk/LMDZ.PLUTO/libf/phypluto/condense_n2.F90

@@ -4 +4 @@
       picen2,psolaralb,pemisurf, &
       pdtc,pdtsrfc,pdpsrf,pduc,pdvc, &
-      pdqc,pdicen2)
+      pdqc,pdicen2,n2frac)
 
   use comgeomfi_h
@@ -71 +71 @@
   REAL pdtsrfc(klon),pdpsrf(klon)
   REAL picen2(klon),psolaralb(klon),pemisurf(klon)
+  REAL n2frac(klon)
 
 
@@ -370 +371 @@
 
          pdtsrfc(ig) = (ztcondsol(ig) - ztsrf(ig))/subtimestep
+
+         zcondices(ig) = zcondices(ig)*n2frac(ig) !N2frac added as ponderation
 
 !        partial sublimation of N2 ice

trunk/LMDZ.PLUTO/libf/phypluto/phyetat0_mod.F90

@@ -23 +23 @@
   use tabfi_mod, only: tabfi,ini_tab_controle_dyn_xios,tab_cntrl_mod
   USE tracer_h, ONLY: noms
-  USE surfdat_h, only: phisfi, albedodat, zmea, zstd, zsig, zgam, zthe
+  USE surfdat_h, only: phisfi, albedodat, zmea, zstd, zsig, zgam, zthe, n2frac
   use iostart, only: nid_start, open_startphy, close_startphy, &
                      get_field, get_var, inquire_field, &
@@ -91 +91 @@
 !
 IF (.not. ALLOCATED(albedodat)) ALLOCATE(albedodat(ngrid))
+IF (.not. ALLOCATED(n2frac)) ALLOCATE(n2frac(ngrid))
 IF (.not. ALLOCATED(phisfi)) ALLOCATE(phisfi(ngrid))
 IF (.not. ALLOCATED(zmea)) ALLOCATE(zmea(ngrid))
@@ -140 +141 @@
 if (is_master) write(*,*) "phyetat0: Bare ground albedo <albedodat> range:", &
              minval(albedodat), maxval(albedodat)
+
+!n2frac
+
+if (startphy_file) then
+  ! Load bare ground albedo:
+  call get_field(nid_start,"n2frac",n2frac,found)
+  if (.not.found) then
+    call abort_physic(modname,"Failed loading <n2frac>",1)
+  endif
+else
+  ! If no startfi file, use fixed value to test
+  n2frac(:)=1
+endif ! of if (startphy_file)
+if (is_master) write(*,*) "phyetat0: n2 fraction of ice  <n2frac> range:", &
+             minval(n2frac), maxval(n2frac)
 
 ! ZMEA

trunk/LMDZ.PLUTO/libf/phypluto/phyredem.F90

@@ -173 +173 @@
 subroutine physdem1(filename,nsoil,ngrid,nlay,nq, &
                     phystep,time,tsurf,tsoil,inertiesoil, &
-                    emis,alb,q2,qsurf)
+                    emis,alb,q2,qsurf,n2frac)
   ! write time-dependent variable to restart file
   use iostart, only : open_restartphy, close_restartphy, &
@@ -192 +192 @@
   real,intent(in) :: emis(ngrid)
   real,intent(in) :: alb(ngrid)
+  real,intent(in) :: n2frac(ngrid)
   real,intent(in) :: inertiesoil(ngrid,nsoil)
   real,intent(in) :: q2(ngrid,nlay+1)
@@ -219 +220 @@
   ! Albedo of the surface
   call put_field(nid_restart,"albedodat","Albedo of bare ground",alb)
+
+  !n2 fraction on the surface
+  call put_field(nid_restart,"n2frac","N2 ice fraction on the surface",n2frac)
 
   ! Planetary Boundary Layer

trunk/LMDZ.PLUTO/libf/phypluto/phys_state_var_mod.F90

@@ -14 +14 @@
       use radcommon_h, only: glat
       USE radinc_h, only : L_NSPECTI, L_NSPECTV,naerkind
-      use surfdat_h, only: phisfi, albedodat,  &
+      use surfdat_h, only: phisfi, albedodat, n2frac, &
                         zmea, zstd, zsig, zgam, zthe
       use turb_mod, only: q2,sensibFlux,wstar,ustar,tstar,hfmax_th,zmax_th
@@ -127 +127 @@
         ALLOCATE(albedo_n2_ice_SPECTV(L_NSPECTV))
         ALLOCATE(albedo_bareground(klon))
+        ALLOCATE(n2frac(klon))
         ALLOCATE(emis(klon))
         ALLOCATE(dtrad(klon,klev))
@@ -251 +252 @@
         DEALLOCATE(phisfi)
         DEALLOCATE(albedodat)
+        DEALLOCATE(n2frac)
         DEALLOCATE(zmea)
         DEALLOCATE(zstd)

trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90

@@ -1397 +1397 @@
                            qsurf(1,igcm_n2),albedo,emis,                &
                            zdtc,zdtsurfc,pdpsrf,zduc,zdvc,              &
-                           zdqc,zdqsc(1,igcm_n2))
+                           zdqc,zdqsc(1,igcm_n2),n2frac)
 
          pdt(1:ngrid,1:nlayer) = pdt(1:ngrid,1:nlayer)+zdtc(1:ngrid,1:nlayer)
@@ -2140 +2140 @@
             call physdem1("restartfi.nc",nsoilmx,ngrid,nlayer,nq, &
                         ptimestep,ztime_restart,tsurf,                &
-                        tsoil,therm_inertia,emis,albedo,q2,qsurf)
+                        tsoil,therm_inertia,emis,albedo,q2,qsurf,n2frac)
          endif ! ngrid
       endif ! is_omp_master

trunk/LMDZ.PLUTO/libf/phypluto/surfdat_h.F90

@@ -3 +3 @@
 
        implicit none
+
+       real,allocatable,dimension(:) :: n2frac ! n2 fraction of ice at the surfacestocked in startfi.nc file.
+!$OMP THREADPRIVATE(n2frac)
 
        real,allocatable,dimension(:) :: albedodat ! albedo of bare ground stocked in startfi.nc file.