Ignore:
Timestamp:
2018-12-03T12:26:33+01:00 (23 months ago)
Author:
dford
Message:

Remove visibility, as now available through FABM, and switch from deprecated get_bulk to get_interior.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/AMM15_v3_6_STABLE_package_collate_bgc_diagnostics/NEMOGCM/NEMO/OPA_SRC/DIA/dia25h.F90

    r10308 r10366  
    2222   USE diatmb 
    2323#if defined key_fabm 
    24    USE trc, ONLY: trn, visib 
     24   USE trc, ONLY: trn 
    2525   USE par_fabm 
    2626   USE st2d_fabm, ONLY: fabm_st2dn 
    27    USE fabm, ONLY: fabm_get_bulk_diagnostic_data, & 
     27   USE fabm, ONLY: fabm_get_interior_diagnostic_data, & 
    2828      &            fabm_get_horizontal_diagnostic_data 
    2929#endif 
     
    5353   REAL(wp),SAVE, ALLOCATABLE,   DIMENSION(:,:,:)   :: fabm_bottom_25h 
    5454   REAL(wp),SAVE, ALLOCATABLE,   DIMENSION(:,:,:)   :: fabm_2d_25h 
    55    REAL(wp),SAVE, ALLOCATABLE,   DIMENSION(:,:,:)   :: fabm_visib_25h 
    5655#endif 
    5756   INTEGER, SAVE :: cnt_25h     ! Counter for 25 hour means 
     
    182181         CALL ctl_stop( 'dia_25h: unable to allocate fabm_2d_25h' )   ;   RETURN 
    183182      ENDIF 
    184       ALLOCATE( fabm_visib_25h(jpi,jpj,jpk), STAT=ierror ) 
    185       IF( ierror > 0 ) THEN 
    186          CALL ctl_stop( 'dia_25h: unable to allocate fabm_visib_25h' )   ;   RETURN 
    187       ENDIF 
    188183#endif  
    189184      ! ------------------------- ! 
     
    216211      END DO 
    217212      DO jn = 1, jp_fabm_3d 
    218          fabm_3d_25h(:,:,:,jn) = fabm_get_bulk_diagnostic_data(model, jn) 
     213         fabm_3d_25h(:,:,:,jn) = fabm_get_interior_diagnostic_data(model, jn) 
    219214      END DO 
    220215      DO jn = 1, jp_fabm_surface 
     
    227222         fabm_2d_25h(:,:,jn) = fabm_get_horizontal_diagnostic_data(model,jn) 
    228223      END DO 
    229       fabm_visib_25h(:,:,:) = visib(:,:,:) 
    230224#endif 
    231225#if defined key_lim3 || defined key_lim2 
     
    333327      END DO 
    334328      DO jn = 1, jp_fabm_3d 
    335          fabm_3d_25h(:,:,:,jn) = fabm_3d_25h(:,:,:,jn) + fabm_get_bulk_diagnostic_data(model, jn) 
     329         fabm_3d_25h(:,:,:,jn) = fabm_3d_25h(:,:,:,jn) + fabm_get_interior_diagnostic_data(model, jn) 
    336330      END DO 
    337331      DO jn = 1, jp_fabm_surface 
     
    344338         fabm_2d_25h(:,:,jn) = fabm_2d_25h(:,:,jn) + fabm_get_horizontal_diagnostic_data(model,jn) 
    345339      END DO 
    346       fabm_visib_25h(:,:,:) = fabm_visib_25h(:,:,:) + visib(:,:,:) 
    347340#endif 
    348341         cnt_25h = cnt_25h + 1 
     
    385378            fabm_bottom_25h(:,:,:)  = fabm_bottom_25h(:,:,:)  / 25.0_wp 
    386379            fabm_2d_25h(:,:,:)      = fabm_2d_25h(:,:,:)      / 25.0_wp 
    387             fabm_visib_25h(:,:,:)   = fabm_visib_25h(:,:,:)   / 25.0_wp 
    388380#endif 
    389381 
     
    427419               CALL iom_put( TRIM(model%horizontal_diagnostic_variables(jn)%name)//"25h", zw2d  ) 
    428420            END DO 
    429             zw3d(:,:,:) = fabm_visib_25h(:,:,:)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 
    430             CALL iom_put( "visib25h", zw3d  ) 
    431421#endif 
    432422 
     
    478468      END DO 
    479469      DO jn = 1, jp_fabm_3d 
    480          fabm_3d_25h(:,:,:,jn) = fabm_get_bulk_diagnostic_data(model, jn) 
     470         fabm_3d_25h(:,:,:,jn) = fabm_get_interior_diagnostic_data(model, jn) 
    481471      END DO 
    482472      DO jn = 1, jp_fabm_surface 
     
    489479         fabm_2d_25h(:,:,jn) = fabm_get_horizontal_diagnostic_data(model,jn) 
    490480      END DO 
    491       fabm_visib_25h(:,:,:) = visib(:,:,:) 
    492481#endif 
    493482            cnt_25h = 1 
Note: See TracChangeset for help on using the changeset viewer.