Changeset 6107 for TOOLS/MOSAIX
- Timestamp:
- 03/24/22 09:44:44 (2 years ago)
- Location:
- TOOLS/MOSAIX
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
TOOLS/MOSAIX/CreateOasisGrids.bash
r6064 r6107 66 66 MpiRun="time ccc_mprun" 67 67 PyRun="time ccc_mprun -n 1" # Needed to force python to run on one process only 68 #module purge 69 source $(ccc_home -u igcmg)/MachineEnvironment/irene/env_irene 68 module purge 70 69 source ${SUBMIT_DIR}/arch.env 71 #module load nco 70 module load nco/4.9.1 71 module load cdo/1.9.5 72 72 module load python3 73 73 module load datadir/igcmg -
TOOLS/MOSAIX/CreateWeightsMask.bash
r6105 r6107 3 3 #MSUB -o Out_WeightsMask # Standard output 4 4 #MSUB -e Out_WeightsMask # Error output 5 #MSUB -n 4# Number of processors5 #MSUB -n 1 # Number of processors 6 6 #MSUB -T 7200 # Time limit (seconds) 7 7 #MSUB -Q normal … … 73 73 #CplModel=ORCA2.3xICO30 74 74 #CplModel=ORCA2.3xICO40 75 CplModel=eORCA1.2xLMD14414275 #CplModel=eORCA1.2xLMD144142 76 76 #CplModel=eORCA1.4.0xLMD144142 77 77 #CplModel=eORCA1.4.2xLMD144142 78 78 #CplModel=eORCA1.2xLMD256256 79 79 #CplModel=eORCA1.2xICO40 80 CplModel=eORCA1.4.2xICO40 80 81 #CplModel=eORCA1.2xICO450 81 82 #CplModel=eORCA025.1xLMD256256 … … 127 128 # Each item in CommandList describes the properties of interpolation weights to generate. 128 129 # White spaces separate analysis. 130 # 129 131 # \!/ No spaces in any analysis \!/ 130 132 # … … 178 180 #CommandList=( Calving ) 179 181 #CommandList=( Grids ) 180 181 if [[ ${Version} = test_runoff_* ]] ; then182 CommandList=( Runoff )183 fi184 182 185 183 ## =========================================================================== … … 281 279 module purge 282 280 source ${SUBMIT_DIR}/arch.env 283 module load nco/4.9.2 281 module load nco/4.9.1 282 module load cdo/1.9.5 284 283 module load python3 285 284 module load datadir/igcmg 286 285 module list 286 cp ${SUBMIT_DIR}/arch.env . ${TMPDIR} 287 287 set -e 288 288 ;; … … 291 291 TMPDIR=${HOME}/Scratch/TMP 292 292 SUBMIT_DIR=$(pwd) 293 MpiRun="/opt/local/bin/mpirun -n 4"293 MpiRun="/opt/local/bin/mpirun -n 1" 294 294 PyRun="time" 295 295 ;; … … 319 319 case ${OCE} in # Periodicity type of ORCA grid 320 320 ( ORCA2.4 ) OcePerio=4.2 ;; # 4 ORCA 6 PALEORCA 321 ( ORCA2* ) OcePerio=4 ;; 321 ( ORCA2* ) OcePerio=4 ;; # 4 ORCA 6 PALEORCA 322 322 ( eORCA1.2 ) OcePerio=6 ;; 323 323 ( eORCA1.4.0 ) OcePerio=6 ;; … … 334 334 cp ${SUBMIT_DIR}/iodef_oce_to_atm.xml . 335 335 336 cp ${R_IN}/OCE/NEMO/${OCE}/${OCE}_coordinates_mask.nc . 336 [[ -f ${R_IN}/OCE/NEMO/${OCE}/${OCE}_coordinates_mask.nc ]] && cp ${R_IN}/OCE/NEMO/${OCE}/${OCE}_coordinates_mask.nc . 337 [[ -f ${R_IN}/OCE/NEMO/${OCE}/GRIDS/${OCE}_coordinates_mask.nc ]] && cp ${R_IN}/OCE/NEMO/${OCE}/GRIDS/${OCE}_coordinates_mask.nc . 338 if [[ ! -f ${OCE}_coordinates_mask.nc ]] ; then 339 exit 1 340 fi 341 337 342 cp ${R_IN}/ATM/GRID/${ATM}_grid.nc . 338 343 -
TOOLS/MOSAIX/src/MOSAIX/interpol.f90
r5154 r6107 105 105 CALL xios_get_domain_attr ("domain_src", lonvalue_1d=lon_src(:,1), latvalue_1d=lat_src(1,:) ) 106 106 lon_src (:,:) = SPREAD ( lon_src(:,1), DIM=2, ncopies=nj_src) 107 lat_src (:,:) = SPREAD ( lat_src(1,:), DIM=1, ncopies=n j_src)107 lat_src (:,:) = SPREAD ( lat_src(1,:), DIM=1, ncopies=ni_src) 108 108 CASE default 109 109 IF ( nj_src == 1 .AND. ni_src == 1) THEN … … 125 125 ALLOCATE ( imask_src (ni_src, nj_src), lmask_src (ni_src, nj_src) ) 126 126 IF ( l_mask_src ) THEN 127 WRITE (UNIT=nout, FMT="('-- ', 1I4.4, ' -- Receive field mask_src ')" ) rank127 WRITE (UNIT=nout, FMT="('-- ', 1I4.4, ' -- Receive field mask_src : ', 2I8)" ) rank, SHAPE (imask_src) 128 128 CALL xios_recv_field ("mask_src", imask_src) 129 129 lmask_src = .FALSE. … … 136 136 !< Read area on the source grid 137 137 IF ( l_use_area ) THEN 138 IF ( l_src_is0D ) ALLOCATE ( area_src ( 1 , 1 ) ) 139 IF ( l_src_is1D ) ALLOCATE ( area_src ( 1 , ni_src) ) 140 IF ( l_src_is2D ) ALLOCATE ( area_src (ni_src, nj_src) ) 141 WRITE (UNIT=nout, FMT="('-- ', 1I4.4, ' -- Receive field area_src')" ) rank 138 ALLOCATE ( area_src (ni_src, nj_src) ) 142 139 CALL xios_recv_field ("area_src", area_src(:,:)) 143 WRITE (UNIT=nout, FMT="('-- ', 1I4.4, ' -- Sum area_src : ', 2E15.3)" ) rank, SUM(area_src) 144 ENDIF 145 140 ENDIF 141 146 142 !< Read characteristics of the destination grid 147 143 WRITE (UNIT=nout, FMT="('-- ', 1I4.4, ' -- Reading domain_dst')" ) rank … … 171 167 ALLOCATE ( imask_dst (ni_dst, nj_dst), lmask_dst (ni_dst, nj_dst) ) 172 168 IF ( l_mask_dst ) THEN 173 WRITE (UNIT=nout, FMT="('-- ', 1I4.4, ' -- Receive field mask_dst ')" ) rank169 WRITE (UNIT=nout, FMT="('-- ', 1I4.4, ' -- Receive field mask_dst : ', 2I8)" ) rank, SHAPE (imask_dst) 174 170 CALL xios_recv_field ("mask_dst", imask_dst) 175 171 lmask_dst = .FALSE. … … 182 178 !< Read area on the destination grid 183 179 IF ( l_use_area ) THEN 184 IF ( l_dst_is0D ) ALLOCATE ( area_dst ( 1 , 1 ) ) 185 IF ( l_dst_is1D ) ALLOCATE ( area_dst ( 1 , ni_dst) ) 186 IF ( l_dst_is2D ) ALLOCATE ( area_dst (ni_dst, nj_dst) ) 187 WRITE (UNIT=nout, FMT="('-- ', 1I4.4, ' -- receive field area_dst')" ) rank 180 ALLOCATE ( area_dst (ni_dst, nj_dst) ) 188 181 CALL xios_recv_field ("area_dst", area_dst(:,:)) 189 WRITE (UNIT=nout, FMT="('-- ', 1I4.4, ' -- Sum area_ dst : ', 2E15.3)" ) rank, SUM(area_dst)182 WRITE (UNIT=nout, FMT="('-- ', 1I4.4, ' -- Sum area_src : ', 2E15.3)" ) rank, SUM(area_src) 190 183 ENDIF 191 184 … … 213 206 CALL xios_close_context_definition () 214 207 208 WRITE (UNIT=nout, FMT="('-- ', 1I4.4, ' -- Update calendar')" ) rank 215 209 CALL xios_update_calendar (1) 216 210 WRITE (UNIT=nout, FMT="('-- ', 1I4.4, ' -- Send field mask_src')" ) rank 217 211 CALL xios_send_field ("mask_src", imask_src) 218 212 219 !!< Creates analytic fields and interpolate213 WRITE (UNIT=nout, FMT="('-- ', 1I4.4, ' -- Creates analytic fields and interpolate')" ) rank 220 214 ALLOCATE ( field_src (ni_src, nj_src) ) 221 215
Note: See TracChangeset
for help on using the changeset viewer.