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 6016 for branches/UKMO/dev_r4650_general_vert_coord_obsoper/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90 – NEMO

Ignore:
Timestamp:
2015-12-07T17:42:57+01:00 (8 years ago)
Author:
kingr
Message:

#1642 bug fixes for general vertical coord obsoper

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/UKMO/dev_r4650_general_vert_coord_obsoper/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90

    r5838 r6016  
    10631063      IF ( ln_t3d .OR. ln_s3d ) THEN 
    10641064         DO jprofset = 1, nprofsets 
    1065             IF( ln_zco .OR. ln_zps ) THEN  
     1065            IF( (.NOT. lk_vvl) .AND. (ln_zco .OR. ln_zps) ) THEN  
     1066               IF(lwp) THEN 
     1067                  WRITE(numout,*) 'dia_obs : calling obs_pro_opt' 
     1068               ENDIF 
    10661069               IF ( ld_enact(jprofset) ) THEN  
    10671070                  CALL obs_pro_opt( prodatqc(jprofset),                     &  
    … …  
    10771080               ENDIF  
    10781081            ELSE 
     1082               IF(lwp) THEN 
     1083                  WRITE(numout,*) 'dia_obs : calling obs_pro_sco_opt' 
     1084               ENDIF 
    10791085               IF ( ld_enact(jprofset) ) THEN  
    10801086                  CALL obs_pro_sco_opt( prodatqc(jprofset),                 &  
    10811087                     &              kstp, jpi, jpj, jpk, nit000, idaystp,   &  
    10821088                     &              tsn(:,:,:,jp_tem), tsn(:,:,:,jp_sal),   &  
    1083                      &              fsdept(:,:,:), tmask, n1dint, n2dint,   &  
     1089                     &              fsdept(:,:,:), fsdepw(:,:,:),           &  
     1090                     &              tmask, n1dint, n2dint,                  &  
    10841091                     &              kdailyavtypes = endailyavtypes )  
    10851092               ELSE  
    … …  
    10871094                     &              kstp, jpi, jpj, jpk, nit000, idaystp,   &  
    10881095                     &              tsn(:,:,:,jp_tem), tsn(:,:,:,jp_sal),   &  
    1089                      &              fsdept(:,:,:), tmask, n1dint, n2dint )  
     1096                     &              fsdept(:,:,:), fsdepw(:,:,:),           & 
     1097                     &              tmask, n1dint, n2dint )  
    10901098               ENDIF  
    10911099            ENDIF 
Note: See TracChangeset for help on using the changeset viewer.