Changeset 4068 for trunk/LMDZ.COMMON/libf/evolution/clim_state_rec.F90

Timestamp:

Feb 12, 2026, 12:32:15 PM (2 weeks ago)

Author:

jbclement

Message:

PEM:
Following r4065, some safeguards forgotten in "io_netcdf.F90" + 'tsurf' is not needed anymore in the "start1D.txt" + few forgotten small updates.
JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/clim_state_rec.F90 (modified) (6 diffs)

trunk/LMDZ.COMMON/libf/evolution/clim_state_rec.F90

@@ -31 +31 @@
 
 !=======================================================================
-SUBROUTINE write_restart(ps4PCM,pa4PCM,preff4PCM,tsurf4PCM,q4PCM,teta4PCM,air_mass4PCM)
+SUBROUTINE write_restart(ps4PCM,pa4PCM,preff4PCM,q4PCM,teta4PCM,air_mass4PCM)
 !-----------------------------------------------------------------------
 ! NAME
@@ -61 +61 @@
 real(dp),                   intent(in) :: pa4PCM, preff4PCM
 real(dp), dimension(:),     intent(in) :: ps4PCM
-real(dp), dimension(:,:),   intent(in) :: tsurf4PCM, teta4PCM, air_mass4PCM
+real(dp), dimension(:,:),   intent(in) :: teta4PCM, air_mass4PCM
 real(dp), dimension(:,:,:), intent(in) :: q4PCM
 
@@ -75 +75 @@
 ! ~~~~~~~~~~~~~
 if (ngrid == 1) then
-    call write_restart1D(ps4PCM,pa4PCM,preff4PCM,tsurf4PCM,q4PCM)
+    call write_restart1D(ps4PCM,pa4PCM,preff4PCM,q4PCM)
     return
 end if
@@ -120 +120 @@
 
 !=======================================================================
-SUBROUTINE write_restart1D(ps4PCM,pa4PCM,preff4PCM,tsurf4PCM,q4PCM)
+SUBROUTINE write_restart1D(ps4PCM,pa4PCM,preff4PCM,q4PCM)
 !-----------------------------------------------------------------------
 ! NAME
@@ -151 +151 @@
 real(dp),                   intent(in) :: pa4PCM, preff4PCM
 real(dp), dimension(:),     intent(in) :: ps4PCM
-real(dp), dimension(:,:),   intent(in) :: tsurf4PCM
 real(dp), dimension(:,:,:), intent(in) :: q4PCM
 
@@ -170 +169 @@
 write(funit,*) 'u', (u_PCM(1,l), l = 1,nlayer)
 write(funit,*) 'v', (v_PCM(1,l), l = 1,nlayer)
-write(funit,*) 'teta', (tsurf4PCM(1,j), j = 1,nslope), (teta_PCM(1,l), l = 1,nlayer)
+write(funit,*) 'teta', (teta_PCM(1,l), l = 1,nlayer)
 close(funit)