Changeset 3380 for trunk/LMDZ.PLUTO/libf
- Timestamp:
- Jun 13, 2024, 3:24:15 PM (6 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto/start2archive.F
r3184 r3380 23 23 24 24 ! use slab_ice_h, only: noceanmx 25 USE ocean_slab_mod, ONLY: nslay25 c USE ocean_slab_mod, ONLY: nslay 26 26 USE ioipsl_getincom, only: getin 27 27 USE planete_mod, only: year_day 28 28 USE mod_const_mpi, ONLY: COMM_LMDZ 29 29 USE control_mod, only: planet_type 30 USE callkeys_mod, ONLY: ok_slab_ocean30 c USE callkeys_mod, ONLY: ok_slab_ocean 31 31 use filtreg_mod, only: inifilr 32 32 USE comvert_mod, ONLY: ap,bp … … 36 36 use phys_state_var_mod, only: phys_state_var_init 37 37 use phyetat0_mod, only: phyetat0 38 use nonoro_gwd_ran_mod, only: du_nonoro_gwd, dv_nonoro_gwd,39 & east_gwstress, west_gwstress38 c use nonoro_gwd_ran_mod, only: du_nonoro_gwd, dv_nonoro_gwd, 39 c & east_gwstress, west_gwstress 40 40 use exner_hyb_m, only: exner_hyb 41 41 implicit none … … 136 136 allocate(qsurf(ngridmx,nqtot)) 137 137 allocate(qsurfS(ip1jmp1,nqtot)) 138 allocate(tslab(ngridmx,nslay)) !Added by SB for slab ocean139 allocate(tslabS(ip1jmp1,nslay)) !Added by SB for slab ocean138 c allocate(tslab(ngridmx,nslay)) !Added by SB for slab ocean 139 c allocate(tslabS(ip1jmp1,nslay)) !Added by SB for slab ocean 140 140 ! other array allocations: 141 141 ! call ini_comsoil_h(ngridmx) ! done via iniphysiq … … 222 222 CALL phyetat0(.true.,ngridmx,llm,fichnom,0,Lmodif,nsoilmx,nqtot, 223 223 . day_ini_fi,timefi, 224 . tsurf,tsoil,emis,q2,qsurf ,224 . tsurf,tsoil,emis,q2,qsurf) 225 225 ! change FF 05/2011 226 . cloudfrac,totalcloudfrac,hice,226 c . cloudfrac,totalcloudfrac,hice, 227 227 ! change BC 05/2014 228 . rnat,pctsrf_sic,tslab,tsea_ice,sea_ice)228 c . rnat,pctsrf_sic,tslab,tsea_ice,sea_ice) 229 229 230 230 … … 330 330 call gr_fi_dyn(llm+1,ngridmx,iip1,jjp1,q2,q2S) 331 331 call gr_fi_dyn(nqtot,ngridmx,iip1,jjp1,qsurf,qsurfS) 332 call gr_fi_dyn(llm,ngridmx,iip1,jjp1,cloudfrac,cloudfracS)333 call gr_fi_dyn(1,ngridmx,iip1,jjp1,hice,hiceS)334 call gr_fi_dyn(1,ngridmx,iip1,jjp1,totalcloudfrac,totalcloudfracS)335 336 call gr_fi_dyn(1,ngridmx,iip1,jjp1,rnat,rnatS)337 call gr_fi_dyn(1,ngridmx,iip1,jjp1,pctsrf_sic,pctsrf_sicS)338 call gr_fi_dyn(1,ngridmx,iip1,jjp1,tsea_ice,tsea_iceS)339 call gr_fi_dyn(1,ngridmx,iip1,jjp1,sea_ice,sea_iceS)340 call gr_fi_dyn(nslay,ngridmx,iip1,jjp1,tslab,tslabS)341 342 call gr_fi_dyn(llm,ngridmx,iip1,jjp1,du_nonoro_gwd,du_nonoro_gwdS)343 call gr_fi_dyn(llm,ngridmx,iip1,jjp1,dv_nonoro_gwd,dv_nonoro_gwdS)344 call gr_fi_dyn(llm,ngridmx,iip1,jjp1,east_gwstress,east_gwstressS)345 call gr_fi_dyn(llm,ngridmx,iip1,jjp1,west_gwstress,west_gwstressS)332 c call gr_fi_dyn(llm,ngridmx,iip1,jjp1,cloudfrac,cloudfracS) 333 c call gr_fi_dyn(1,ngridmx,iip1,jjp1,hice,hiceS) 334 c call gr_fi_dyn(1,ngridmx,iip1,jjp1,totalcloudfrac,totalcloudfracS) 335 336 c call gr_fi_dyn(1,ngridmx,iip1,jjp1,rnat,rnatS) 337 c call gr_fi_dyn(1,ngridmx,iip1,jjp1,pctsrf_sic,pctsrf_sicS) 338 c call gr_fi_dyn(1,ngridmx,iip1,jjp1,tsea_ice,tsea_iceS) 339 c call gr_fi_dyn(1,ngridmx,iip1,jjp1,sea_ice,sea_iceS) 340 c call gr_fi_dyn(nslay,ngridmx,iip1,jjp1,tslab,tslabS) 341 342 c call gr_fi_dyn(llm,ngridmx,iip1,jjp1,du_nonoro_gwd,du_nonoro_gwdS) 343 c call gr_fi_dyn(llm,ngridmx,iip1,jjp1,dv_nonoro_gwd,dv_nonoro_gwdS) 344 c call gr_fi_dyn(llm,ngridmx,iip1,jjp1,east_gwstress,east_gwstressS) 345 c call gr_fi_dyn(llm,ngridmx,iip1,jjp1,west_gwstress,west_gwstressS) 346 346 c======================================================================= 347 347 c Info pour controler … … 490 490 c Ecriture du champs cloudfrac,hice,totalcloudfrac 491 491 c----------------------------------------------------------------------- 492 call write_archive(nid,ntime,'hice',493 & 'Height of oceanic ice','m',2,hiceS)494 call write_archive(nid,ntime,'totalcloudfrac',495 & 'Total cloud Fraction','',2,totalcloudfracS)496 call write_archive(nid,ntime,'cloudfrac'497 & ,'Cloud fraction','',3,cloudfracS)492 c call write_archive(nid,ntime,'hice', 493 c & 'Height of oceanic ice','m',2,hiceS) 494 c call write_archive(nid,ntime,'totalcloudfrac', 495 c & 'Total cloud Fraction','',2,totalcloudfracS) 496 c call write_archive(nid,ntime,'cloudfrac' 497 c & ,'Cloud fraction','',3,cloudfracS) 498 498 499 499 c----------------------------------------------------------------------- … … 506 506 IF(ierr.EQ.0) THEN 507 507 508 write(*,*) "Use slab-ocean ?"509 ok_slab_ocean=.false. ! default value510 call getin("ok_slab_ocean",ok_slab_ocean)511 write(*,*) "ok_slab_ocean = ",ok_slab_ocean512 513 if(ok_slab_ocean) then514 call write_archive(nid,ntime,'rnat'515 & ,'rnat','',2,rnatS)516 call write_archive(nid,ntime,'pctsrf_sic'517 & ,'pctsrf_sic','',2,pctsrf_sicS)518 call write_archive(nid,ntime,'sea_ice'519 & ,'sea_ice','',2,sea_iceS)520 call write_archive(nid,ntime,'tslab'521 & ,'tslab','',-2,tslabS)522 call write_archive(nid,ntime,'tsea_ice'523 & ,'tsea_ice','',2,tsea_iceS)524 endif !ok_slab_ocean508 c write(*,*) "Use slab-ocean ?" 509 c ok_slab_ocean=.false. ! default value 510 c call getin("ok_slab_ocean",ok_slab_ocean) 511 c write(*,*) "ok_slab_ocean = ",ok_slab_ocean 512 513 c if(ok_slab_ocean) then 514 c call write_archive(nid,ntime,'rnat' 515 c & ,'rnat','',2,rnatS) 516 c call write_archive(nid,ntime,'pctsrf_sic' 517 c & ,'pctsrf_sic','',2,pctsrf_sicS) 518 c call write_archive(nid,ntime,'sea_ice' 519 c & ,'sea_ice','',2,sea_iceS) 520 c call write_archive(nid,ntime,'tslab' 521 c & ,'tslab','',-2,tslabS) 522 c call write_archive(nid,ntime,'tsea_ice' 523 c & ,'tsea_ice','',2,tsea_iceS) 524 c endif !ok_slab_ocean 525 525 526 526 ENDIF ! of IF(ierr.EQ.0) 527 527 528 528 ! Non-orographic gavity waves 529 call write_archive(nid,ntime,"du_nonoro_gwd",530 & "Zonal wind tendency due to GW",'m.s-1',3,du_nonoro_gwdS)531 call write_archive(nid,ntime,"dv_nonoro_gwd",532 & "Meridional wind tendency due to GW",'m.s-1',533 & 3,dv_nonoro_gwdS)534 call write_archive(nid,ntime,"east_gwstress",535 & "Eastward stress profile due to GW",'kg.m-1.s-2',536 & 3,east_gwstressS)537 call write_archive(nid,ntime,"west_gwstress",538 & "Westward stress profile due to GW",'kg.m-1.s-2',539 & 3,west_gwstressS)529 c call write_archive(nid,ntime,"du_nonoro_gwd", 530 c & "Zonal wind tendency due to GW",'m.s-1',3,du_nonoro_gwdS) 531 c call write_archive(nid,ntime,"dv_nonoro_gwd", 532 c & "Meridional wind tendency due to GW",'m.s-1', 533 c & 3,dv_nonoro_gwdS) 534 c call write_archive(nid,ntime,"east_gwstress", 535 c & "Eastward stress profile due to GW",'kg.m-1.s-2', 536 c & 3,east_gwstressS) 537 c call write_archive(nid,ntime,"west_gwstress", 538 c & "Westward stress profile due to GW",'kg.m-1.s-2', 539 c & 3,west_gwstressS) 540 540 541 541 c-----------------------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.