Changeset 5451 for LMDZ6/trunk/libf/phylmd
- Timestamp:
- Dec 23, 2024, 6:08:12 PM (45 hours ago)
- Location:
- LMDZ6/trunk/libf/phylmd
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/iophys.F90
r5390 r5451 110 110 111 111 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 112 SUBROUTINE iophys_ini(timestep) 112 SUBROUTINE iophys_ini(timestep,nlev) 113 USE dimphy, ONLY: klev 113 114 USE mod_phys_lmdz_para, ONLY: is_mpi_root 114 115 USE vertical_layers_mod, ONLY: presnivs 115 116 USE regular_lonlat_mod, ONLY: lon_reg, lat_reg 116 USE dimphy, ONLY: klev117 117 USE mod_grid_phy_lmdz, ONLY: klon_glo 118 118 USE time_phylmdz_mod, ONLY : annee_ref, day_ref, day_ini … … 139 139 ! ------------- 140 140 141 integer, intent(in) :: nlev 142 real, intent(in) :: timestep 143 141 144 real pi 142 145 INTEGER nlat_eff 143 146 INTEGER jour0,mois0,an0 144 REAL t imestep,t0147 REAL t0 145 148 CHARACTER(len=20) :: calendrier 149 integer ilev 150 real coord_vert(nlev) 146 151 147 152 ! Arguments: … … 178 183 print*,'iophys_ini annee_ref day_ref',annee_ref,day_ref,day_ini,calend,t0 179 184 180 185 if ( nlev == klev ) then 186 coord_vert=presnivs 187 print*,'ON EST LA ' 188 else 189 do ilev=1,nlev 190 coord_vert(ilev)=ilev 191 enddo 192 endif 193 print*,'nlev=',nlev 194 print*,'coord_vert',coord_vert 181 195 call iotd_ini('phys.nc', & 182 size(lon_reg),nlat_eff,klev,lon_reg(:)*180./pi,lat_reg*180./pi,presnivs,jour0,mois0,an0,t0,timestep,calendrier) 196 size(lon_reg),nlat_eff,nlev,lon_reg(:)*180./pi,lat_reg*180./pi,coord_vert,jour0,mois0,an0,t0,timestep,calendrier) 197 ! SUBROUTINE iotd_ini(fichnom,iim,jjm,llm,prlon,prlat,pcoordv,jour0,mois0,an0,t0,timestep,calendrier) 198 ! ------- 183 199 ENDIF 184 200 !$OMP END MASTER … … 216 232 217 233 SUBROUTINE iotd_ecrit_seq(nom,lllm,titre,unite,px) 234 !call iotd_ecrit_seq('f0',1,'f0 in thermcell_plume_6A',' ',f0(1:ngrid)) 235 218 236 USE iotd_mod_h 219 237 … … 230 248 integer i,j,l,ijl 231 249 250 !print*,'iotd_ecrit_seq ,nom,lllm,titre,unite,px',nom,lllm,titre,unite,px 232 251 allocate(zx(imax,jmax,lllm)) 233 252 -
LMDZ6/trunk/libf/phylmd/lmdz_wake.f90
r5400 r5451 358 358 IF (CPPKEY_IOPHYS_WK) THEN 359 359 IF (phys_sub) THEN 360 call iophys_ini(dtimesub )360 call iophys_ini(dtimesub,klev) 361 361 ELSE 362 call iophys_ini(dtime )362 call iophys_ini(dtime,klev) 363 363 ENDIF 364 364 END IF -
LMDZ6/trunk/libf/phylmd/output_physiqex_mod.f90
r5268 r5451 58 58 59 59 !$OMP MASTER 60 CALL iophys_ini(pdtphys )60 CALL iophys_ini(pdtphys,klev) 61 61 !$OMP END MASTER 62 62 !$OMP BARRIER -
LMDZ6/trunk/libf/phylmd/yamada_c.F90
r5390 r5451 138 138 CALL getin_p('iflag_tke_diff',iflag_tke_diff) 139 139 allocate(l0(klon)) 140 #define IOPHYS141 #ifdef IOPHYS142 ! call iophys_ini(timestep)143 #endif144 140 firstcall=.false. 145 141 endif 146 142 147 143 IF (ngrid<=0) RETURN ! Bizarre : on n a pas ce probeleme pour coef_diff_turb 148 149 #ifdef IOPHYS150 if (okiophys) then151 call iophys_ecrit('q2i',klev,'q2 debut my','m2/s2',q2(:,1:klev))152 call iophys_ecrit('kmi',klev,'Kz debut my','m/s2',km(:,1:klev))153 endif154 #endif155 144 156 145 nlay=klev
Note: See TracChangeset
for help on using the changeset viewer.