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 12918 for branches/UKMO/AMM15_v3_6_STABLE_package_collate_PS44/NEMOGCM/NEMO/OPA_SRC/DIA/dia25h.F90 – NEMO

Ignore:
Timestamp:
2020-05-13T14:39:24+02:00 (4 years ago)
Author:
petesykes
Message:

Land masking for operational diagnostics

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/UKMO/AMM15_v3_6_STABLE_package_collate_PS44/NEMOGCM/NEMO/OPA_SRC/DIA/dia25h.F90

    r10390 r12918  
    267267      INTEGER                          ::   i_steps                               ! no of timesteps per hour 
    268268      REAL(wp), DIMENSION(jpi,jpj    ) ::   zw2d, un_dm, vn_dm                    ! temporary workspace 
     269      REAL(wp), DIMENSION(jpi,jpj    ) ::   zw2d_inst, zw2d_temp 
    269270      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zw3d                                  ! temporary workspace 
    270271      REAL(wp), DIMENSION(jpi,jpj,3)   ::   zwtmb                                 ! temporary workspace 
    271272      INTEGER                          ::   iyear0, nimonth0,iday0                ! start year,imonth,day 
     273      INTEGER :: jkbot, jj, ji 
    272274 
    273275      !!---------------------------------------------------------------------- 
    … …  
    388390            zw3d(:,:,:) = rinsitu_t_25h(:,:,:)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 
    389391            CALL iom_put("tempis25h", zw3d)   ! in-situ temperature 
    390             zw2d(:,:) = insitu_bot_25h(:,:)*tmask(:,:,1) + zmdi*(1.0-tmask(:,:,1)) 
    391             CALL iom_put("tempisbot25h", zw2d) ! bottom in-situ temperature 
    392             zw2d(:,:) = temp_bot_25h(:,:)*tmask(:,:,1) + zmdi*(1.0-tmask(:,:,1)) 
    393             CALL iom_put("temperbot25h",zw2d) ! bottom potential temperature 
     392            DO jj = 1, jpj 
     393               DO ji = 1, jpi 
     394                  jkbot = mbkt(ji,jj) 
     395                  zw2d_inst(ji,jj) = insitu_bot_25h(ji,jj)*tmask(ji,jj,jkbot) + zmdi*(1.0-tmask(ji,jj,jkbot)) 
     396                  zw2d_temp(ji,jj) = temp_bot_25h(ji,jj)*tmask(ji,jj,jkbot) + zmdi*(1.0-tmask(ji,jj,jkbot)) 
     397               ENDDO 
     398            ENDDO 
     399            CALL iom_put("tempisbot25h", zw2d_inst) ! bottom in-situ temperature 
     400            CALL iom_put("temperbot25h",zw2d_temp) ! bottom potential temperature 
    394401            zw3d(:,:,:) = sn_25h(:,:,:)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 
    395402            CALL iom_put( "salin25h", zw3d  )   ! salinity 
Note: See TracChangeset for help on using the changeset viewer.