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 13540 for NEMO/branches/2020/r12377_ticket2386/src/SAO/sao_read.F90 – NEMO

Ignore:
Timestamp:
2020-09-29T12:41:06+02:00 (4 years ago)
Author:
andmirek
Message:

Ticket #2386: update to latest trunk

Location:
NEMO/branches/2020/r12377_ticket2386
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/r12377_ticket2386

    • Property svn:externals
      •  

        old new  
        33^/utils/build/mk@HEAD         mk 
        44^/utils/tools@HEAD            tools 
        5 ^/vendors/AGRIF/dev@HEAD      ext/AGRIF 
         5^/vendors/AGRIF/dev_r12970_AGRIF_CMEMS      ext/AGRIF 
        66^/vendors/FCM@HEAD            ext/FCM 
        77^/vendors/IOIPSL@HEAD         ext/IOIPSL 
        88 
        99# SETTE 
        10 ^/utils/CI/sette@HEAD         sette 
         10^/utils/CI/sette@13507        sette 
  • NEMO/branches/2020/r12377_ticket2386/src/SAO/sao_read.F90

    r10069 r13540  
    1010   USE netcdf 
    1111   USE oce,     ONLY: tsn, sshn 
    12    USE dom_oce, ONLY: nlci, nlcj, nimpp, njmpp, tmask 
     12   USE dom_oce, ONLY: nimpp, njmpp, tmask 
    1313   USE par_oce, ONLY: jpi, jpj, jpk 
    1414   ! 
     
    9494         IF (ifcst .LE. ntimes) THEN 
    9595            ! Allocate temporary temperature array 
    96             ALLOCATE(temp_tn(nlci,nlcj,jpk)) 
    97             ALLOCATE(temp_sn(nlci,nlcj,jpk)) 
    98             ALLOCATE(temp_sshn(nlci,nlcj)) 
     96            ALLOCATE(temp_tn(jpi,jpj,jpk)) 
     97            ALLOCATE(temp_sn(jpi,jpj,jpk)) 
     98            ALLOCATE(temp_sshn(jpi,jpj)) 
    9999 
    100100            ! Set temp_tn, temp_sn to 0. 
     
    104104 
    105105            ! Create start and count arrays 
    106             start_n = (/ nimpp, njmpp, 1,  ifcst /) 
    107             count_n = (/ nlci,  nlcj,  jpk, 1     /) 
    108             start_s = (/ nimpp, njmpp, ifcst /) 
    109             count_s = (/ nlci,  nlcj,  1     /) 
     106            start_n = (/ nimpp, njmpp,      1, ifcst /) 
     107            count_n = (/   jpi,   jpj, jpk, 1        /) 
     108            start_s = (/ nimpp, njmpp        , ifcst /) 
     109            count_s = (/   jpi,   jpj,      1        /) 
    110110 
    111111            ! Read information into temporary arrays 
     
    138138 
    139139            ! Mask out missing data index 
    140             tsn(1:nlci,1:nlcj,1:jpk,1) = temp_tn(:,:,:) * tmask(1:nlci,1:nlcj,1:jpk) 
    141             tsn(1:nlci,1:nlcj,1:jpk,2) = temp_sn(:,:,:) * tmask(1:nlci,1:nlcj,1:jpk) 
    142             sshn(1:nlci,1:nlcj)        = temp_sshn(:,:) * tmask(1:nlci,1:nlcj,1) 
    143  
    144             ! Remove halo from tmask, tsn, sshn to prevent double obs counting 
    145             IF (jpi > nlci) THEN 
    146                 tmask(nlci+1:,:,:) = 0 
    147                 tsn(nlci+1:,:,:,1) = 0 
    148                 tsn(nlci+1:,:,:,2) = 0 
    149                 sshn(nlci+1:,:) = 0 
    150             END IF 
    151             IF (jpj > nlcj) THEN 
    152                 tmask(:,nlcj+1:,:) = 0 
    153                 tsn(:,nlcj+1:,:,1) = 0 
    154                 tsn(:,nlcj+1:,:,2) = 0 
    155                 sshn(:,nlcj+1:) = 0 
    156             END IF 
    157  
     140            tsn(1:jpi,1:jpj,1:jpk,1) = temp_tn(:,:,:) * tmask(1:jpi,1:jpj,1:jpk) 
     141            tsn(1:jpi,1:jpj,1:jpk,2) = temp_sn(:,:,:) * tmask(1:jpi,1:jpj,1:jpk) 
     142            sshn(1:jpi,1:jpj)        = temp_sshn(:,:) * tmask(1:jpi,1:jpj,1) 
     143             
    158144            ! Deallocate arrays 
    159145            DEALLOCATE(temp_tn, temp_sn, temp_sshn) 
Note: See TracChangeset for help on using the changeset viewer.