Changeset 6809


Ignore:
Timestamp:
2016-07-19T16:30:01+02:00 (5 years ago)
Author:
kingr
Message:

Removed gcx from assim background write. Variable=0, and is not read by assim. Added if def key around allocation of en for assim bkg write.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_25hr_mean_assim_bkg/NEMOGCM/NEMO/OPA_SRC/ASM/asmbkg.F90

    r6802 r6809  
    6363   REAL(wp),SAVE, ALLOCATABLE,   DIMENSION(:,:,:) ::   en_tavg 
    6464#endif 
    65    REAL(wp),SAVE, ALLOCATABLE,   DIMENSION(:,:)   ::   sshn_tavg, gcx_tavg 
     65   REAL(wp),SAVE, ALLOCATABLE,   DIMENSION(:,:)   ::   sshn_tavg 
    6666   REAL(wp),SAVE :: numtimes_tavg     ! No of times to average over 
    6767 
     
    123123         ENDIF 
    124124         sshn_tavg=0 
     125#if defined key_zdftke 
    125126         ALLOCATE( en_tavg(jpi,jpj,jpk), STAT=ierror ) 
    126127         IF( ierror > 0 ) THEN 
     
    128129         ENDIF 
    129130         en_tavg=0 
     131#endif 
    130132         ALLOCATE( avt_tavg(jpi,jpj,jpk), STAT=ierror ) 
    131133         IF( ierror > 0 ) THEN 
    132134            CALL ctl_stop( 'asm_wri_bkg: unable to allocate avt_tavg' )   ;   RETURN 
    133135         ENDIF 
    134          ALLOCATE( gcx_tavg(jpi,jpj), STAT=ierror ) 
    135          IF( ierror > 0 ) THEN 
    136             CALL ctl_stop( 'asm_wri_bkg: unable to allocate gcx_tavg' )   ;   RETURN 
    137          ENDIF 
    138          gcx_tavg=0     
    139136         avt_tavg=0 
    140137          
     
    154151         un_tavg(:,:,:)        = un_tavg(:,:,:) + un(:,:,:) / numtimes_tavg 
    155152         vn_tavg(:,:,:)        = vn_tavg(:,:,:) + vn(:,:,:) / numtimes_tavg 
    156 !         gcx_tavg(:,:)         = gcx_tavg(:,:)  + gcx(:,:) / numtimes_tavg 
    157153         avt_tavg(:,:,:)        = avt_tavg(:,:,:) + avt(:,:,:) / numtimes_tavg 
    158154#if defined key_zdftke 
     
    199195               CALL iom_rstput( kt, nitavgbkg_r, inum, 'en'     , en_tavg ) 
    200196#endif 
    201                CALL iom_rstput( kt, nitavgbkg_r, inum, 'gcx'    , gcx_tavg) 
    202197               CALL iom_rstput( kt, nitavgbkg_r, inum, 'avt'    , avt_tavg) 
    203198               ! 
     
    223218               CALL iom_rstput( kt, nitbkg_r, inum, 'en'     , en                ) 
    224219#endif 
    225                CALL iom_rstput( kt, nitbkg_r, inum, 'gcx'    , gcx               ) 
    226220               CALL iom_rstput( kt, nitbkg_r, inum, 'avt'    , avt               ) 
    227221               ! 
Note: See TracChangeset for help on using the changeset viewer.