Changeset 2758

Timestamp:

Jan 10, 2017, 3:36:36 PM (8 years ago)

Author:

jyg

Message:

Bug fix in dyn1d

File:

• LMDZ5/trunk/libf/phylmd/dyn1d/mod_1D_cases_read2.F90 (modified) (9 diffs)

LMDZ5/trunk/libf/phylmd/dyn1d/mod_1D_cases_read2.F90

@@ -178 +178 @@
 
         print*,'Allocations OK'
-        call read_cas(nid,nlev_cas,nt_cas                                       &
+        call read_cas2(nid,nlev_cas,nt_cas                                       &
      &     ,z_cas,plev_cas,t_cas,q_cas,rh_cas,th_cas,rv_cas,u_cas,v_cas         &
      &     ,ug_cas,vg_cas,vitw_cas,du_cas,hu_cas,vu_cas,dv_cas,hv_cas,vv_cas    &
@@ -382 +382 @@
 END MODULE mod_1D_cases_read2
 !=====================================================================
-      subroutine read_cas(nid,nlevel,ntime                          &
+      subroutine read_cas2(nid,nlevel,ntime                          &
      &     ,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,                     &
@@ -423 +423 @@
 
        do i=1,nbvar3d
-         print *,'Dans read_cas, on va lire ',nid,i,name_var(i)
+         print *,'Dans read_cas2, on va lire ',nid,i,name_var(i)
        enddo
        do i=1,nbvar3d
@@ -433 +433 @@
        enddo
        do i=1,nbvar3d
-         print *,'Dans read_cas, on va lire ',var3didin(i),name_var(i)
+         print *,'Dans read_cas2, on va lire ',var3didin(i),name_var(i)
          if(i.LE.35) then
 #ifdef NC_DOUBLE
@@ -440 +440 @@
          ierr = NF_GET_VAR_REAL(nid,var3didin(i),resul)
 #endif
-         print *,'Dans read_cas, on a lu ',ierr,var3didin(i),name_var(i)
+         print *,'Dans read_cas2, on a lu ',ierr,var3didin(i),name_var(i)
          if(ierr/=NF_NOERR) then
             print *,'Pb a la lecture de cas.nc: ',name_var(i)
@@ -451 +451 @@
          ierr = NF_GET_VAR_REAL(nid,var3didin(i),resul1)
 #endif
-         print *,'Dans read_cas, on a lu ',ierr,var3didin(i),name_var(i)
+         print *,'Dans read_cas2, on a lu ',ierr,var3didin(i),name_var(i)
          if(ierr/=NF_NOERR) then
             print *,'Pb a la lecture de cas.nc: ',name_var(i)
@@ -501 +501 @@
 
          return 
-         end subroutine read_cas
+         end subroutine read_cas2
 !======================================================================
       subroutine read2_cas(nid,nlevel,ntime,                                       &
@@ -686 +686 @@
          end subroutine read2_cas
 !======================================================================
-        SUBROUTINE interp_case_time(day,day1,annee_ref                &
+        SUBROUTINE interp_case_time2(day,day1,annee_ref                &
 !    &         ,year_cas,day_cas,nt_cas,pdt_forc,nlev_cas      &
      &         ,nt_cas,nlev_cas                                       &
@@ -887 +887 @@
 
         return
-        END SUBROUTINE interp_case_time
+        END SUBROUTINE interp_case_time2
 
 !**********************************************************************************************