Changeset 3916 for trunk/LMDZ.PLUTO/libf

Timestamp:

Sep 18, 2025, 10:13:14 AM (3 months ago)

Author:

tbertrand

Message:

Pluto PCM :
Small fix to initialize n2frac in startfi
TB

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• dyn1d/rcm1d.F (modified) (3 diffs)
• phyetat0_mod.F90 (modified) (1 diff)

trunk/LMDZ.PLUTO/libf/phypluto/dyn1d/rcm1d.F

@@ -10 +10 @@
      &                     zmea, zstd, zsig, zgam, zthe,
      &                     emissiv, emisice, iceradius,
-     &                     dtemisice
+     &                     dtemisice, n2frac
       use comdiurn_h, only: sinlat, coslat, sinlon, coslon
       use comsoil_h, only: nsoilmx, layer, mlayer, inertiedat, volcapa
@@ -851 +851 @@
       write(*,*) " inertia = ",inertiedat(1,1)
 
+c  Initialisation n2frac
+c  --------------------------------------
+      n2frac(1)=1. ! default value for n2frac
+
 ! Initialize soil properties and temperature
 ! ------------------------------------------
@@ -920 +924 @@
       call physdem1("startfi.nc",nsoilmx,1,llm,nq,
      &                dtphys,time,
-     &                tsurf,tsoil,inertiedat,emis,albedodat,q2,qsurf)
+     &                tsurf,tsoil,inertiedat,emis,albedodat,q2,qsurf,
+     &                n2frac)
 
 c=======================================================================

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

@@ -148 +148 @@
   call get_field(nid_start,"n2frac",n2frac,found)
   if (.not.found) then
-    call abort_physic(modname,"Failed loading <n2frac>",1)
+    n2frac(:)=1.
+    !call abort_physic(modname,"Failed loading <n2frac>",1)
   endif
 else
   ! If no startfi file, use fixed value to test
-  n2frac(:)=1
+  n2frac(:)=1.
 endif ! of if (startphy_file)
 if (is_master) write(*,*) "phyetat0: n2 fraction of ice  <n2frac> range:", &