Changeset 3228 for trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto

Timestamp:

Feb 20, 2024, 4:46:35 PM (22 months ago)

Author:

tbertrand

Message:

Pluto GCM:
Cleaning code and adding extra options for newstart.F

File:

• trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto/newstart.F (modified) (3 diffs)

trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto/newstart.F

@@ -120 +120 @@
       real mugaz ! molar mass of the atmosphere
 
-      integer ierr 
+      integer ierr,iref 
 
 c Variables on the new grid along scalar points 
@@ -585 +585 @@
       write(*,*) 'qs=x : give a uniform value to a surface tracer'
       write(*,*) 'q=profile    : specify a profile for a tracer'
+      write(*,*) 'subsoil_all : set seasonal subsurface thermal inertia'
+      write(*,*) 'diurnal_TI : set diurnal subsurface thermal inertia'
 
         write(*,*)
@@ -824 +826 @@
 
 
+
+c       subsoil_all : initialize subsurface thermal inertia
+c       --------------------------------------------------
+        else if (modif(1:len_trim(modif)) .eq. 'subsoil_all') then
+
+          write(*,*) 'New value for subsoil thermal inertia  ?'
+ 104      read(*,*,iostat=ierr) ith_bb
+          if(ierr.ne.0) goto 104
+          write(*,*) 'thermal inertia (new value):',ith_bb
+
+          write(*,*)'At which depth (in m.) does the ice layer begin?'
+          write(*,*)'(here, the deepest soil layer extends down to:'
+     &              ,layer(1),' - ',layer(nsoilmx),')'
+          write(*,*)'write 0 for uniform value for all subsurf levels?'
+          ierr=1
+          do while (ierr.ne.0)
+           read(*,*,iostat=ierr) val2
+           write(*,*)'val2 in subsoil_all:',val2,'ierr=',ierr
+           if (ierr.eq.0) then ! got a value, but do a sanity check
+             if(val2.gt.layer(nsoilmx)) then
+               write(*,*)'Depth should be less than ',layer(nsoilmx)
+               ierr=1
+             endif
+             if(val2.lt.layer(1)) then
+              if(val2.eq.0) then
+               write(*,*)'Thermal inertia set for all subsurface layers'
+               ierr=0
+              else
+               write(*,*)'Depth should be more than ',layer(1)
+               ierr=1
+              endif
+             endif
+           endif
+          enddo ! of do while
+
+          ! find the reference index iref the depth corresponds to
+          if(val2.eq.0) then
+           iref=1
+           write(*,*)'Level selected is first level: ',layer(iref),' m'
+          else
+           do isoil=1,nsoilmx-1
+            if ((val2.gt.layer(isoil)).and.(val2.lt.layer(isoil+1)))
+     &       then
+             iref=isoil+1
+             write(*,*)'Level selected : ',layer(isoil+1),' m'
+             exit
+            endif
+           enddo
+          endif
+
+          DO ig=1,ngridmx
+             DO j=iref,nsoilmx
+                   ithfi(ig,j)=ith_bb
+             ENDDO
+          ENDDO
+
+          CALL gr_fi_dyn(nsoilmx,ngridmx,iip1,jjp1,ithfi,ith)
+
+c       diurnal_TI : choice of thermal inertia values (global)
+c       ----------------------------------------------------------------
+        else if (modif(1:len_trim(modif)) .eq. 'diurnal_TI') then
+
+         write(*,*) 'New value for diurnal thermal inertia  ?'
+ 106     read(*,*,iostat=ierr) ith_bb
+         if(ierr.ne.0) goto 106
+         write(*,*) 'Diurnal thermal inertia (new value):',ith_bb
+
+         write(*,*)'At which depth (in m.) does the ice layer ends?'
+         write(*,*)'(currently, the soil layer 1 and nsoil are:'
+     &              ,layer(1),' - ',layer(nsoilmx),')'
+         ierr=1
+         do while (ierr.ne.0)
+          read(*,*,iostat=ierr) val2
+          write(*,*)'val2 in diurnal_TI:',val2,'ierr=',ierr
+          if (ierr.eq.0) then ! got a value, but do a sanity check
+            if(val2.gt.layer(nsoilmx)) then
+              write(*,*)'Depth should be less than ',layer(nsoilmx)
+              ierr=1
+            endif
+            if(val2.lt.layer(1)) then
+              write(*,*)'Depth should be more than ',layer(1)
+              ierr=1
+            endif
+          endif
+         enddo ! of do while
+
+           ! find the reference index iref the depth corresponds to
+         do isoil=1,nsoilmx-1
+            !write(*,*)'isoil, ',isoil,val2
+            !write(*,*)'lay(i),lay(i+1):',layer(isoil),layer(isoil+1),' m'
+            if ((val2.gt.layer(isoil)).and.(val2.lt.layer(isoil+1)))
+     &       then
+             iref=isoil+1
+             write(*,*)'Level selected : ',layer(isoil+1),' m'
+             exit
+            endif
+         enddo
+
+         DO ig=1,ngridmx
+             DO j=1,iref
+                   ithfi(ig,j)=ith_bb
+             ENDDO
+         ENDDO
+
+         CALL gr_fi_dyn(nsoilmx,ngridmx,iip1,jjp1,ithfi,ith)
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
         endif