Changeset 2307 for LMDZ5/trunk/libf/phylmd/mod_1D_cases_read.F90

Timestamp:

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

Author:

Laurent Fairhead

Message:

Modifications to run 1D Bomex case
MPL

File:

• LMDZ5/trunk/libf/phylmd/mod_1D_cases_read.F90 (modified) (18 diffs)

LMDZ5/trunk/libf/phylmd/mod_1D_cases_read.F90

@@ +1 @@
+!
+! $Id$
+!
 MODULE mod_1D_cases_read
 
@@ -11 +14 @@
         real heure_ini_cas
         real day_ju_ini_cas   ! Julian day of case first day
-        parameter (year_ini_cas=2011)
-        parameter (mth_ini_cas=10)
-        parameter (day_ini_cas=1)  ! 10 = 10Juil2006
+!       parameter (year_ini_cas=2011)
+        parameter (year_ini_cas=1969)
+!       parameter (mth_ini_cas=10)
+        parameter (mth_ini_cas=6)
+!       parameter (day_ini_cas=1)  ! 10 = 10Juil2006
+        parameter (day_ini_cas=24)  ! 24 = 24 juin 1969
         parameter (heure_ini_cas=0.) !0h en secondes
         real pdt_cas
@@ -43 +49 @@
         real, allocatable::  vitw_cas(:,:)
         real, allocatable::  ug_cas(:,:),vg_cas(:,:)
-        real, allocatable::  lat_cas(:),sens_cas(:),ts_cas(:)
+        real, allocatable::  lat_cas(:),sens_cas(:),ts_cas(:),ustar_cas(:)
 
 !champs interpoles
@@ -70 +76 @@
 
 
-        real lat_prof_cas,sens_prof_cas,ts_prof_cas
+        real lat_prof_cas,sens_prof_cas,ts_prof_cas,ustar_prof_cas
       
 
@@ -143 +149 @@
         allocate(ug_cas(nlev_cas,nt_cas))
         allocate(vg_cas(nlev_cas,nt_cas))
-        allocate(lat_cas(nt_cas),sens_cas(nt_cas),ts_cas(nt_cas))
+        allocate(lat_cas(nt_cas),sens_cas(nt_cas),ts_cas(nt_cas),ustar_cas(nt_cas))
 
 
@@ -175 +181 @@
      &     ,ug_cas,vg_cas,vitw_cas,du_cas,hu_cas,vu_cas,dv_cas,hv_cas,vv_cas &
      &     ,dt_cas,dtrad_cas,ht_cas,vt_cas,dq_cas,hq_cas,vq_cas             &
-     &     ,dth_cas,hth_cas,vth_cas,dr_cas,hr_cas,vr_cas,sens_cas,lat_cas,ts_cas)
+     &     ,dth_cas,hth_cas,vth_cas,dr_cas,hr_cas,vr_cas,sens_cas,lat_cas,ts_cas,ustar_cas)
 
 
@@ -203 +209 @@
         deallocate(ug_cas)
         deallocate(vg_cas)
-        deallocate(lat_cas,sens_cas,ts_cas)
+        deallocate(lat_cas,sens_cas,ts_cas,ustar_cas)
 
 !champs interpoles
@@ -241 +247 @@
      &     ,zz,pp,temp,qv,rh,theta,rv,u,v,ug,vg,w,                   &
      &     du,hu,vu,dv,hv,vv,dt,dtrad,ht,vt,dq,hq,vq,                     &
-     &     dth,hth,vth,dr,hr,vr,sens,flat,ts)
+     &     dth,hth,vth,dr,hr,vr,sens,flat,ts,ustar)
 
 !program reading forcing of the case study
@@ -265 +271 @@
       real dth(nlevel,ntime),hth(nlevel,ntime),vth(nlevel,ntime)
       real dr(nlevel,ntime),hr(nlevel,ntime),vr(nlevel,ntime)
-      real flat(ntime),sens(ntime),ts(ntime)
+      real flat(ntime),sens(ntime),ts(ntime),ustar(ntime)
 
 
       integer nid, ierr,rid
       integer nbvar3d
-      parameter(nbvar3d=34)
+      parameter(nbvar3d=39)
       integer var3didin(nbvar3d)
 
@@ -478 +484 @@
            stop 'ts'
          endif
