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 15144 for NEMO/branches/UKMO/NEMO_4.0.4_generic_obs/src/OCE/OBS/obs_oper.F90 – NEMO

Ignore:
Timestamp:
2021-07-26T16:46:48+02:00 (3 years ago)
Author:
dford
Message:

Initial implementation of generic OBS interface.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0.4_generic_obs/src/OCE/OBS/obs_oper.F90

    r15089 r15144  
    178178         ! Initialize daily mean for first timestep of the day 
    179179         IF ( idayend == 1 .OR. kt == 0 ) THEN 
    180             DO_3D( 1, 1, 1, 1, 1, jpk ) 
    181                prodatqc%vdmean(ji,jj,jk,kvar) = 0.0 
    182             END_3D 
    183          ENDIF 
    184  
    185          DO_3D( 1, 1, 1, 1, 1, jpk ) 
    186             ! Increment field 1 for computing daily mean 
    187             prodatqc%vdmean(ji,jj,jk,kvar) = prodatqc%vdmean(ji,jj,jk,kvar) & 
    188                &                           + pvar(ji,jj,jk) 
    189          END_3D 
     180            DO jk = 1, jpk 
     181               DO jj = 1, jpj 
     182                  DO ji = 1, jpi 
     183                     prodatqc%vdmean(ji,jj,jk,kvar) = 0.0 
     184                  END DO 
     185               END DO 
     186            END DO 
     187         ENDIF 
     188 
     189         DO jk = 1, jpk 
     190            DO jj = 1, jpj 
     191               DO ji = 1, jpi 
     192                  ! Increment field 1 for computing daily mean 
     193                  prodatqc%vdmean(ji,jj,jk,kvar) = prodatqc%vdmean(ji,jj,jk,kvar) & 
     194                     &                           + pvar(ji,jj,jk) 
     195               END DO 
     196            END DO 
     197         END DO 
    190198 
    191199         ! Compute the daily mean at the end of day 
     
    194202            IF (lwp) WRITE(numout,*) 'Calculating prodatqc%vdmean on time-step: ',kt 
    195203            CALL FLUSH(numout) 
    196             DO_3D( 1, 1, 1, 1, 1, jpk ) 
    197                prodatqc%vdmean(ji,jj,jk,kvar) = prodatqc%vdmean(ji,jj,jk,kvar) & 
    198                   &                           * zdaystp 
    199             END_3D 
     204            DO jk = 1, jpk 
     205               DO jj = 1, jpj 
     206                  DO ji = 1, jpi 
     207                     prodatqc%vdmean(ji,jj,jk,kvar) = prodatqc%vdmean(ji,jj,jk,kvar) & 
     208                        &                           * zdaystp 
     209                  END DO 
     210               END DO 
     211            END DO 
    200212         ENDIF 
    201213 
Note: See TracChangeset for help on using the changeset viewer.