Changeset 1208 for trunk/LMDZ.MARS/libf/phymars/lect_start_archive.F

Timestamp:

Mar 12, 2014, 4:30:18 PM (11 years ago)

Author:

tnavarro

Message:

added tauscaling in startfi + moved start_archive routines from dyn3d to phymars

File:

• trunk/LMDZ.MARS/libf/phymars/lect_start_archive.F (moved) (moved from trunk/LMDZ.MARS/libf/dyn3d/lect_start_archive.F) (9 diffs)

trunk/LMDZ.MARS/libf/phymars/lect_start_archive.F

@@ -2 +2 @@
      &     date,tsurf,tsoil,emis,q2,
      &     t,ucov,vcov,ps,co2ice,h,phisold_newgrid,
-     &     q,qsurf,surfith,nid)
+     &     q,qsurf,tauscaling,surfith,nid)
 c=======================================================================
 c
@@ -109 +109 @@
       REAL emis(ngrid)
       REAL q2(ngrid,nlayer+1),qsurf(ngrid,nqtot)
+      REAL tauscaling(ngrid) ! dust conversion factor
 c     REAL phisfi(ngrid)
 
@@ -136 +137 @@
       real co2iceS(iip1,jjp1),emisS(iip1,jjp1)
       REAL q2S(iip1,jjp1,llm+1),qsurfS(iip1,jjp1,nqtot)
+      real tauscalingS(iip1,jjp1)
 
       real ptotal, co2icetotal
@@ -163 +165 @@
       real, dimension(:,:), allocatable :: emisold
       real, dimension(:,:,:,:), allocatable :: qold
+      real, dimension(:,:), allocatable :: tauscalingold
 
       real tab_cntrl(100)
@@ -337 +340 @@
       allocate(mlayerold(nsoilold))
       allocate(qsurfold(imold+1,jmold+1,nqtot))
+      allocate(tauscalingold(imold+1,jmold+1))
 
       allocate(var (imold+1,jmold+1,llm))
@@ -662 +666 @@
 
 c-----------------------------------------------------------------------
-c 5.1 Lecture des champs 2D (co2ice, emis,ps,tsurf,Tg[10], q2surf)
+c 5.1 Lecture des champs 2D (co2ice, emis,ps,tsurf,Tg[10], q2surf, tauscaling)
 c-----------------------------------------------------------------------
  
@@ -742 +746 @@
          PRINT*, "lect_start_archive: Failed loading <q2surf>"
          CALL abort
+      ENDIF
+c
+      ierr = NF_INQ_VARID (nid, "tauscaling", nvarid)
+      IF (ierr .NE. NF_NOERR) THEN
+         PRINT*, "lect_start_archive: <tauscaling> not in file"
+         tauscalingold(:,:) = -1
+      ELSE
+#ifdef NC_DOUBLE
+        ierr = NF_GET_VARA_DOUBLE(nid, nvarid,start,count,tauscalingold)
+#else
+        ierr = NF_GET_VARA_REAL(nid, nvarid,start,count,tauscalingold)
+#endif
+        IF (ierr .NE. NF_NOERR) THEN
+           PRINT*, "lect_start_archive: Failed loading <tauscaling>"
+           PRINT*, NF_STRERROR(ierr)
+           CALL abort
+        ENDIF
       ENDIF
 c
@@ -1019 +1040 @@
       call interp_horiz (emisold,emiss,imold,jmold,iim,jjm,1,
      &                   rlonuold,rlatvold,rlonu,rlatv)
-      call gr_dyn_fi (1,iim+1,jjm+1,ngrid,emiss,emis)
+c Dust conversion factor
+      call interp_horiz (tauscalingold,tauscalings,imold,jmold,iim,jjm,
+     &                   1,rlonuold,rlatvold,rlonu,rlatv)
+      call gr_dyn_fi (1,iim+1,jjm+1,ngrid,tauscalings,tauscaling)
 c     write(46,*) 'emis',emis
 c-----------------------------------------------------------------------
@@ -1378 +1402 @@
       deallocate(mlayerold)
       deallocate(qsurfold)
+      deallocate(tauscalingold)
       deallocate(var,varp1)