Changeset 10234


Ignore:
Timestamp:
2018-10-26T16:19:35+02:00 (2 years ago)
Author:
cbricaud
Message:

correc a bug for MPP decompostion for coarsened grid for limited area model

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2018/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/CRS/crsdom.F90

    r10210 r10234  
    189189            ENDDO 
    190190      END SELECT 
    191       WRITE(narea+1000-1,*)"end glam_crs gphi_crs ",p_glam_crs(1,2),p_gphi_crs(1,2) 
    192191 
    193192      ! Retroactively add back the boundary halo cells. 
    194 !????      CALL crs_lbc_lnk( p_gphi_crs, cd_type, 1.0 ) 
    195 !????      CALL crs_lbc_lnk( p_glam_crs, cd_type, 1.0 ) 
     193      !CALL crs_lbc_lnk( p_gphi_crs, cd_type, 1.0 ) 
     194      !CALL crs_lbc_lnk( p_glam_crs, cd_type, 1.0 ) 
    196195      WRITE(narea+1000-1,*)"end1 glam_crs gphi_crs ",p_glam_crs(1,2),p_gphi_crs(1,2) 
    197196      ! 
     
    11351134      SELECT CASE ( jperio ) 
    11361135 
    1137       CASE ( 0, 1 )   ! limited area 
     1136      CASE ( 0, 1 ) 
     1137 
     1138      !exp for a regional configuration 
     1139      ! i=1 => ghost cell   
    11381140 
    11391141              !--!--!--!--!--!--!--!--!--! 
     
    11441146      !IBI36               1  2  3  4  5   
    11451147      !IBI12  !    1   !   2    !   3    !     
    1146  
    1147  
     1148  
    11481149      jpiglo_crs   = INT( (jpiglo - 4) / nn_factx ) + 4 
    11491150      jpjglo_crs   = INT( (jpjglo - 4) / nn_facty ) + 4 
    11501151 
    1151       CASE ( 3, 4 )    ! global  
     1152      jpj_crs = ( jpjglo_crs -1  - 2 * jprecj + (jpnj-1) ) / jpnj + 2 * jprecj 
     1153      IF( MOD( jpjglo_crs   - 2 * jprecj + (jpnj-1) , jpnj ) .NE.  0._wp )THEN 
     1154         jpj_crs = jpj_crs + 1 
     1155      ENDIF 
     1156 
     1157      CASE ( 3, 4 )    !   
    11521158      jpiglo_crs   = INT( (jpiglo - 2) / nn_factx ) + 2 
    11531159      jpjglo_crs   = INT( (jpjglo - MOD(jpjglo, nn_facty)) / nn_facty ) + 2 
     1160 
     1161      jpj_crs = ( jpjglo_crs -1  - 2 * jprecj + (jpnj-1) ) / jpnj + 2 * jprecj 
    11541162 
    11551163      END SELECT 
     
    11591167 
    11601168      jpi_crs = ( jpiglo_crs   - 2 * jpreci + (jpni-1) ) / jpni + 2 * jpreci 
    1161       jpj_crs = ( jpjglo_crsm1 - 2 * jprecj + (jpnj-1) ) / jpnj + 2 * jprecj 
    11621169        
    11631170      jpi_crsm1   = jpi_crs - 1 
     
    13181325                  nldi_crs=2 
    13191326                  mis2_crs(ijis-1) = mie2_crs(ijis-1) -1 
     1327               !CASE DEFAULT 
     1328               !  WRITE(narea+8000-1,*)"WRONG VALUE FOR iistart ",ii_start 
    13201329            END SELECT 
    13211330 
     
    15551564                     nldj_crs=2 
    15561565                     mjs2_crs(ijjs-1) = mje2_crs(ijjs-1) -1 
    1557                   CASE DEFAULT 
    1558                      WRITE(narea+8000-1,*)"WRONG VALUE FOR iistart ",ii_start 
     1566                  !CASE DEFAULT 
     1567                  !   WRITE(narea+8000-1,*)"WRONG VALUE FOR iistart ",ii_start 
    15591568               END SELECT 
    15601569 
     
    18951904 
    18961905      !!---------------------------------------------------------------- 
    1897  
    18981906      mikt_crs(:,:)=1 
    18991907 
Note: See TracChangeset for help on using the changeset viewer.