Changeset 9449


Ignore:
Timestamp:
2018-03-30T17:37:02+02:00 (3 years ago)
Author:
smasson
Message:

dev_merge_2017: agrif bugfix for non-constant jpi/jpj + some cleaning…

Location:
branches/2017/dev_merge_2017/NEMOGCM/NEMO
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/NST_SRC/agrif_user.F90

    r9169 r9449  
    1919   USE dom_oce 
    2020   USE nemogcm 
     21   USE mppini 
    2122   !! 
    2223   IMPLICIT NONE 
     
    2425   ! 
    2526   IF( .NOT. Agrif_Root() ) THEN 
    26       jpni  = Agrif_Parent(jpni) 
    27       jpnj  = Agrif_Parent(jpnj) 
    28       jpnij = Agrif_Parent(jpnij) 
    29       jpiglo = nbcellsx + 2 + 2*nbghostcells 
    30       jpjglo = nbcellsy + 2 + 2*nbghostcells 
    31       jpi    = ( jpiglo-2*nn_hls + (jpni-1+0) ) / jpni + 2*nn_hls 
    32       jpj    = ( jpjglo-2*nn_hls + (jpnj-1+0) ) / jpnj + 2*nn_hls 
    33 ! JC: change to allow for different vertical levels 
    34 !     jpk is already set 
    35 !     keep it jpk possibly different from jpkglo which  
    36 !     hold parent grid vertical levels number (set earlier) 
    37 !      jpk     = jpkglo  
    38       jpim1  = jpi-1  
    39       jpjm1  = jpj-1  
    40       jpkm1  = MAX( 1, jpk-1 )                                          
    41       jpij   = jpi*jpj  
    42       nperio = 0 
    43       jperio = 0 
     27      ! no more static variables 
     28!!$! JC: change to allow for different vertical levels 
     29!!$!     jpk is already set 
     30!!$!     keep it jpk possibly different from jpkglo which  
     31!!$!     hold parent grid vertical levels number (set earlier) 
     32!!$!      jpk     = jpkglo  
    4433   ENDIF 
    4534   ! 
     
    6352   !!---------------------------------------------------------------------- 
    6453   ! 
    65 !!gm  I think this is now useless ...   nn_cfg & cn_cfg are set to -999999 and "UNKNOWN"  
    66 !!gm                                    when reading the AGRIF domain configuration file 
    67    IF( cn_cfg == 'orca' ) THEN 
    68       IF ( nn_cfg == 2 .OR. nn_cfg == 025 .OR. nn_cfg == 05  .OR. nn_cfg == 4 ) THEN 
    69          nn_cfg = -1    ! set special value for nn_cfg on fine grids 
    70          cn_cfg = "default" 
    71       ENDIF 
    72    ENDIF 
    73    !                    !* Specific fine grid Initializations 
    74    ln_tradmp = .FALSE.        ! no tracer damping on fine grids 
    75    ! 
    76    ln_bdy    = .FALSE.        ! no open boundary on fine grids 
    77  
    7854   CALL nemo_init       !* Initializations of each fine grid 
    7955 
  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/BDY/bdyini.F90

    r9169 r9449  
    8585      IF(lwm) WRITE ( numond, nambdy ) 
    8686 
     87      IF( .NOT. Agrif_Root() ) ln_bdy = .FALSE.   ! forced for Agrif children 
     88       
    8789      ! ----------------------------------------- 
    8890      ! unstructured open boundaries use control 
  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/DOM/domvvl.F90

    r9404 r9449  
    147147      CALL dom_vvl_interpol( e3v_n(:,:,:), e3vw_n(:,:,:), 'VW' )  ! from V to UW 
    148148      CALL dom_vvl_interpol( e3v_b(:,:,:), e3vw_b(:,:,:), 'VW' ) 
     149 
     150      ! We need to define e3[tuv]_a for AGRIF initialisation (should not be a problem for the restartability...) 
     151      e3t_a(:,:,:) = e3t_n(:,:,:) 
     152      e3u_a(:,:,:) = e3u_n(:,:,:) 
     153      e3v_a(:,:,:) = e3v_n(:,:,:) 
    149154      ! 
    150155      !                    !==  depth of t and w-point  ==!   (set the isf depth as it is in the initial timestep) 
  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/LBC/mppini.F90

    r9446 r9449  
    106106   !!---------------------------------------------------------------------- 
    107107 
     108 
    108109   SUBROUTINE mpp_init 
    109110      !!---------------------------------------------------------------------- 
     
    159160      ! If dimensions of processor grid weren't specified in the namelist file 
    160161      ! then we calculate them here now that we have our communicator size 
    161       IF( jpni < 1 .OR. jpnj < 1 ) THEN 
    162          IF( Agrif_Root() )   CALL mpp_init_partition( mppsize ) 
    163       ENDIF 
     162      IF( jpni < 1 .OR. jpnj < 1 )   CALL mpp_init_partition( mppsize ) 
    164163      ! 
    165164#if defined key_agrif 
    166165      IF( jpnij /= jpni*jpnj ) CALL ctl_stop( 'STOP', 'Cannot remove land proc with AGRIF' ) 
    167166#endif 
    168  
    169167      ! 
    170168      ALLOCATE(  nfiimpp(jpni,jpnj), nfipproc(jpni,jpnj), nfilcit(jpni,jpnj) ,    & 
     
    183181       
    184182      ! 
    185 #if defined key_agrif 
    186183      IF( .NOT. Agrif_Root() ) THEN       ! AGRIF children: specific setting (cf. agrif_user.F90) 
    187          jpiglo  = nbcellsx + 2 + 2*nbghostcells 
    188          jpjglo  = nbcellsy + 2 + 2*nbghostcells 
    189          jpimax  = ( jpiglo-2*nn_hls + (jpni-1+0) ) / jpni + 2*nn_hls 
    190          jpjmax  = ( jpjglo-2*nn_hls + (jpnj-1+0) ) / jpnj + 2*nn_hls 
    191          jperio  = 0 
    192          ln_use_jattr = .false. 
    193       ENDIF 
     184         IF( jpiglo /= nbcellsx + 2 + 2*nbghostcells )   & 
     185            CALL ctl_stop( 'STOP', 'mpp_init: Agrif children requires jpiglo == nbcellsx + 2 + 2*nbghostcells' ) 
     186         IF( jpjglo /= nbcellsy + 2 + 2*nbghostcells )   & 
     187            CALL ctl_stop( 'STOP', 'mpp_init: Agrif children requires jpjglo == nbcellsy + 2 + 2*nbghostcells' ) 
     188         IF( ln_use_jattr )   CALL ctl_stop( 'STOP', 'mpp_init:Agrif children requires ln_use_jattr = .false. ' ) 
     189      ENDIF 
     190 
     191#if defined key_nemocice_decomp 
     192      jpimax = ( nx_global+2-2*nn_hls + (jpni-1) ) / jpni + 2*nn_hls    ! first  dim. 
     193      jpjmax = ( ny_global+2-2*nn_hls + (jpnj-1) ) / jpnj + 2*nn_hls    ! second dim.  
     194#else 
     195      jpimax = ( jpiglo - 2*nn_hls + (jpni-1) ) / jpni + 2*nn_hls    ! first  dim. 
     196      jpjmax = ( jpjglo - 2*nn_hls + (jpnj-1) ) / jpnj + 2*nn_hls    ! second dim. 
    194197#endif 
    195  
    196       IF( Agrif_Root() ) THEN       ! AGRIF mother: specific setting from jpni and jpnj 
    197 #if defined key_nemocice_decomp 
    198          jpimax = ( nx_global+2-2*nn_hls + (jpni-1) ) / jpni + 2*nn_hls    ! first  dim. 
    199          jpjmax = ( ny_global+2-2*nn_hls + (jpnj-1) ) / jpnj + 2*nn_hls    ! second dim.  
    200 #else 
    201          jpimax = ( jpiglo - 2*nn_hls + (jpni-1) ) / jpni + 2*nn_hls    ! first  dim. 
    202          jpjmax = ( jpjglo - 2*nn_hls + (jpnj-1) ) / jpnj + 2*nn_hls    ! second dim. 
    203 #endif 
    204       ENDIF 
    205198 
    206199      ! 
     
    449442      ii_noea(:) = -1 
    450443      ii_nowe(:) = -1  
    451       DO jarea = 1, jpnij 
    452          ii = iin(jarea) 
    453          ij = ijn(jarea) 
     444      DO jproc = 1, jpnij 
     445         ii = iin(jproc) 
     446         ij = ijn(jproc) 
    454447         IF( 0 <= ioso(ii,ij) .AND. ioso(ii,ij) <= (jpni*jpnj-1) ) THEN 
    455448            iiso = 1 + MOD( ioso(ii,ij) , jpni ) 
    456449            ijso = 1 +      ioso(ii,ij) / jpni 
    457             ii_noso(jarea) = ipproc(iiso,ijso) 
     450            ii_noso(jproc) = ipproc(iiso,ijso) 
    458451         ENDIF 
    459452         IF( 0 <= iowe(ii,ij) .AND. iowe(ii,ij) <= (jpni*jpnj-1) ) THEN 
    460453          iiwe = 1 + MOD( iowe(ii,ij) , jpni ) 
    461454          ijwe = 1 +      iowe(ii,ij) / jpni 
    462           ii_nowe(jarea) = ipproc(iiwe,ijwe) 
     455          ii_nowe(jproc) = ipproc(iiwe,ijwe) 
    463456         ENDIF 
    464457         IF( 0 <= ioea(ii,ij) .AND. ioea(ii,ij) <= (jpni*jpnj-1) ) THEN 
    465458            iiea = 1 + MOD( ioea(ii,ij) , jpni ) 
    466459            ijea = 1 +      ioea(ii,ij) / jpni 
    467             ii_noea(jarea)= ipproc(iiea,ijea) 
     460            ii_noea(jproc)= ipproc(iiea,ijea) 
    468461         ENDIF 
    469462         IF( 0 <= iono(ii,ij) .AND. iono(ii,ij) <= (jpni*jpnj-1) ) THEN 
    470463            iino = 1 + MOD( iono(ii,ij) , jpni ) 
    471464            ijno = 1 +      iono(ii,ij) / jpni 
    472             ii_nono(jarea)= ipproc(iino,ijno) 
     465            ii_nono(jproc)= ipproc(iino,ijno) 
    473466         ENDIF 
    474467      END DO 
     
    501494      ! Save maximum number of levels in jpkglo, then define all vertical grids with this number. 
    502495      ! Suppress once vertical online interpolation is ok 
    503       IF(.NOT.Agrif_Root())   jpkglo = Agrif_Parent( jpkglo ) 
     496!!$      IF(.NOT.Agrif_Root())   jpkglo = Agrif_Parent( jpkglo ) 
    504497#endif 
    505498      jpim1 = jpi-1                                            ! inner domain indices 
  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/par_oce.F90

    r9019 r9449  
    4141   INTEGER       ::   jpkglo           !: 3nd    -                  -    --> k levels 
    4242 
    43 #if defined key_agrif 
    44  
    45 !!gm  BUG ?   I'm surprised by the calculation below of nbcellsx and nbcellsy before jpiglo,jpjglo  
    46 !!gm                           has been assigned to a value.... 
    47 !!gm 
    48  
    4943   ! global domain size for AGRIF     !!! * total AGRIF computational domain * 
     44   INTEGER, PUBLIC            ::   nbug_in_agrif_conv_do_not_remove_or_modify = 1 - 1 
    5045   INTEGER, PUBLIC, PARAMETER ::   nbghostcells = 1                             !: number of ghost cells 
    51    INTEGER, PUBLIC            ::   nbcellsx     = jpiglo - 2 - 2*nbghostcells   !: number of cells in i-direction 
    52    INTEGER, PUBLIC            ::   nbcellsy     = jpjglo - 2 - 2*nbghostcells   !: number of cells in j-direction 
    53 #endif 
     46   INTEGER, PUBLIC            ::   nbcellsx   ! = jpiglo - 2 - 2*nbghostcells   !: number of cells in i-direction 
     47   INTEGER, PUBLIC            ::   nbcellsy   ! = jpjglo - 2 - 2*nbghostcells   !: number of cells in j-direction 
    5448 
    5549   ! local domain size                !!! * local computational domain * 
    56    INTEGER, PUBLIC ::   jpi   ! = ( jpiglo-2*nn_hls + (jpni-1) ) / jpni + 2*nn_hls   !: first  dimension 
    57    INTEGER, PUBLIC ::   jpj   ! = ( jpjglo-2*nn_hls + (jpnj-1) ) / jpnj + 2*nn_hls   !: second dimension 
    58    INTEGER, PUBLIC ::   jpk   ! = jpkglo 
     50   INTEGER, PUBLIC ::   jpi   !                                                    !: first  dimension 
     51   INTEGER, PUBLIC ::   jpj   !                                                    !: second dimension 
     52   INTEGER, PUBLIC ::   jpk   ! = jpkglo                                           !: third  dimension 
    5953   INTEGER, PUBLIC ::   jpim1 ! = jpi-1                                            !: inner domain indices 
    6054   INTEGER, PUBLIC ::   jpjm1 ! = jpj-1                                            !:   -     -      - 
    6155   INTEGER, PUBLIC ::   jpkm1 ! = jpk-1                                            !:   -     -      - 
    6256   INTEGER, PUBLIC ::   jpij  ! = jpi*jpj                                          !:  jpi x jpj 
    63    INTEGER, PUBLIC ::   jpimax! = maximum jpi across all areas  
    64    INTEGER, PUBLIC ::   jpjmax! = maximum jpj across all areas 
     57   INTEGER, PUBLIC ::   jpimax! = ( jpiglo-2*nn_hls + (jpni-1) ) / jpni + 2*nn_hls !: maximum jpi 
     58   INTEGER, PUBLIC ::   jpjmax! = ( jpjglo-2*nn_hls + (jpnj-1) ) / jpnj + 2*nn_hls !: maximum jpj 
    6559 
    6660   !!--------------------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.