Ignore:
Timestamp:
2014-06-26T12:22:40+02:00 (7 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:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2014/dev_r4621_NOC4_BDY_VERT_INTERP/NEMOGCM/NEMO/OPA_SRC/BDY/bdyini.F90

    r4354 r4694  
    103103         &             cn_ice_lim, nn_ice_lim_dta,                           & 
    104104#endif 
    105          &             ln_vol, nn_volctl, nn_rimwidth 
     105         &             ln_vol, nn_volctl, nn_rimwidth, nb_jpk_bdy 
    106106      !! 
    107107      NAMELIST/nambdy_index/ ctypebdy, nbdyind, nbdybeg, nbdyend 
     
    380380          IF(lwp) WRITE(numout,*) 
    381381        ENDIF 
     382        IF( nb_jpk_bdy > 0 ) THEN 
     383           IF(lwp) WRITE(numout,*) '*** open boundary will be interpolate in the vertical onto the native grid ***' 
     384        ELSE 
     385           IF(lwp) WRITE(numout,*) '*** open boundary will be read straight onto the native grid without vertical interpolation ***' 
    382386     ENDIF 
    383387 
     
    505509            &      nbrdta(jpbdta, jpbgrd, nb_bdy) ) 
    506510 
    507          ALLOCATE( dta_global(jpbdtau, 1, jpk) ) 
    508          IF ( icount>0 ) ALLOCATE( dta_global2(jpbdtas, nrimmax, jpk) ) 
     511         IF( jpk_bdy>0 ) THEN 
     512            ALLOCATE( dta_global(jpbdtau, 1, jpk_bdy) ) 
     513            ALLOCATE( dta_global_z(jpbdtau, 1, jpk_bdy) ) 
     514         ELSE 
     515            ALLOCATE( dta_global(jpbdtau, 1, jpk) ) 
     516            ALLOCATE( dta_global_z(jpbdtau, 1, jpk) ) 
     517         ENDIF 
     518 
     519         IF ( icount>0 ) THEN 
     520            IF( jpk_bdy>0 ) THEN 
     521               ALLOCATE( dta_global2(jpbdtas, nrimmax, jpk_bdy) ) 
     522               ALLOCATE( dta_global2_z(jpbdtas, nrimmax, jpk_bdy) ) 
     523            ELSE 
     524               ALLOCATE( dta_global2(jpbdtas, nrimmax, jpk) ) 
     525               ALLOCATE( dta_global2_z(jpbdtas, nrimmax, jpk) ) 
     526            ENDIF 
    509527         !  
    510528      ENDIF 
Note: See TracChangeset for help on using the changeset viewer.