New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 7521 for branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC – NEMO

Ignore:
Timestamp:
2016-12-29T16:20:07+01:00 (7 years ago)
Author:
cbricaud
Message:

correc bugs in crs branch

File:
1 edited

Legend:

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

    r7520 r7521  
    139139         CASE ( 'T' ) 
    140140            DO jj =  nldj_crs, nlej_crs 
    141                ijj = mjs_crs(jj) + 1 
     141               ijj = mjs_crs(jj) + + INT(0.5*nfacty(jj)) 
    142142               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)) 
    145144                  p_gphi_crs(ji,jj) = p_gphi(iji,ijj) 
    146145                  p_glam_crs(ji,jj) = p_glam(iji,ijj) 
     
    149148         CASE ( 'U' ) 
    150149            DO jj =  nldj_crs, nlej_crs 
    151                ijj = mjs_crs(jj) + 1 
     150               ijj = mjs_crs(jj) + INT(0.5*nfacty(jj)) 
    152151               DO ji = nldi_crs, nlei_crs 
    153152                  iji = mie_crs(ji) 
     
    161160               ijj = mje_crs(jj) 
    162161               DO ji = nldi_crs, nlei_crs 
    163                   iji = mis_crs(ji) + 1 
     162                  iji = mis_crs(ji) + INT(0.5*nfactx(ji)) 
    164163                  p_gphi_crs(ji,jj) = p_gphi(iji,ijj) 
    165164                  p_glam_crs(ji,jj) = p_glam(iji,ijj) 
     
    212211      INTEGER :: ji, jj, jk     ! dummy loop indices 
    213212      INTEGER :: ijis,ijie,ijjs,ijje 
    214       INTEGER :: ji1, jj1 
     213      INTEGER :: i1, j1 
    215214   
    216215      !!----------------------------------------------------------------   
     
    227226               ijje = mje_crs(jj) 
    228227 
     228               i1=INT(0.5*nfactx(ji)) 
     229               j1=INT(0.5*nfacty(jj)) 
     230 
    229231               ! Only for a factro 3 coarsening 
    230232               SELECT CASE ( cd_type ) 
    231233                   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) 
    234236                   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)  
    237239 
    238240                   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)  
    241243                   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)  
    244246               END SELECT 
    245247            ENDDO 
     
    11691171         ! check 
    11701172         !========================================================================== 
    1171          !WRITE(narea+8000-1,*)"nfipproc(ji,jj),narea :",nfipproc(iproci,iprocj),narea 
    1172          !WRITE(narea+8000-1,*)"proc i,j ",iproci,iprocj 
    1173          !WRITE(narea+8000-1,*)"jpni  jpnj jpnij ",jpni,jpnj,jpnij 
    1174          !WRITE(narea+8000-1,*)"nperio jperio ",nperio,jperio 
    1175          !WRITE(narea+8000-1,*)"nowe noea",nowe,noea 
    1176          !WRITE(narea+8000-1,*)"noso nono",noso,nono 
    1177          !WRITE(narea+8000-1,*)"nbondi nbondj ",nbondi,nbondj 
    1178          !WRITE(narea+8000-1,*)"jpiglo jpjglo ",jpiglo,jpjglo 
    1179          !WRITE(narea+8000-1,*)"jpi jpj ",jpi,jpj 
    1180          !WRITE(narea+8000-1,*)"nbondi nbondj",nbondi,nbondj 
    1181          !WRITE(narea+8000-1,*)"nimpp njmpp ",nimpp,njmpp 
    1182          !WRITE(narea+8000-1,*)"loc jpi nldi,nlei,nlci ",jpi, nldi        ,nlei         ,nlci 
    1183          !WRITE(narea+8000-1,*)"glo jpi nldi,nlei      ",jpi, nldi+nimpp-1,nlei+nimpp-1 
    1184          !WRITE(narea+8000-1,*)"loc jpj nldj,nlej,nlcj ",jpj, nldj        ,nlej         ,nlcj 
    1185          !WRITE(narea+8000-1,*)"glo jpj nldj,nlej      ",jpj, nldj+njmpp-1,nlej+njmpp-1 
    1186          !WRITE(narea+8000-1,*)"jpiglo_crs jpjglo_crs ",jpiglo_crs,jpjglo_crs 
    1187          !WRITE(narea+8000-1,*)"jpi_crs jpj_crs ",jpi_crs,jpj_crs 
    1188          !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) 
    11901192 
    11911193         !========================================================================== 
     
    14761478            ENDIF 
    14771479            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) 
    14781481         ENDDO 
    14791482 
     
    15181521         ! check 
    15191522         !========================================================================== 
    1520          !WRITE(narea+8000-1,*)"loc crs jpi nldi,nlei,nlci ",jpi_crs, nldi_crs            ,nlei_crs             ,nlci_crs 
    1521          !WRITE(narea+8000-1,*)"glo crs jpi nldi,nlei      ",jpi_crs, nldi_crs+nimpp_crs-1,nlei_crs+nimpp_crs-1 
    1522          !WRITE(narea+8000-1,*)"loc crs jpj nldj,nlej,nlcj ",jpj_crs, nldj_crs            ,nlej_crs             ,nlcj_crs 
    1523          !WRITE(narea+8000-1,*)"glo crs jpj nldj,nlej      ",jpj_crs, nldj_crs+njmpp_crs-1,nlej_crs+njmpp_crs-1 
     1523         !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 
    15241527 
    15251528         !========================================================================== 
Note: See TracChangeset for help on using the changeset viewer.