+
+      ierr=NF_INQ_VARID(nid,"ustar",var3didin(35))
+         if(ierr/=NF_NOERR) then
+           write(*,*) NF_STRERROR(ierr)
+           stop 'ustar'
+         endif
  
 #ifdef NC_DOUBLE
@@ -853 +865 @@
          endif
 !          write(*,*)'lecture ts ok',ts
+
+#ifdef NC_DOUBLE
+         ierr = NF_GET_VAR_DOUBLE(nid,var3didin(35),ustar)
+#else
+         ierr = NF_GET_VAR_REAL(nid,var3didin(35),ustar)
+#endif
+         if(ierr/=NF_NOERR) then
+            write(*,*) NF_STRERROR(ierr)
+            stop "getvarup"
+         endif
+!         write(*,*)'lecture ustar ok',ustar
+
 
          return 
@@ -862 +886 @@
      &         ,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         &
@@ -868 +892 @@
      &         ,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)
+     &         ,hq_prof_cas,vq_prof_cas,lat_prof_cas,sens_prof_cas,ustar_prof_cas)
           
 
@@ -903 +927 @@
         real lat_cas(nt_cas)
         real sens_cas(nt_cas)
+        real ustar_cas(nt_cas)
 
 ! outputs:
@@ -915 +940 @@
         real dtrad_prof_cas(nlev_cas)
         real dq_prof_cas(nlev_cas),hq_prof_cas(nlev_cas),vq_prof_cas(nlev_cas)
-        real lat_prof_cas,sens_prof_cas,ts_prof_cas
+        real lat_prof_cas,sens_prof_cas,ts_prof_cas,ustar_prof_cas
 ! local:
         integer it_cas1, it_cas2,k
@@ -921 +946 @@
 
 
-        print*,'Check time',day1,day_ini_cas,day_ini_cas+90
+        print*,'Check time',day1,day_ini_cas,day_ini_cas+1
 
         if ((forcing_type.eq.10).and.(1.eq.1)) then
 ! Check that initial day of the simulation consistent with the case:
-       if (annee_ref.ne.2011) then
-        print*,'Pour CINDY, annee_ref doit etre 2011'
+       if (annee_ref.ne.1969) then
+!       print*,'Pour CINDY, annee_ref doit etre 2011'
+        print*,'Pour BOMEX, annee_ref doit etre 1969'
         print*,'Changer annee_ref dans run.def'
         stop
        endif
-       if (annee_ref.eq.2011 .and. day1.lt.day_ini_cas) then
-        print*,'CINDY a débuté le 1 octobre 2011',day1,day_ini_cas
+       if (annee_ref.eq.1969 .and. day1.lt.day_ini_cas) then
+!       print*,'CINDY a debute le 1 octobre 2011',day1,day_ini_cas
+        print*,'BOMEX a debute le 24 juin 1969',day1,day_ini_cas
         print*,'Changer dayref dans run.def'
         stop
        endif
-       if (annee_ref.eq.2011 .and. day1.gt.day_ini_cas+90) then
-        print*,'CINDY a fini le 31 decembre'
+       if (annee_ref.eq.1969 .and. day1.gt.day_ini_cas+1) then
+!       print*,'CINDY a fini le 31 decembre'
+        print*,'BOMEX a fini le 24 juin'
         print*,'Changer dayref ou nday dans run.def',day1,day_ini_cas+90
         stop
@@ -950 +978 @@
        endif
        if (annee_ref.eq.2006 .and. day1.lt.day_ini_cas) then
-        print*,'AMMA a débuté le 10 juillet 2006',day1,day_ini_cas
+        print*,'AMMA a debute le 10 juillet 2006',day1,day_ini_cas
         print*,'Changer dayref dans run.def'
         stop
@@ -999 +1027 @@
        sens_prof_cas = sens_cas(it_cas2)                                     &
      &          -frac*(sens_cas(it_cas2)-sens_cas(it_cas1))
-       ts_prof_cas = ts_cas(it_cas2)                                     &
+       ts_prof_cas = ts_cas(it_cas2)                                         &
      &          -frac*(ts_cas(it_cas2)-ts_cas(it_cas1))
+       ustar_prof_cas = ustar_cas(it_cas2)                                   &
+     &          -frac*(ustar_cas(it_cas2)-ustar_cas(it_cas1))
 
        do k=1,nlev_cas