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 3211 for branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/OBS/obs_grid.F90 – NEMO

Ignore:
Timestamp:
2011-12-11T16:00:26+01:00 (13 years ago)
Author:
spickles2
Message:

Stephen Pickles, 11 Dec 2011

Commit to bring the rest of the DCSE NEMO development branch
in line with the latest development version. This includes
array index re-ordering of all OPA_SRC/.

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/OBS/obs_grid.F90

    r2715 r3211  
    8585      & grid_search_file    ! file name head for grid search lookup  
    8686 
     87   !! * Control permutation of array indices 
     88#  include "dom_oce_ftrans.h90" 
     89 
    8790   !!---------------------------------------------------------------------- 
    8891   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    … …  
    128131         ELSE 
    129132            IF ( cdgrid == 'T' ) THEN 
    130                CALL obs_grd_bruteforce( jpi, jpj, jpiglo, jpjglo, & 
    131                   &                             nldi, nlei,nldj,  nlej,   & 
    132                   &                             nproc, jpnij,             & 
    133                   &                             glamt, gphit, tmask,      & 
    134                   &                             kobsin, plam, pphi,       & 
    135                   &                             kobsi, kobsj, kproc ) 
     133               CALL obs_grd_bruteforce( jpi, jpj, jpiglo, jpjglo,        & 
     134                  &                     nldi, nlei,nldj,  nlej,          & 
     135                  &                     nproc, jpnij,                    & 
     136#if defined key_z_first 
     137                  &                     glamt, gphit, tmask_1(:,:),      & 
     138#else 
     139                  &                     glamt, gphit, tmask(:,:,1),      & 
     140#endif 
     141                  &                     kobsin, plam, pphi,              & 
     142                  &                     kobsi, kobsj, kproc ) 
    136143            ELSEIF ( cdgrid == 'U' ) THEN 
    137                CALL obs_grd_bruteforce( jpi, jpj, jpiglo, jpjglo, & 
    138                   &                             nldi, nlei,nldj,  nlej,   & 
    139                   &                             nproc, jpnij,             & 
    140                   &                             glamu, gphiu, umask,      & 
    141                   &                             kobsin, plam, pphi,       & 
    142                   &                             kobsi, kobsj, kproc ) 
     144               CALL obs_grd_bruteforce( jpi, jpj, jpiglo, jpjglo,        & 
     145                  &                     nldi, nlei,nldj,  nlej,          & 
     146                  &                     nproc, jpnij,                    & 
     147#if defined key_z_first 
     148                  &                     glamu, gphiu, umask_1(:,:),      & 
     149#else 
     150                  &                     glamu, gphiu, umask(:,:,1),      & 
     151#endif 
     152                  &                     kobsin, plam, pphi,              & 
     153                  &                     kobsi, kobsj, kproc ) 
    143154            ELSEIF ( cdgrid == 'V' ) THEN 
    144                CALL obs_grd_bruteforce( jpi, jpj, jpiglo, jpjglo, & 
    145                   &                             nldi, nlei,nldj,  nlej,   & 
    146                   &                             nproc, jpnij,             & 
    147                   &                             glamv, gphiv, vmask,      & 
    148                   &                             kobsin, plam, pphi,       & 
    149                   &                             kobsi, kobsj, kproc ) 
     155               CALL obs_grd_bruteforce( jpi, jpj, jpiglo, jpjglo,        & 
     156                  &                     nldi, nlei,nldj,  nlej,          & 
     157                  &                     nproc, jpnij,                    & 
     158#if defined key_z_first 
     159                  &                     glamv, gphiv, vmask_1(:,:),      & 
     160#else 
     161                  &                     glamv, gphiv, vmask(:,:,1),      & 
     162#endif 
     163                  &                     kobsin, plam, pphi,              & 
     164                  &                     kobsi, kobsj, kproc ) 
    150165            ELSEIF ( cdgrid == 'F' ) THEN 
    151                CALL obs_grd_bruteforce( jpi, jpj, jpiglo, jpjglo, & 
    152                   &                             nldi, nlei,nldj,  nlej,   & 
    153                   &                             nproc, jpnij,             & 
    154                   &                             glamf, gphif, fmask,      & 
    155                   &                             kobsin, plam, pphi,       & 
    156                   &                             kobsi, kobsj, kproc ) 
     166               CALL obs_grd_bruteforce( jpi, jpj, jpiglo, jpjglo,        & 
     167                  &                     nldi, nlei,nldj,  nlej,          & 
     168                  &                     nproc, jpnij,                    & 
     169#if defined key_z_first 
     170                  &                     glamf, gphif, fmask_1(:,:),      & 
     171#else 
     172                  &                     glamf, gphif, fmask(:,:,1),      & 
     173#endif 
     174                  &                     kobsin, plam, pphi,              & 
     175                  &                     kobsi, kobsj, kproc ) 
    157176            ELSE 
    158177               CALL ctl_stop( 'Grid not supported' ) 
    … …  
    283302               zlamg(mig(ji),mjg(jj)) = glamt(ji,jj) 
    284303               zphig(mig(ji),mjg(jj)) = gphit(ji,jj) 
     304#if defined key_z_first 
     305               zmskg(mig(ji),mjg(jj)) = tmask_1(ji,jj) 
     306#else 
    285307               zmskg(mig(ji),mjg(jj)) = tmask(ji,jj,1) 
     308#endif 
    286309            END DO 
    287310         END DO 
    … …  
    295318               zlamg(ji,jj) = glamt(ji,jj) 
    296319               zphig(ji,jj) = gphit(ji,jj) 
     320#if defined key_z_first 
     321               zmskg(ji,jj) = tmask_1(ji,jj) 
     322#else 
    297323               zmskg(ji,jj) = tmask(ji,jj,1) 
     324#endif 
    298325            END DO 
    299326         END DO 
    … …  
    813840            END DO 
    814841             
    815             CALL obs_grd_bruteforce( jpi, jpj, jpiglo, jpjglo,  & 
    816                &                     nldi, nlei,nldj,  nlej,    & 
    817                &                     nproc, jpnij,              & 
    818                &                     glamt, gphit, tmask,      & 
    819                &                     nlons*nlats, lonsi, latsi, & 
     842            CALL obs_grd_bruteforce( jpi, jpj, jpiglo, jpjglo,   & 
     843               &                     nldi, nlei,nldj,  nlej,     & 
     844               &                     nproc, jpnij,               & 
     845               &                     glamt, gphit, tmask(:,:,1), & 
     846               &                     nlons*nlats, lonsi, latsi,  & 
    820847               &                     ixposi, iyposi, iproci ) 
    821848             
Note: See TracChangeset for help on using the changeset viewer.