Changeset 2307 for LMDZ5/trunk/libf/phylmd/dyn1d

Timestamp:

Jun 23, 2015, 1:09:31 PM (10 years ago)

Author:

Laurent Fairhead

Message:

Modifications to run 1D Bomex case
MPL

Location:

LMDZ5/trunk/libf/phylmd/dyn1d

Files:

• 1D_interp_cases.h (modified) (3 diffs)
• 1D_read_forc_cases.h (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/dyn1d/1D_interp_cases.h

@@ +1 @@
+!
+! $Id$
+!
 !---------------------------------------------------------------------
 ! Interpolation forcing in time and onto model levels
@@ -615 +618 @@
      &       ,vitw_cas,du_cas,hu_cas,vu_cas           &
      &       ,dv_cas,hv_cas,vv_cas,dt_cas,ht_cas,vt_cas,dtrad_cas             &
-     &       ,dq_cas,hq_cas,vq_cas,lat_cas,sens_cas                 &
+     &       ,dq_cas,hq_cas,vq_cas,lat_cas,sens_cas,ustar_cas       &
      &       ,ts_prof_cas,plev_prof_cas,t_prof_cas,q_prof_cas,u_prof_cas,v_prof_cas         &
      &       ,ug_prof_cas,vg_prof_cas,vitw_prof_cas,du_prof_cas,hu_prof_cas,vu_prof_cas     &
      &       ,dv_prof_cas,hv_prof_cas,vv_prof_cas,dt_prof_cas,ht_prof_cas,vt_prof_cas       &
      &       ,dtrad_prof_cas,dq_prof_cas,hq_prof_cas,vq_prof_cas,lat_prof_cas               &
-     &       ,sens_prof_cas)
+     &       ,sens_prof_cas,ustar_prof_cas)
 
              ts_cur = ts_prof_cas 
@@ -746 +749 @@
         if (tend_rayo.eq.1) then
            dt_cooling(l) = dtrad_mod_cas(l)
+!          print *,'dt_cooling=',dt_cooling(l)
         else
            dt_cooling(l) = 0.0

LMDZ5/trunk/libf/phylmd/dyn1d/1D_read_forc_cases.h

@@ +1 @@
+!
+! $Id$
+!
 !----------------------------------------------------------------------
 ! forcing_les = .T. : Impose a constant cooling 
@@ -738 +741 @@
      &         ,ug_cas,vg_cas,vitw_cas,du_cas,hu_cas,vu_cas           &
      &         ,dv_cas,hv_cas,vv_cas,dt_cas,ht_cas,vt_cas,dtrad_cas             &
-     &         ,dq_cas,hq_cas,vq_cas,lat_cas,sens_cas                 &
+     &         ,dq_cas,hq_cas,vq_cas,lat_cas,sens_cas,ustar_cas       &
      &         ,ts_prof_cas,plev_prof_cas,t_prof_cas,q_prof_cas,u_prof_cas,v_prof_cas         &
      &         ,ug_prof_cas,vg_prof_cas,vitw_prof_cas,du_prof_cas,hu_prof_cas,vu_prof_cas     &
      &         ,dv_prof_cas,hv_prof_cas,vv_prof_cas,dt_prof_cas,ht_prof_cas,vt_prof_cas,dtrad_prof_cas       &
-     &         ,dq_prof_cas,hq_prof_cas,vq_prof_cas,lat_prof_cas,sens_prof_cas)
+     &         ,dq_prof_cas,hq_prof_cas,vq_prof_cas,lat_prof_cas,sens_prof_cas,ustar_prof_cas)
 
 ! vertical interpolation using TOGA interpolation routine:
@@ -778 +781 @@
        d_u_adv(l) = (hv_mod_cas(l)+vv_mod_cas(l))
       enddo     
-       
+
+! In case fluxes are imposed
+       IF (ok_flux_surf) THEN
+       fsens=sens_prof_cas
+       flat=lat_prof_cas
+       ENDIF
+       IF (ok_prescr_ust) THEN
+       ust=ustar_prof_cas
+       print *,'ust=',ust
+       ENDIF
+
       endif !forcing_case
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!