Changeset 5451 for LMDZ6/trunk/libf/phylmd/iophys.F90

Timestamp:

Dec 23, 2024, 6:08:12 PM (33 hours ago)

Author:

fhourdin

Message:

Concerns iophys_ini

File:

• LMDZ6/trunk/libf/phylmd/iophys.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/iophys.F90

@@ -110 +110 @@
 
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-      SUBROUTINE iophys_ini(timestep)
+      SUBROUTINE iophys_ini(timestep,nlev)
+      USE dimphy, ONLY: klev
       USE mod_phys_lmdz_para, ONLY: is_mpi_root
       USE vertical_layers_mod, ONLY: presnivs
       USE regular_lonlat_mod, ONLY: lon_reg, lat_reg
-      USE dimphy, ONLY: klev
       USE mod_grid_phy_lmdz, ONLY: klon_glo
       USE time_phylmdz_mod, ONLY : annee_ref, day_ref, day_ini
@@ -139 +139 @@
 !   -------------
 
+integer, intent(in) :: nlev
+real, intent(in) :: timestep
+
 real pi
 INTEGER nlat_eff
 INTEGER jour0,mois0,an0
-REAL timestep,t0
+REAL t0
 CHARACTER(len=20) :: calendrier
+integer ilev
+real coord_vert(nlev)
 
 !   Arguments:
@@ -178 +183 @@
 print*,'iophys_ini annee_ref day_ref',annee_ref,day_ref,day_ini,calend,t0
 
-
+if ( nlev == klev ) then
+     coord_vert=presnivs
+print*,'ON EST LA '
+else
+     do ilev=1,nlev
+        coord_vert(ilev)=ilev
+     enddo
+endif
+print*,'nlev=',nlev
+print*,'coord_vert',coord_vert
 call iotd_ini('phys.nc', &
-size(lon_reg),nlat_eff,klev,lon_reg(:)*180./pi,lat_reg*180./pi,presnivs,jour0,mois0,an0,t0,timestep,calendrier)
+size(lon_reg),nlat_eff,nlev,lon_reg(:)*180./pi,lat_reg*180./pi,coord_vert,jour0,mois0,an0,t0,timestep,calendrier)
+ !  SUBROUTINE iotd_ini(fichnom,iim,jjm,llm,prlon,prlat,pcoordv,jour0,mois0,an0,t0,timestep,calendrier)
+!   -------
     ENDIF
 !$OMP END MASTER
@@ -216 +232 @@
 
       SUBROUTINE iotd_ecrit_seq(nom,lllm,titre,unite,px)
+!call iotd_ecrit_seq('f0',1,'f0 in thermcell_plume_6A',' ',f0(1:ngrid))
+
         USE iotd_mod_h
 
@@ -230 +248 @@
       integer i,j,l,ijl
 
+      !print*,'iotd_ecrit_seq ,nom,lllm,titre,unite,px',nom,lllm,titre,unite,px
       allocate(zx(imax,jmax,lllm))