Ignore:
Timestamp:
04/14/15 17:41:48 (9 years ago)
Author:
omamce
Message:

O.M. : handling old ORCA1 and extended ORCA1 grids

in closea.F90

File:
1 edited

Legend:

Unmodified
Added
Removed
  • CONFIG/UNIFORM/v6/IPSLCM6/SOURCES/NEMO/closea.F90

    r2223 r2485  
    7777      !!---------------------------------------------------------------------- 
    7878      INTEGER ::   jc               ! dummy loop indices 
    79       REAL(wp)::   ztmp    
     79      REAL(wp)::   ztmp 
     80      INTEGER ::   isrow 
    8081      !!---------------------------------------------------------------------- 
    8182       
     
    9697         !                                           ! ======================= 
    9798         CASE ( 1 )                                  ! ORCA_R1 configuration 
    98             !                                        ! ======================= 
     99            ! 
     100            !! This dirty section will be suppressed by simplification process: all this will come back in input files 
     101            !! Currently these hard-wired indices relate to the original (pre-v3.6) configuration which had a grid-size of 362x292. 
     102            !! This grid has been extended southwards for use with the under ice-shelf options (isf) introduced in v3.6. The original 
     103            !! domain can still be used optionally if the isf code is not activated. An adjustment (isrow) is made to the hard-wired 
     104            !! indices if the extended domain (362x332) is being used. 
     105            !! ======================= 
     106            IF     ( jpjglo == 292 ) THEN  !  ORCA1 R1 ­ Using pre-v3.6 files or adjusted start row from isf-extended grid 
     107               isrow = 0 
     108            ELSEIF ( jpjglo == 332 ) THEN  !  ORCA1 R1 - Using full isf­extended domain. 
     109               isrow = 40                  !    - Adjust j­indices to account for more southerly starting latitude 
     110            ENDIF 
     111             
    99112            ncsnr(1)   = 1    ; ncstt(1)   = 0           ! Caspian Sea 
    100             ncsi1(1)   = 332  ; ncsj1(1)   = 203           ! spread over the globe 
    101             ncsi2(1)   = 344  ; ncsj2(1)   = 235 
     113            ncsi1(1)   = 332  ; ncsj1(1)   = 203 + isrow   ! spread over the globe 
     114            ncsi2(1)   = 344  ; ncsj2(1)   = 235 + isrow 
    102115            ncsir(1,1) = 1    ; ncsjr(1,1) = 1 
    103116            !  
    104117            ncsnr(2)   = 4    ; ncstt(2)   = 2           ! Great North American Lakes 
    105             ncsi1(2)   = 198  ; ncsj1(2)   = 209           ! put at St Laurent mouth 
    106             ncsi2(2)   = 213  ; ncsj2(2)   = 223 
    107             ncsir(1,2) = 225  ; ncsjr(1,2) = 220 
    108             ncsir(1,2) = 225  ; ncsjr(1,2) = 221 
    109             ncsir(1,2) = 226  ; ncsjr(1,2) = 220 
    110             ncsir(1,2) = 226  ; ncsjr(1,2) = 221 
     118            ncsi1(2)   = 198  ; ncsj1(2)   = 209 + isrow   ! put at St Laurent mouth 
     119            ncsi2(2)   = 213  ; ncsj2(2)   = 223 + isrow 
     120            ncsir(1,2) = 225  ; ncsjr(1,2) = 220 + isrow 
     121            ncsir(1,2) = 225  ; ncsjr(1,2) = 221 + isrow 
     122            ncsir(1,2) = 226  ; ncsjr(1,2) = 220 + isrow 
     123            ncsir(1,2) = 226  ; ncsjr(1,2) = 221 + isrow 
    111124 
    112125            !                                         
Note: See TracChangeset for help on using the changeset viewer.