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 5931 for branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcsed_medusa.F90 – NEMO

Ignore:
Timestamp:
2015-11-26T17:28:18+01:00 (8 years ago)
Author:
jpalmier
Message:

JPALM --26-11-2015 -- Update MEDUSA diagnostics with iom_use

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcsed_medusa.F90

    r5841 r5931  
    2323   !! AXY (10/02/09) 
    2424   USE iom 
     25    USE trc_nam_dia         ! JPALM 13-11-2015 -- if iom_use for diag 
     26    USE trc_nam_iom_medusa  ! JPALM 13-11-2015 -- if iom_use for diag 
    2527 
    2628   USE lbclnk 
     
    9092 
    9193      CHARACTER (len=25) :: charout 
     94       
     95      !! JPALM - 26-11-2015 -add iom_use for diagnostic 
     96       REAL(wp), POINTER, DIMENSION(:,:  ) :: zw2d 
    9297      !!--------------------------------------------------------------------- 
    93  
     98      !! 
     99      IF( lk_iomput) THEN   
     100           IF( med_diag%DSED%dgsave ) THEN 
     101               CALL wrk_alloc( jpi, jpj,      zw2d ) 
     102                zw2d(:,:)      = 0.0      !! 
     103           ENDIF 
     104      ENDIF 
     105       
    94106      !! AXY (10/02/09) 
    95107      jnt = 1 
     
    165177               trbio(ji,jj,jk,8) = ztra 
    166178# endif 
    167                IF( ln_diatrc ) & 
    168                   &  trc2d(ji,jj,8) = trc2d(ji,jj,8) + ztra * fse3t(ji,jj,jk) * 86400. 
     179               IF (lk_iomput .AND. .NOT. ln_diatrc) THEN 
     180                     IF( med_diag%DSED%dgsave ) THEN 
     181                         zw2d(ji,jj) = zw2d(ji,jj) + ztra * fse3t(ji,jj,jk) * 86400 
     182                      ENDIF    
     183               ELSE IF( ln_diatrc )  THEN 
     184                    trc2d(ji,jj,8) = trc2d(ji,jj,8) + ztra * fse3t(ji,jj,jk) * 86400 
     185               ENDIF     
     186                 
     187                  . 
    169188            END DO 
    170189         END DO 
     
    175194# endif 
    176195      IF( ln_diatrc ) CALL lbc_lnk( trc2d(:,:,8), 'T', 1. )      ! Lateral boundary conditions on trc2d 
    177 # if defined key_iomput 
    178       CALL iom_put( "DSED",trc2d(:,:,8) ) 
    179 # endif 
    180  
     196      !! 
     197      IF (lk_iomput .AND. .NOT. ln_diatrc) THEN 
     198           IF( med_diag%DSED%dgsave ) THEN 
     199                CALL iom_put( "DSED"  ,  zw2d) 
     200                CALL wrk_dealloc( jpi, jpj,    zw2d  ) 
     201            ENDIF 
     202      ELSE IF (lk_iomput .AND. ln_diatrc)  THEN     
     203          CALL iom_put( "DSED",trc2d(:,:,8) ) 
     204      ENDIF 
     205      !! 
    181206# if defined key_roam 
    182207 
Note: See TracChangeset for help on using the changeset viewer.