- Timestamp:
- Jan 11, 2021, 11:24:08 PM (4 years ago)
- Location:
- LMDZ6/branches/Ocean_skin
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Ocean_skin
-
LMDZ6/branches/Ocean_skin/libf/phylmd/dyn1d/1D_read_forc_cases.h
r3605 r3798 27 27 & ,u_cas,v_cas,ug_cas,vg_cas & 28 28 & ,temp_nudg_cas,qv_nudg_cas,u_nudg_cas,v_nudg_cas & 29 & ,vitw_cas,omega_cas, du_cas,hu_cas,vu_cas &29 & ,vitw_cas,omega_cas,tke_cas,du_cas,hu_cas,vu_cas & 30 30 & ,dv_cas,hv_cas,vv_cas,dt_cas,ht_cas,vt_cas,dtrad_cas & 31 31 & ,dq_cas,hq_cas,vq_cas,dth_cas,hth_cas,vth_cas,lat_cas,sens_cas,ustar_cas & 32 & ,uw_cas,vw_cas,q1_cas,q2_cas,tke _cas &32 & ,uw_cas,vw_cas,q1_cas,q2_cas,tkes_cas & 33 33 ! 34 34 & ,ts_prof_cas,ps_prof_cas,plev_prof_cas,t_prof_cas,theta_prof_cas,thv_prof_cas & … … 36 36 & ,u_prof_cas,v_prof_cas,ug_prof_cas,vg_prof_cas & 37 37 & ,temp_nudg_prof_cas,qv_nudg_prof_cas,u_nudg_prof_cas,v_nudg_prof_cas & 38 & ,vitw_prof_cas,omega_prof_cas 38 & ,vitw_prof_cas,omega_prof_cas,tke_prof_cas & 39 39 & ,du_prof_cas,hu_prof_cas,vu_prof_cas & 40 40 & ,dv_prof_cas,hv_prof_cas,vv_prof_cas,dt_prof_cas,ht_prof_cas,vt_prof_cas & 41 41 & ,dtrad_prof_cas,dq_prof_cas,hq_prof_cas,vq_prof_cas & 42 42 & ,dth_prof_cas,hth_prof_cas,vth_prof_cas,lat_prof_cas & 43 & ,sens_prof_cas,ustar_prof_cas,uw_prof_cas,vw_prof_cas,q1_prof_cas,q2_prof_cas,tke _prof_cas)43 & ,sens_prof_cas,ustar_prof_cas,uw_prof_cas,vw_prof_cas,q1_prof_cas,q2_prof_cas,tkes_prof_cas) 44 44 45 45 do l = 1, nlev_cas … … 49 49 ! vertical interpolation using interpolation routine: 50 50 ! write(*,*)'avant interp vert', t_prof 51 CALL interp2_case_vertical_std(play, nlev_cas,plev_prof_cas &51 CALL interp2_case_vertical_std(play,plev,nlev_cas,plev_prof_cas & 52 52 & ,t_prof_cas,theta_prof_cas,thv_prof_cas,thl_prof_cas & 53 53 & ,qv_prof_cas,ql_prof_cas,qi_prof_cas,u_prof_cas,v_prof_cas & … … 55 55 & ,temp_nudg_prof_cas,qv_nudg_prof_cas,u_nudg_prof_cas,v_nudg_prof_cas & 56 56 57 & ,vitw_prof_cas,omega_prof_cas 57 & ,vitw_prof_cas,omega_prof_cas,tke_prof_cas & 58 58 & ,du_prof_cas,hu_prof_cas,vu_prof_cas,dv_prof_cas,hv_prof_cas,vv_prof_cas & 59 59 & ,dt_prof_cas,ht_prof_cas,vt_prof_cas,dtrad_prof_cas,dq_prof_cas,hq_prof_cas,vq_prof_cas & … … 63 63 & ,u_mod_cas,v_mod_cas,ug_mod_cas,vg_mod_cas & 64 64 & ,temp_nudg_mod_cas,qv_nudg_mod_cas,u_nudg_mod_cas,v_nudg_mod_cas & 65 & ,w_mod_cas,omega_mod_cas 65 & ,w_mod_cas,omega_mod_cas,tke_mod_cas & 66 66 & ,du_mod_cas,hu_mod_cas,vu_mod_cas,dv_mod_cas,hv_mod_cas,vv_mod_cas & 67 67 & ,dt_mod_cas,ht_mod_cas,vt_mod_cas,dtrad_mod_cas,dq_mod_cas,hq_mod_cas,vq_mod_cas & … … 70 70 71 71 ! initial and boundary conditions : 72 ! 72 ! tsurf = ts_prof_cas 73 73 psurf = ps_prof_cas 74 ts_cur = ts_prof_cas 74 !EV tg instead of ts_cur 75 tg = ts_prof_cas 76 print*, 'tg=', tg 77 75 78 do l = 1, llm 76 79 temp(l) = t_mod_cas(l) … … 95 98 d_u_adv(l) = du_mod_cas(l)+hu_mod_cas(l)+vu_mod_cas(l) 96 99 d_v_adv(l) = dv_mod_cas(l)+hv_mod_cas(l)+vv_mod_cas(l) 100 enddo 97 101 98 ! print*,'d_t_adv ',d_t_adv(1:20)*86400102 ! Etienne pour initialisation de TKE 99 103 100 enddo 104 do l=1,llm+1 105 pbl_tke(:,l,:)=tke_mod_cas(l) 106 enddo 101 107 102 108 ! Faut-il multiplier par -1 ? (MPL 20160713) … … 108 114 IF (ok_prescr_ust) THEN 109 115 ust=ustar_prof_cas 110 print *,'ust=',ust111 116 ENDIF 112 117 118 113 119 endif !forcing_SCM
Note: See TracChangeset
for help on using the changeset viewer.