Changeset 1892 for trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan/lect_start_archive.F

Timestamp:

Jan 10, 2018, 6:21:35 PM (7 years ago)

Author:

jvatant

Message:

Making chemistry handling more flexible - step 3
+ Added the horizontal interpolation for chemistry fields
+ Next step : deal with vertical interp !
--JVO

File:

• trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan/lect_start_archive.F (modified) (6 diffs)

trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan/lect_start_archive.F

@@ -10 +10 @@
 !      USE control_mod
 ! to use  'getin'
+      USE callkeys_mod, only: callchim
       USE comvert_mod, ONLY: ap,bp,aps,bps,preff
       USE comconst_mod, ONLY: kappa,g,pi
@@ -256 +257 @@
       ierr= NF_INQ_DIMID(nid,"upper_chemistry_layers",dimid)
       ierr= NF_INQ_DIMLEN(nid,dimid,nlaykimold)
-
-      ! NB : The vertical regriding, if needed cannot be done here since the new
-      ! pressure grid is only computed at the end of newstart
-      ! Here we will just do the horizontal interpolation on scalar grid.
 
 ! 1.3 Report dimensions
@@ -307 +304 @@
       
       allocate(preskimold(nlaykimold))
-      allocate(ykim_upS(44,iip1,jjp1))
-      allocate(ykim_upoldS(44,imold+1,jmold+1))
+      allocate(ykim_upS(44,iip1,jjp1,nlaykimold))
+      allocate(ykim_upoldS(44,imold+1,jmold+1,nlaykimold))
+      allocate(ykim_up_oldv(44,ngrid,nlaykimold))
 
       allocate(var (imold+1,jmold+1,llm))
@@ -830 +828 @@
 !      endif
 
-c-----------------------------------------------------------------------
-c 5.3   Lecture des champs 3D (t,u,v, q2atm,q)
+
+c-----------------------------------------------------------------------
+c 5.3 Read 3D upper chemistry fields, if needed
+c-----------------------------------------------------------------------  
+       
+      start=(/1,1,1,memo/)
+      count=(/imold+1,jmold+1,nlaykimold,1/)
+    
+c NB : The sanity check on callchim is on H_up but could be on any chem. field
+c as we assume we can't do incomplete chemistry - JVO 18
+    
+      PRINT*, "lect_start_archive: loading upper chemistry fields..."
+      
+      ierr=NF_INQ_VARID(nid,"H_up",nvarid)
+      
+      IF (ierr .NE. NF_NOERR) THEN
+      
+         PRINT*, "lect_start_archive: Le champ <H_up> est absent..."        
+         IF (callchim) THEN
+           PRINT*, "... mais callchim=.TRUE. dans callphys.def !"
+           PRINT*, "Verifiez start_archive.nc ou desactivez callchim !"
+           CALL abort
+         ELSE
+           PRINT*, '... je suppose que les autres champs aussi et je
+     &passerai donc mon chemin pour tout ce qui concerne la chimie !'
+           WRITE (*,*)
+         ENDIF
+         
+      ELSE
+      
+        IF (.not.callchim) THEN
+          PRINT*, "lect_start_archive: Le champ <H_up> est present..."
+          PRINT*, "... mais callchim=.FALSE. dans callphys.def !"
+          PRINT*, "Veuillez activer callchim pour gerer ce champ !"
+          CALL abort
+        ELSE
+#ifdef NC_DOUBLE
+          ierr = NF_GET_VARA_DOUBLE(nid,nvarid,start,count,
+     &                             ykim_upoldS(1,:,:,:))
+#else
+          ierr = NF_GET_VARA_REAL(nid,nvarid,start,count,
+     &                            ykim_upoldS(1,:,:,:))
+#endif
+          IF (ierr .NE. NF_NOERR) THEN
+             PRINT*, "lect_start_archive: Lecture echouee pour <H_up>"
+             CALL abort
+          ENDIF          
+          ! Then read all the others by their name if needed
+          CALL read_startarch_kim(nid,start,count)              
+         ENDIF
+         
+      ENDIF  ! if ierr.ne.nf_no_err
+
+c-----------------------------------------------------------------------
+c 5.4   Lecture des champs 3D (t,u,v, q2atm,q)
 c-----------------------------------------------------------------------
 
@@ -1161 +1212 @@
       ! Reshape tsoilS to scalar grid as tsoil
        call gr_dyn_fi (nsoilmx,iim+1,jjm+1,ngrid,tsoilS,tsoil)
-
-c-----------------------------------------------------------------------
-c 6.3 Variable 3d :
+       
+c-----------------------------------------------------------------------
+c 6.3 Upper chemistry 3d fields:
+c
+c NB : The vertical regriding, if needed cannot be done here since the
+c new pressure grid is only computed at the end of newstart.
+c Here we will just do the horizontal interpolation on scalar grid.
+c
+c-----------------------------------------------------------------------
+
+      ! Do the horizontal interpolation
+      DO i=1,44
+        call interp_horiz(ykim_upoldS(i,:,:,:),ykim_upS(i,:,:,:),
+     &                  imold,jmold,iim,jjm,nlaykimold,
+     &                   rlonuold,rlatvold,rlonu,rlatv)
+
+        ! Reshape ykim_upS to scalar grid as ykim_up_oldv
+        call gr_dyn_fi(nlaykimold,iim+1,jjm+1,ngrid,
+     &                 ykim_upS(i,:,:,:),ykim_up_oldv(i,:,:))
+      ENDDO
+
+
+c-----------------------------------------------------------------------
+c 6.4 Variable 3d :
 c-----------------------------------------------------------------------
       
@@ -1328 +1400 @@
       deallocate(qsurfold)
       deallocate(var,varp1)
+      
+      deallocate(ykim_upS)
+      deallocate(ykim_upoldS)
 
 !      write(*,*)'lect_start_archive: END'