Changeset 4584 for LMDZ6/branches/Portage_acc/libf/phylmd/dyn1d
- Timestamp:
- Jun 23, 2023, 4:10:41 PM (15 months ago)
- Location:
- LMDZ6/branches/Portage_acc/libf/phylmd/dyn1d
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Portage_acc/libf/phylmd/dyn1d/mod_1D_cases_read2.F90
r4446 r4584 1135 1135 ! endif 1136 1136 timeit=(day-day_ju_ini_cas)*86400 1137 print *,'day=',day1138 print *,'day_ju_ini_cas=',day_ju_ini_cas1139 print *,'pdt_cas=',pdt_cas1140 print *,'timeit=',timeit1141 print *,'nt_cas=',nt_cas1137 !print *,'day=',day 1138 !print *,'day_ju_ini_cas=',day_ju_ini_cas 1139 !print *,'pdt_cas=',pdt_cas 1140 !print *,'timeit=',timeit 1141 !print *,'nt_cas=',nt_cas 1142 1142 1143 1143 ! Determine the closest observation times: … … 1155 1155 time_cas1=(it_cas1-1)*pdt_cas 1156 1156 time_cas2=(it_cas2-1)*pdt_cas 1157 print *,'it_cas1,it_cas2,time_cas1,time_cas2=',it_cas1,it_cas2,time_cas1,time_cas21157 !print *,'it_cas1,it_cas2,time_cas1,time_cas2=',it_cas1,it_cas2,time_cas1,time_cas2 1158 1158 1159 1159 if (it_cas1 .gt. nt_cas) then … … 1350 1350 ! endif 1351 1351 timeit=(day-day_ju_ini_cas)*86400 1352 print *,'day=',day1353 print *,'day_ju_ini_cas=',day_ju_ini_cas1354 print *,'pdt_cas=',pdt_cas1355 print *,'timeit=',timeit1356 print *,'nt_cas=',nt_cas1352 !print *,'day=',day 1353 !print *,'day_ju_ini_cas=',day_ju_ini_cas 1354 !print *,'pdt_cas=',pdt_cas 1355 !print *,'timeit=',timeit 1356 !print *,'nt_cas=',nt_cas 1357 1357 1358 1358 ! Determine the closest observation times: … … 1370 1370 time_cas1=(it_cas1-1)*pdt_cas 1371 1371 time_cas2=(it_cas2-1)*pdt_cas 1372 print *,'timeit,pdt_cas,nt_cas=',timeit,pdt_cas,nt_cas1373 print *,'it_cas1,it_cas2,time_cas1,time_cas2=',it_cas1,it_cas2,time_cas1,time_cas21372 !print *,'timeit,pdt_cas,nt_cas=',timeit,pdt_cas,nt_cas 1373 !print *,'it_cas1,it_cas2,time_cas1,time_cas2=',it_cas1,it_cas2,time_cas1,time_cas2 1374 1374 1375 1375 if (it_cas1 .gt. nt_cas) then … … 1403 1403 t_prof_cas(k) = t_cas(k,it_cas2) & 1404 1404 -frac*(t_cas(k,it_cas2)-t_cas(k,it_cas1)) 1405 print *,'k,frac,plev_cas1,plev_cas2=',k,frac,plev_cas(k,it_cas1),plev_cas(k,it_cas2)1405 !print *,'k,frac,plev_cas1,plev_cas2=',k,frac,plev_cas(k,it_cas1),plev_cas(k,it_cas2) 1406 1406 theta_prof_cas(k) = theta_cas(k,it_cas2) & 1407 1407 -frac*(theta_cas(k,it_cas2)-theta_cas(k,it_cas1)) -
LMDZ6/branches/Portage_acc/libf/phylmd/dyn1d/old_lmdz1d.F90
r4446 r4584 13 13 ftsol, beta_aridity, pbl_tke, pctsrf, radsol, rain_fall, snow_fall, ratqs, & 14 14 rnebcon, rugoro, sig1, w01, solaire_etat0, sollw, sollwdown, & 15 solsw, solswfdiff, t_ancien, q_ancien, u_ancien, v_ancien, &15 solsw, solswfdiff, t_ancien, q_ancien, u_ancien, v_ancien, rneb_ancien, & 16 16 wake_delta_pbl_TKE, delta_tsurf, wake_fip, wake_pe, & 17 17 wake_deltaq, wake_deltat, wake_s, wake_dens, & … … 920 920 u_ancien(1,:)=u(:) 921 921 v_ancien(1,:)=v(:) 922 rneb_ancien(1,:)=0. 922 923 923 924 u10m=0. … … 945 946 !------------------------------------------------------------------------ 946 947 !Al1 =============== restart option ========================== 948 iflag_physiq=0 949 call getin('iflag_physiq',iflag_physiq) 950 947 951 if (.not.restart) then 948 952 iflag_pbl = 5 … … 1248 1252 !! Increment state variables 1249 1253 !!!!!!!!!!!!!!!!!!!!!!!!!!!!! 1250 print*,'OLDLMDZ1D AANT'1251 1254 IF (flag_inhib_forcing == 0) then ! if tendency of forcings should be added 1252 print*,'OLDLMDZ1D ARES' , forcing_sandu1253 1255 1254 1256 ! pour les cas sandu et astex, on reclacule u,v,q,temp et teta dans 1D_nudge_sandu_astex.h -
LMDZ6/branches/Portage_acc/libf/phylmd/dyn1d/scm.F90
r4446 r4584 9 9 ftsol, beta_aridity, pbl_tke, pctsrf, radsol, rain_fall, snow_fall, ratqs, & 10 10 rnebcon, rugoro, sig1, w01, solaire_etat0, sollw, sollwdown, & 11 solsw, solswfdiff, t_ancien, q_ancien, u_ancien, v_ancien, &11 solsw, solswfdiff, t_ancien, q_ancien, u_ancien, v_ancien, rneb_ancien, & 12 12 wake_delta_pbl_TKE, delta_tsurf, wake_fip, wake_pe, & 13 13 wake_deltaq, wake_deltat, wake_s, wake_dens, & … … 678 678 u_ancien(1,:)=u(:) 679 679 v_ancien(1,:)=v(:) 680 rneb_ancien(1,:)=0. 680 681 681 682 u10m=0. … … 703 704 !------------------------------------------------------------------------ 704 705 !Al1 =============== restart option ====================================== 706 iflag_physiq=0 707 call getin('iflag_physiq',iflag_physiq) 708 705 709 if (.not.restart) then 706 710 iflag_pbl = 5
Note: See TracChangeset
for help on using the changeset viewer.