Changeset 5755


Ignore:
Timestamp:
2015-09-22T15:42:08+02:00 (5 years ago)
Author:
flavoni
Message:

remove hard coded reduction of scale factors

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/NEMO/OPA_SRC/DOM/domhgr.F90

    r5737 r5755  
    144144         ! 
    145145         !                                                ! ===================== 
    146          IF( cp_cfg == "orca" .AND. jp_cfg == 2 ) THEN    ! ORCA R2 configuration 
    147             !                                             ! ===================== 
    148             IF( nn_cla == 0 ) THEN 
    149                ! 
    150                ii0 = 139   ;   ii1 = 140        ! Gibraltar Strait (e2u = 20 km) 
    151                ij0 = 102   ;   ij1 = 102   ;   e2u( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  20.e3 
    152                IF(lwp) WRITE(numout,*) 
    153                IF(lwp) WRITE(numout,*) '             orca_r2: Gibraltar    : e2u reduced to 20 km' 
    154                ! 
    155                ii0 = 160   ;   ii1 = 160        ! Bab el Mandeb (e2u = 18 km) 
    156                ij0 =  88   ;   ij1 =  88   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  18.e3 
    157                                                e2u( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  30.e3 
    158                IF(lwp) WRITE(numout,*) 
    159                IF(lwp) WRITE(numout,*) '             orca_r2: Bab el Mandeb: e2u reduced to 30 km' 
    160                IF(lwp) WRITE(numout,*) '                                     e1v reduced to 18 km' 
    161             ENDIF 
    162  
    163             ii0 = 145   ;   ii1 = 146        ! Danish Straits (e2u = 10 km) 
    164             ij0 = 116   ;   ij1 = 116   ;   e2u( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  10.e3 
    165             IF(lwp) WRITE(numout,*) 
    166             IF(lwp) WRITE(numout,*) '             orca_r2: Danish Straits : e2u reduced to 10 km' 
    167             ! 
    168          ENDIF 
    169          ! 
    170          !                                                ! ===================== 
    171          IF( cp_cfg == "orca" .AND. jp_cfg == 1 ) THEN    ! ORCA R1 configuration 
    172             !                                             ! ===================== 
    173             ! This dirty section will be suppressed by simplification process: all this will come back in input files 
    174             ! Currently these hard-wired indices relate to configuration with 
    175             ! extend grid (jpjglo=332)  which had a grid-size of 362x292. 
    176             !  
    177             isrow = 332 - jpjglo 
    178             ! 
    179             ii0 = 282           ;   ii1 = 283        ! Gibraltar Strait (e2u = 20 km) 
    180             ij0 = 241 - isrow   ;   ij1 = 241 - isrow   ;   e2u( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  20.e3 
    181             IF(lwp) WRITE(numout,*) 
    182             IF(lwp) WRITE(numout,*) '             orca_r1: Gibraltar : e2u reduced to 20 km' 
    183  
    184             ii0 = 314           ;   ii1 = 315        ! Bhosporus Strait (e2u = 10 km) 
    185             ij0 = 248 - isrow   ;   ij1 = 248 - isrow   ;   e2u( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  10.e3 
    186             IF(lwp) WRITE(numout,*) 
    187             IF(lwp) WRITE(numout,*) '             orca_r1: Bhosporus : e2u reduced to 10 km' 
    188  
    189             ii0 =  44           ;   ii1 =  44        ! Lombok Strait (e1v = 13 km) 
    190             ij0 = 164 - isrow   ;   ij1 = 165 - isrow   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  13.e3 
    191             IF(lwp) WRITE(numout,*) 
    192             IF(lwp) WRITE(numout,*) '             orca_r1: Lombok : e1v reduced to 10 km' 
    193  
    194             ii0 =  48           ;   ii1 =  48        ! Sumba Strait (e1v = 8 km) [closed from bathy_11 on] 
    195             ij0 = 164 - isrow   ;   ij1 = 165 - isrow   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  8.e3 
    196             IF(lwp) WRITE(numout,*) 
    197             IF(lwp) WRITE(numout,*) '             orca_r1: Sumba : e1v reduced to 8 km' 
    198  
    199             ii0 =  53           ;   ii1 =  53        ! Ombai Strait (e1v = 13 km) 
    200             ij0 = 164 - isrow   ;   ij1 = 165 - isrow   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 13.e3 
    201             IF(lwp) WRITE(numout,*) 
    202             IF(lwp) WRITE(numout,*) '             orca_r1: Ombai : e1v reduced to 13 km' 
    203  
    204             ii0 =  56           ;   ii1 =  56        ! Timor Passage (e1v = 20 km) 
    205             ij0 = 164 - isrow   ;   ij1 = 145 - isrow   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 20.e3 
    206             IF(lwp) WRITE(numout,*) 
    207             IF(lwp) WRITE(numout,*) '             orca_r1: Timor Passage : e1v reduced to 20 km' 
    208  
    209             ii0 =  55           ;   ii1 =  55        ! West Halmahera Strait (e1v = 30 km) 
    210             ij0 = 181 - isrow   ;   ij1 = 182 - isrow   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 30.e3 
    211             IF(lwp) WRITE(numout,*) 
    212             IF(lwp) WRITE(numout,*) '             orca_r1: W Halmahera : e1v reduced to 30 km' 
    213  
    214             ii0 =  58           ;   ii1 =  58        ! East Halmahera Strait (e1v = 50 km) 
    215             ij0 = 181 - isrow   ;   ij1 = 182 - isrow   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 50.e3 
    216             IF(lwp) WRITE(numout,*) 
    217             IF(lwp) WRITE(numout,*) '             orca_r1: E Halmahera : e1v reduced to 50 km' 
    218             ! 
    219          ENDIF 
    220          ! 
    221          !                                                ! ====================== 
    222          IF( cp_cfg == "orca" .AND. jp_cfg == 05 ) THEN   ! ORCA R05 configuration 
    223             !                                             ! ====================== 
    224             ii0 = 563   ;   ii1 = 564        ! Gibraltar Strait (e2u = 20 km) 
    225             ij0 = 327   ;   ij1 = 327   ;   e2u( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  20.e3 
    226             IF(lwp) WRITE(numout,*) 
    227             IF(lwp) WRITE(numout,*) '             orca_r05: Reduced e2u at the Gibraltar Strait' 
    228             ! 
    229             ii0 = 627   ;   ii1 = 628        ! Bosphore Strait (e2u = 10 km) 
    230             ij0 = 343   ;   ij1 = 343   ;   e2u( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  10.e3 
    231             IF(lwp) WRITE(numout,*) 
    232             IF(lwp) WRITE(numout,*) '             orca_r05: Reduced e2u at the Bosphore Strait' 
    233             ! 
    234             ii0 =  93   ;   ii1 =  94        ! Sumba Strait (e2u = 40 km) 
    235             ij0 = 232   ;   ij1 = 232   ;   e2u( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  40.e3 
    236             IF(lwp) WRITE(numout,*) 
    237             IF(lwp) WRITE(numout,*) '             orca_r05: Reduced e2u at the Sumba Strait' 
    238             ! 
    239             ii0 = 103   ;   ii1 = 103        ! Ombai Strait (e2u = 15 km) 
    240             ij0 = 232   ;   ij1 = 232   ;   e2u( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  15.e3 
    241             IF(lwp) WRITE(numout,*) 
    242             IF(lwp) WRITE(numout,*) '             orca_r05: Reduced e2u at the Ombai Strait' 
    243             ! 
    244             ii0 =  15   ;   ii1 =  15        ! Palk Strait (e2u = 10 km) 
    245             ij0 = 270   ;   ij1 = 270   ;   e2u( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  10.e3 
    246             IF(lwp) WRITE(numout,*) 
    247             IF(lwp) WRITE(numout,*) '             orca_r05: Reduced e2u at the Palk Strait' 
    248             ! 
    249             ii0 =  87   ;   ii1 =  87        ! Lombok Strait (e1v = 10 km) 
    250             ij0 = 232   ;   ij1 = 233   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  10.e3 
    251             IF(lwp) WRITE(numout,*) 
    252             IF(lwp) WRITE(numout,*) '             orca_r05: Reduced e1v at the Lombok Strait' 
    253             ! 
    254             ! 
    255             ii0 = 662   ;   ii1 = 662        ! Bab el Mandeb (e1v = 25 km) 
    256             ij0 = 276   ;   ij1 = 276   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  25.e3 
    257             IF(lwp) WRITE(numout,*) 
    258             IF(lwp) WRITE(numout,*) '             orca_r05: Reduced e1v at the Bab el Mandeb' 
    259             ! 
    260          ENDIF 
    261           
    262             !                       ! create 'coordinate_e1e2u_v.nc' file  that contains 
    263             !                       ! reduced scale factor in some strait but full e1e2u and e1e2v surfaces          
    264146            IF( ie1e2u_v == 0 )   CALL dom_wri_coordinate 
    265147            ! 
Note: See TracChangeset for help on using the changeset viewer.