Changeset 4694 for branches/2014/dev_r4621_NOC4_BDY_VERT_INTERP/NEMOGCM/NEMO/OPA_SRC/BDY/bdyini.F90

Timestamp:

2014-06-26T12:22:40+02:00 (10 years ago)

Author:

jamesharle

Message:

Update of fldread to handle depth information in BDY files and addition of an interpolation routine. Updated BDY code to handle T/S BDY interpolation on the fly. Conservative remapping of U/V still to be coded. Not compiled or test yet.

File:

• branches/2014/dev_r4621_NOC4_BDY_VERT_INTERP/NEMOGCM/NEMO/OPA_SRC/BDY/bdyini.F90 (modified) (3 diffs)

branches/2014/dev_r4621_NOC4_BDY_VERT_INTERP/NEMOGCM/NEMO/OPA_SRC/BDY/bdyini.F90

@@ -103 +103 @@
          &             cn_ice_lim, nn_ice_lim_dta,                           &
 #endif
-         &             ln_vol, nn_volctl, nn_rimwidth
+         &             ln_vol, nn_volctl, nn_rimwidth, nb_jpk_bdy
       !!
       NAMELIST/nambdy_index/ ctypebdy, nbdyind, nbdybeg, nbdyend
@@ -380 +380 @@
           IF(lwp) WRITE(numout,*)
         ENDIF
+        IF( nb_jpk_bdy > 0 ) THEN
+           IF(lwp) WRITE(numout,*) '*** open boundary will be interpolate in the vertical onto the native grid ***'
+        ELSE
+           IF(lwp) WRITE(numout,*) '*** open boundary will be read straight onto the native grid without vertical interpolation ***'
      ENDIF
 
@@ -505 +509 @@
             &      nbrdta(jpbdta, jpbgrd, nb_bdy) )
 
-         ALLOCATE( dta_global(jpbdtau, 1, jpk) )
-         IF ( icount>0 ) ALLOCATE( dta_global2(jpbdtas, nrimmax, jpk) )
+         IF( jpk_bdy>0 ) THEN
+            ALLOCATE( dta_global(jpbdtau, 1, jpk_bdy) )
+            ALLOCATE( dta_global_z(jpbdtau, 1, jpk_bdy) )
+         ELSE
+            ALLOCATE( dta_global(jpbdtau, 1, jpk) )
+            ALLOCATE( dta_global_z(jpbdtau, 1, jpk) )
+         ENDIF
+
+         IF ( icount>0 ) THEN
+            IF( jpk_bdy>0 ) THEN
+               ALLOCATE( dta_global2(jpbdtas, nrimmax, jpk_bdy) )
+               ALLOCATE( dta_global2_z(jpbdtas, nrimmax, jpk_bdy) )
+            ELSE
+               ALLOCATE( dta_global2(jpbdtas, nrimmax, jpk) )
+               ALLOCATE( dta_global2_z(jpbdtas, nrimmax, jpk) )
+            ENDIF
          ! 
       ENDIF