Changeset 5451 for LMDZ6/trunk

Timestamp:

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

Author:

fhourdin

Message:

Concerns iophys_ini

Location:

LMDZ6/trunk/libf/phylmd

Files:

• iophys.F90 (modified) (5 diffs)
• lmdz_wake.f90 (modified) (1 diff)
• output_physiqex_mod.f90 (modified) (1 diff)
• yamada_c.F90 (modified) (1 diff)

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))
 

LMDZ6/trunk/libf/phylmd/lmdz_wake.f90

@@ -358 +358 @@
 IF (CPPKEY_IOPHYS_WK) THEN
   IF (phys_sub) THEN
-    call iophys_ini(dtimesub)
+    call iophys_ini(dtimesub,klev)
   ELSE
-    call iophys_ini(dtime)
+    call iophys_ini(dtime,klev)
   ENDIF
 END IF

LMDZ6/trunk/libf/phylmd/output_physiqex_mod.f90

@@ -58 +58 @@
 
    !$OMP MASTER
-   CALL iophys_ini(pdtphys)
+   CALL iophys_ini(pdtphys,klev)
    !$OMP END MASTER
    !$OMP BARRIER

LMDZ6/trunk/libf/phylmd/yamada_c.F90

@@ -138 +138 @@
         CALL getin_p('iflag_tke_diff',iflag_tke_diff)
         allocate(l0(klon))
-#define IOPHYS
-#ifdef IOPHYS
-!        call iophys_ini(timestep)
-#endif
         firstcall=.false.
       endif
 
    IF (ngrid<=0) RETURN ! Bizarre : on n a pas ce probeleme pour coef_diff_turb
-
-#ifdef IOPHYS
-if (okiophys) then
-call iophys_ecrit('q2i',klev,'q2 debut my','m2/s2',q2(:,1:klev))
-call iophys_ecrit('kmi',klev,'Kz debut my','m/s2',km(:,1:klev))
-endif
-#endif
 
       nlay=klev