Changeset 2683 for LMDZ5/trunk/libf/phylmd/dyn1d/1D_read_forc_cases.h
- Timestamp:
- Oct 24, 2016, 6:46:30 PM (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/dyn1d/1D_read_forc_cases.h
r2672 r2683 367 367 fich_dice='dice_driver.nc' 368 368 call read_dice(fich_dice,nlev_dice,nt_dice & 369 & ,zz_dice,plev_dice,t h_dice,qv_dice,u_dice,v_dice,o3_dice &369 & ,zz_dice,plev_dice,t_dice,qv_dice,u_dice,v_dice,o3_dice & 370 370 & ,shf_dice,lhf_dice,lwup_dice,swup_dice,tg_dice,ustar_dice& 371 371 & ,psurf_dice,ug_dice,vg_dice,ht_dice,hq_dice & … … 376 376 !champs initiaux: 377 377 do k=1,nlev_dice 378 t h_dicei(k)=th_dice(k)378 t_dicei(k)=t_dice(k) 379 379 qv_dicei(k)=qv_dice(k) 380 380 u_dicei(k)=u_dice(k) … … 405 405 406 406 CALL interp_dice_vertical(play,nlev_dice,nt_dice,plev_dice & 407 & ,t h_dicei,qv_dicei,u_dicei,v_dicei,o3_dicei &407 & ,t_dicei,qv_dicei,u_dicei,v_dicei,o3_dicei & 408 408 & ,ht_dicei,hq_dicei,hu_dicei,hv_dicei,w_dicei,omega_dicei& 409 & ,t h_mod,qv_mod,u_mod,v_mod,o3_mod &409 & ,t_mod,qv_mod,u_mod,v_mod,o3_mod & 410 410 & ,ht_mod,hq_mod,hu_mod,hv_mod,w_mod,omega_mod,mxcalc) 411 411 … … 425 425 do l = 1, llm 426 426 ! Ligne du dessous ?? decommenter si on lit theta au lieu de temp 427 428 !temp(l) = t_mod(l)427 ! temp(l) = th_mod(l)*(play(l)/pzero)**rkappa 428 temp(l) = t_mod(l) 429 429 q(l,1) = qv_mod(l) 430 430 q(l,2) = 0.0 … … 476 476 !--------------------------------------------------------------------- 477 477 478 !!!! Si la temperature de surface n'est pas impos ée:478 !!!! Si la temperature de surface n'est pas impos??e: 479 479 480 480 if (forcing_gabls4) then … … 531 531 write(*,*) 'SST initiale mxcalc: ',tsurf,mxcalc 532 532 do l = 1, llm 533 ! Ligne du dessous àdecommenter si on lit theta au lieu de temp533 ! Ligne du dessous ?? decommenter si on lit theta au lieu de temp 534 534 ! temp(l) = th_mod(l)*(play(l)/pzero)**rkappa 535 535 temp(l) = t_mod(l) … … 909 909 endif !forcing_case 910 910 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 911
Note: See TracChangeset
for help on using the changeset viewer.