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

Timestamp:

Jun 25, 2024, 1:54:02 PM (18 months ago)

Author:

tbertrand

Message:

LMDZ.PLUTO
resolving some issues in the code for 3D runs
TB

File:

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

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

@@ -149 +149 @@
       real tsurf_bb,tsurf_bb2
       real ptoto,pcap,patm,airetot,ptotn,patmn,psea,geop
-      real tempsoil(22),levsoil(22)
+      real tempsoil(24),levsoil(24)
       character*1 yes
       logical :: flagtset=.false. ,  flagps0=.false.
@@ -872 +872 @@
                ENDDO
              ENDDO
-             write(*,*) 'SURFACE value of tracer ',trim(tname(iq)),
-     &                   ' ? (kg/m2)'
-             read(*,*) val
-             DO ig=1,ngridmx
-                 qsurf(ig,iq)=val
-             ENDDO
+c             write(*,*) 'SURFACE value of tracer ',trim(tname(iq)),
+c     &                   ' ? (kg/m2)'
+c             read(*,*) val
+c             DO ig=1,ngridmx
+c                 qsurf(ig,iq)=val
+c             ENDDO
              
 c       qs=x : initialise surface tracer manually 
@@ -1730 +1730 @@
      &              (phisfi(ig).ge.val6) .and.
      &              (phisfi(ig).le.val11) .and.
-     &              (qsurf(ig,igcm_n2).le.0.001) ) then 
+     &              (qsurf(ig,igcm_n2).gt.0.001) ) then 
 
 !                    DO i=1,ngridmx
@@ -1736 +1736 @@
 !     &              (lonfi(i).eq.0.) ) then
 !
-                         tsurf(ig)=50. 
+                         tsurf(ig)=34.7 
                          !qsurf(ig,igcm_ch4_ice)=qsurf(i,igcm_ch4_ice)
 !
                          DO l=1,nsoilmx
-                           tsoil(ig,l)=50. !tsoil(i,l)
+                           tsoil(ig,l)=34.7 !tsoil(i,l)
                          ENDDO 
                        !ENDIF
@@ -2897 +2897 @@
 
          open(333,file='./tsoil_180_30',form='formatted')
-         do i=1,22
+         do i=1,24
            read(333,*) levsoil(i), tempsoil(i)
          enddo
@@ -3735 +3735 @@
                       tsurf(ig)=field_input(ig)
                       do l=1,nsoilmx
-                         tsoil(ig,l) = field_input(ig)
-                        !tsoil(ig,l) = field_inputs(ig,l)
+                         !tsoil(ig,l) = field_input(ig)
+                         tsoil(ig,l) = field_inputs(ig,l)
                       end do
                    endif