Changeset 3483 for trunk/LMDZ.PLUTO/libf/phypluto/surfprop.F90

Timestamp:

Oct 24, 2024, 1:10:51 PM (4 weeks ago)

Author:

afalco

Message:

Pluto PCM: Added thermal inertia input when startphy_file == false.
Renamed thermal inertia variables to therm_inertia.
AF

File:

• trunk/LMDZ.PLUTO/libf/phypluto/surfprop.F90 (modified) (6 diffs)

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

@@ -1 @@
-      subroutine surfprop(ngrid,nq,pfra,qsurf,tsurface,     &
-                     tidat,capcal,adjust,dist,albedo,emis,fluold, &
-                     ptimestep,zls)
+      subroutine surfprop(ngrid,nq,pfra,qsurf,tsurface,       &
+                     capcal,adjust,dist,fluold,ptimestep,zls, &
+                     albedo,emis,therm_inertia)
 
     !   use comgeomfi_h, only:
@@ -43 +43 @@
 !     albedo(ngrid)
 !     emis(ngrid)
+!     therm_inertia(ngrid,nsoilmx)
 !
 !     Both
@@ -56 +57 @@
 !     Arguments
 
+
       INTEGER ngrid, nq
+      REAL,INTENT(IN) :: pfra(ngrid)
       REAL,INTENT(IN) :: qsurf(ngrid,nq)
+      REAL,INTENT(IN) :: tsurface(ngrid)
       REAL,INTENT(IN) :: fluold(ngrid,nq)
       REAL,INTENT(IN) :: ptimestep
       REAL,INTENT(IN) :: zls
-      REAL,INTENT(IN) :: tsurface(ngrid)
       REAL,INTENT(IN) :: capcal(ngrid)
       REAL,INTENT(IN) :: adjust
@@ -67 +70 @@
       REAL,INTENT(OUT) :: albedo(ngrid,L_NSPECTV)
       REAL,INTENT(OUT) :: emis(ngrid)
-      REAL,INTENT(OUT) :: tidat(ngrid,nsoilmx)
-      REAL,INTENT(IN) :: pfra(ngrid)
+      REAL,INTENT(OUT) :: therm_inertia(ngrid,nsoilmx) ! therm_inertia = inertiedat
 !-----------------------------------------------------------------------
 !     Local variables
@@ -448 +450 @@
            do isoil=0,nsoilmx-1
               if (mlayer(isoil).le.emin) then ! diurnal TI
-                   tidat(ig,isoil+1)=tid
+                   therm_inertia(ig,isoil+1)=tid
               else if (isoil.gt.0.and.(mlayer(isoil).gt.emin).and.(mlayer(isoil-1).lt.emin)) then ! still diurnal TI
-                   tidat(ig,isoil+1)=tid
+                   therm_inertia(ig,isoil+1)=tid
               else if ((mlayer(isoil).gt.emin).and.(mlayer(isoil).le.emax)) then ! TI N2
-                   tidat(ig,isoil+1)=ITN2
+                   therm_inertia(ig,isoil+1)=ITN2
               else if ((mlayer(isoil).gt.emin).and.(mlayer(isoil).le.emax2)) then
-                   tidat(ig,isoil+1)=ITCH4
+                   therm_inertia(ig,isoil+1)=ITCH4
               else
-                   tidat(ig,isoil+1)=ITH2O
+                   therm_inertia(ig,isoil+1)=ITH2O
               endif
 
@@ -465 +467 @@
 
         DO ig=1,ngrid
-           tidat(ig,:)=inertiedat(ig,:)
+           therm_inertia(ig,:)=inertiedat(ig,:)
         ENDDO