Changeset 10239


Ignore:
Timestamp:
2018-10-26T17:43:47+02:00 (19 months ago)
Author:
cmao
Message:

Copy in changes to dev_r5003_MERCATOR6_CRS at r10234

Location:
branches/UKMO/dev_r10171_test_crs_AMM7/NEMOGCM/NEMO/OPA_SRC/CRS
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r10171_test_crs_AMM7/NEMOGCM/NEMO/OPA_SRC/CRS/crs.F90

    r10207 r10239  
    326326#elif defined key_traldf_c1d 
    327327      ALLOCATE( ahtt_crs(        jpk) , ahtu_crs(        jpk) , ahtv_crs(        jpk) , ahtw_crs(        jpk) , & 
     328#else 
     329      ALLOCATE( & 
    328330#endif 
    329331              & r_fact_lap_crs(jpi_crs,jpj_crs,jpk) , STAT=ierr(13) ) 
  • branches/UKMO/dev_r10171_test_crs_AMM7/NEMOGCM/NEMO/OPA_SRC/CRS/crsdom.F90

    r10216 r10239  
    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              !--!--!--!--!--!--!--!--!--! 
     
    11431145              !--!--!--!--!--!--!--!--!--! 
    11441146      !IBI36               1  2  3  4  5   
    1145       !IBI12  !    1   !   2    !   3    !  
    1146        
    1147  
     1147      !IBI12  !    1   !   2    !   3    !     
     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 
    13221331         ENDIF 
    1323  
     1332          
    13241333         !---------------------------------------------------------------------------------------------- 
    13251334         ! I-4 compute nimpp_crs 
     
    13681377            mie_crs(ji) = mie2_crs(mig_crs(ji)) - nimpp + 1 
    13691378            nfactx(ji)  = mie_crs(ji)-mis_crs(ji)+1 
    1370          ENDDO          
     1379         ENDDO 
    13711380 
    13721381         !----------------------------------------------------------------------------- 
     
    14561465         END SELECT 
    14571466 
    1458  
     1467          
    14591468         !WRITE(narea+1000-1,*)"loc crs jpi nldi,nlei,nlci ",jpi_crs, nldi_crs            ,nlei_crs             ,nlci_crs 
    14601469         !CALL FLUSH(narea+1000-1) 
     
    15341543         !     global: special case when the domain on the north is land:  
    15351544         !             compute mjs2_crs and mje2_crs for the first cell of the local domain  
    1536          !     local :  special case for southern domain   
     1545         !     local :  special case for southern domain  
    15371546         !--------------------------------------------------------------------------------------------- 
    15381547         nldj_crs = 2 
     
    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.