- Timestamp:
- 2016-12-29T16:20:07+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/CRS/crsdom.F90
r7520 r7521 139 139 CASE ( 'T' ) 140 140 DO jj = nldj_crs, nlej_crs 141 ijj = mjs_crs(jj) + 1141 ijj = mjs_crs(jj) + + INT(0.5*nfacty(jj)) 142 142 DO ji = nldi_crs, nlei_crs 143 iji = mis_crs(ji) + 1 144 IF( ijj .GT. jpj )WRITE(narea+8000-1,*)"BUG ijj ",jj,mjs_crs(jj);CALL FLUSH(narea+8000-1) 143 iji = mis_crs(ji) + INT(0.5*nfactx(ji)) 145 144 p_gphi_crs(ji,jj) = p_gphi(iji,ijj) 146 145 p_glam_crs(ji,jj) = p_glam(iji,ijj) … … 149 148 CASE ( 'U' ) 150 149 DO jj = nldj_crs, nlej_crs 151 ijj = mjs_crs(jj) + 1150 ijj = mjs_crs(jj) + INT(0.5*nfacty(jj)) 152 151 DO ji = nldi_crs, nlei_crs 153 152 iji = mie_crs(ji) … … 161 160 ijj = mje_crs(jj) 162 161 DO ji = nldi_crs, nlei_crs 163 iji = mis_crs(ji) + 1162 iji = mis_crs(ji) + INT(0.5*nfactx(ji)) 164 163 p_gphi_crs(ji,jj) = p_gphi(iji,ijj) 165 164 p_glam_crs(ji,jj) = p_glam(iji,ijj) … … 212 211 INTEGER :: ji, jj, jk ! dummy loop indices 213 212 INTEGER :: ijis,ijie,ijjs,ijje 214 INTEGER :: ji1, jj1213 INTEGER :: i1, j1 215 214 216 215 !!---------------------------------------------------------------- … … 227 226 ijje = mje_crs(jj) 228 227 228 i1=INT(0.5*nfactx(ji)) 229 j1=INT(0.5*nfacty(jj)) 230 229 231 ! Only for a factro 3 coarsening 230 232 SELECT CASE ( cd_type ) 231 233 CASE ( 'T' ) 232 p_e1_crs(ji,jj) = REAL(nn_factx,wp)*p_e1(ijis+ 1,ijjs+1)233 p_e2_crs(ji,jj) = REAL(nn_facty,wp)*p_e2(ijis+ 1,ijjs+1)234 p_e1_crs(ji,jj) = REAL(nn_factx,wp)*p_e1(ijis+i1,ijjs+j1) 235 p_e2_crs(ji,jj) = REAL(nn_facty,wp)*p_e2(ijis+i1,ijjs+j1) 234 236 CASE ( 'U' ) 235 p_e1_crs(ji,jj) = REAL(nn_factx,wp)*p_e1(ijis+ 1,ijjs+1)236 p_e2_crs(ji,jj) = REAL(nn_facty,wp)*p_e2(ijie ,ijjs+1)237 p_e1_crs(ji,jj) = REAL(nn_factx,wp)*p_e1(ijis+i1,ijjs+j1) 238 p_e2_crs(ji,jj) = REAL(nn_facty,wp)*p_e2(ijie ,ijjs+j1) 237 239 238 240 CASE ( 'V' ) 239 p_e1_crs(ji,jj) = REAL(nn_factx,wp)*p_e1(ijis+ 1,ijje)240 p_e2_crs(ji,jj) = REAL(nn_facty,wp)*p_e2(ijis+ 1,ijjs+1)241 p_e1_crs(ji,jj) = REAL(nn_factx,wp)*p_e1(ijis+i1,ijje ) 242 p_e2_crs(ji,jj) = REAL(nn_facty,wp)*p_e2(ijis+i1,ijjs+j1) 241 243 CASE ( 'F' ) 242 p_e1_crs(ji,jj) = REAL(nn_factx,wp)*p_e1(ijis+ 1,ijje)243 p_e2_crs(ji,jj) = REAL(nn_facty,wp)*p_e2(ijie ,ijjs+1)244 p_e1_crs(ji,jj) = REAL(nn_factx,wp)*p_e1(ijis+i1,ijje ) 245 p_e2_crs(ji,jj) = REAL(nn_facty,wp)*p_e2(ijie ,ijjs+j1) 244 246 END SELECT 245 247 ENDDO … … 1169 1171 ! check 1170 1172 !========================================================================== 1171 !WRITE(narea+ 8000-1,*)"nfipproc(ji,jj),narea :",nfipproc(iproci,iprocj),narea1172 !WRITE(narea+ 8000-1,*)"proc i,j ",iproci,iprocj1173 !WRITE(narea+ 8000-1,*)"jpni jpnj jpnij ",jpni,jpnj,jpnij1174 !WRITE(narea+ 8000-1,*)"nperio jperio ",nperio,jperio1175 !WRITE(narea+ 8000-1,*)"nowe noea",nowe,noea1176 !WRITE(narea+ 8000-1,*)"noso nono",noso,nono1177 !WRITE(narea+ 8000-1,*)"nbondi nbondj ",nbondi,nbondj1178 !WRITE(narea+ 8000-1,*)"jpiglo jpjglo ",jpiglo,jpjglo1179 !WRITE(narea+ 8000-1,*)"jpi jpj ",jpi,jpj1180 !WRITE(narea+ 8000-1,*)"nbondi nbondj",nbondi,nbondj1181 !WRITE(narea+ 8000-1,*)"nimpp njmpp ",nimpp,njmpp1182 !WRITE(narea+ 8000-1,*)"loc jpi nldi,nlei,nlci ",jpi, nldi ,nlei ,nlci1183 !WRITE(narea+ 8000-1,*)"glo jpi nldi,nlei ",jpi, nldi+nimpp-1,nlei+nimpp-11184 !WRITE(narea+ 8000-1,*)"loc jpj nldj,nlej,nlcj ",jpj, nldj ,nlej ,nlcj1185 !WRITE(narea+ 8000-1,*)"glo jpj nldj,nlej ",jpj, nldj+njmpp-1,nlej+njmpp-11186 !WRITE(narea+ 8000-1,*)"jpiglo_crs jpjglo_crs ",jpiglo_crs,jpjglo_crs1187 !WRITE(narea+ 8000-1,*)"jpi_crs jpj_crs ",jpi_crs,jpj_crs1188 !WRITE(narea+ 8000-1,*)"glamt gphit ",glamt(1,1),gphit(jpi,jpj),glamt(jpi,jpj),gphit(jpi,jpj)1189 !WRITE(narea+ 8000-1,*)"min max tmask ",MINVAL(tmask),MAXVAL(tmask)1173 !WRITE(narea+1000-1,*)"nfipproc(ji,jj),narea :",nfipproc(iproci,iprocj),narea 1174 !WRITE(narea+1000-1,*)"proc i,j ",iproci,iprocj 1175 !WRITE(narea+1000-1,*)"jpni jpnj jpnij ",jpni,jpnj,jpnij 1176 !WRITE(narea+1000-1,*)"nperio jperio ",nperio,jperio 1177 !WRITE(narea+1000-1,*)"nowe noea",nowe,noea 1178 !WRITE(narea+1000-1,*)"noso nono",noso,nono 1179 !WRITE(narea+1000-1,*)"nbondi nbondj ",nbondi,nbondj 1180 !WRITE(narea+1000-1,*)"jpiglo jpjglo ",jpiglo,jpjglo 1181 !WRITE(narea+1000-1,*)"jpi jpj ",jpi,jpj 1182 !WRITE(narea+1000-1,*)"nbondi nbondj",nbondi,nbondj 1183 !WRITE(narea+1000-1,*)"nimpp njmpp ",nimpp,njmpp 1184 !WRITE(narea+1000-1,*)"loc jpi nldi,nlei,nlci ",jpi, nldi ,nlei ,nlci 1185 !WRITE(narea+1000-1,*)"glo jpi nldi,nlei ",jpi, nldi+nimpp-1,nlei+nimpp-1 1186 !WRITE(narea+1000-1,*)"loc jpj nldj,nlej,nlcj ",jpj, nldj ,nlej ,nlcj 1187 !WRITE(narea+1000-1,*)"glo jpj nldj,nlej ",jpj, nldj+njmpp-1,nlej+njmpp-1 1188 !WRITE(narea+1000-1,*)"jpiglo_crs jpjglo_crs ",jpiglo_crs,jpjglo_crs 1189 !WRITE(narea+1000-1,*)"jpi_crs jpj_crs ",jpi_crs,jpj_crs 1190 !WRITE(narea+1000-1,*)"glamt gphit ",glamt(1,1),gphit(jpi,jpj),glamt(jpi,jpj),gphit(jpi,jpj) 1191 !WRITE(narea+1000-1,*)"min max tmask ",MINVAL(tmask),MAXVAL(tmask) 1190 1192 1191 1193 !========================================================================== … … 1476 1478 ENDIF 1477 1479 nfacty(jj) = mje_crs(jj)-mjs_crs(jj)+1 1480 !WRITE(narea+1000-1,*)"test J",jj,mjg_crs(jj),mjs_crs(jj),mje_crs(jj),mjs_crs(jj)+njmpp-1,mje_crs(jj)+njmpp-1,nfacty(jj) 1478 1481 ENDDO 1479 1482 … … 1518 1521 ! check 1519 1522 !========================================================================== 1520 !WRITE(narea+ 8000-1,*)"loc crs jpi nldi,nlei,nlci ",jpi_crs, nldi_crs ,nlei_crs ,nlci_crs1521 !WRITE(narea+ 8000-1,*)"glo crs jpi nldi,nlei ",jpi_crs, nldi_crs+nimpp_crs-1,nlei_crs+nimpp_crs-11522 !WRITE(narea+ 8000-1,*)"loc crs jpj nldj,nlej,nlcj ",jpj_crs, nldj_crs ,nlej_crs ,nlcj_crs1523 !WRITE(narea+ 8000-1,*)"glo crs jpj nldj,nlej ",jpj_crs, nldj_crs+njmpp_crs-1,nlej_crs+njmpp_crs-11523 !WRITE(narea+1000-1,*)"loc crs jpi nldi,nlei,nlci ",jpi_crs, nldi_crs ,nlei_crs ,nlci_crs 1524 !WRITE(narea+1000-1,*)"glo crs jpi nldi,nlei ",jpi_crs, nldi_crs+nimpp_crs-1,nlei_crs+nimpp_crs-1 1525 !WRITE(narea+1000-1,*)"loc crs jpj nldj,nlej,nlcj ",jpj_crs, nldj_crs ,nlej_crs ,nlcj_crs 1526 !WRITE(narea+1000-1,*)"glo crs jpj nldj,nlej ",jpj_crs, nldj_crs+njmpp_crs-1,nlej_crs+njmpp_crs-1 1524 1527 1525 1528 !==========================================================================
Note: See TracChangeset
for help on using the changeset viewer.