Changeset 9022


Ignore:
Timestamp:
2017-12-13T17:38:29+01:00 (4 years ago)
Author:
glong
Message:

Removed old dyn_vrt_dia subrouting and associated code. Updated trdvor to use new io and output vertical integrals as well as vertical averages.

Location:
branches/UKMO/GO6_dyn_vrt_diag/NEMOGCM
Files:
16 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/GO6_dyn_vrt_diag/NEMOGCM/CONFIG/SHARED/field_def.xml

    r8300 r9022  
    499499        <field id="tke"          long_name="Turbulent Kinetic Energy"                  unit="m2/s2"    /> 
    500500        <field id="tke_niw"      long_name="Turbulent Kinetic Energy due to near-inertial wave breaking"  unit="m2/s2" /> 
    501         <!-- variables available with dyn_vrt_dia - available with TODO FLAGS and TODO FLAGS --> 
    502         <field id="dia_vor_int-hpg"    long_name="vertical integral vorticity tendencies for hpg"  standard_name="vertically_integrated_vorticity_tendencies_hpg"   unit="s-2"      grid_ref="grid_W_2D" /> 
    503         <field id="dia_vor_mean-hpg"    long_name="vertical mean vorticity tendencies for hpg"  standard_name="vertical_mean_vorticity_tendencies_hpg"   unit="s-2"      grid_ref="grid_W_2D" /> 
    504         <field id="dia_vor_int-keg"    long_name="vertical integral vorticity tendencies for keg"  standard_name="vertically_integrated_vorticity_tendencies_keg"   unit="s-2"      grid_ref="grid_W_2D" /> 
    505         <field id="dia_vor_mean-keg"    long_name="vertical mean vorticity tendencies for keg"  standard_name="vertical_mean_vorticity_tendencies_keg"   unit="s-2"      grid_ref="grid_W_2D" /> 
    506         <field id="dia_vor_int-ldf"    long_name="vertical integral vorticity tendencies for ldf"  standard_name="vertically_integrated_vorticity_tendencies_ldf"   unit="s-2"      grid_ref="grid_W_2D" /> 
    507         <field id="dia_vor_mean-ldf"    long_name="vertical mean vorticity tendencies for ldf"  standard_name="vertical_mean_vorticity_tendencies_ldf"   unit="s-2"      grid_ref="grid_W_2D" /> 
    508         <field id="dia_vor_int-spg"    long_name="vertical integral vorticity tendencies for spg"  standard_name="vertically_integrated_vorticity_tendencies_spg"   unit="s-2"      grid_ref="grid_W_2D" /> 
    509         <field id="dia_vor_mean-spg"    long_name="vertical mean vorticity tendencies for spg"  standard_name="vertical_mean_vorticity_tendencies_spg"   unit="s-2"      grid_ref="grid_W_2D" /> 
    510         <field id="dia_vor_int-rvo"    long_name="vertical integral vorticity tendencies for rvo"  standard_name="vertically_integrated_vorticity_tendencies_rvo"   unit="s-2"      grid_ref="grid_W_2D" /> 
    511         <field id="dia_vor_mean-rvo"    long_name="vertical mean vorticity tendencies for rvo"  standard_name="vertical_mean_vorticity_tendencies_rvo"   unit="s-2"      grid_ref="grid_W_2D" /> 
    512         <field id="dia_vor_int-pvo"    long_name="vertical integral vorticity tendencies for pvo"  standard_name="vertically_integrated_vorticity_tendencies_pvo"   unit="s-2"      grid_ref="grid_W_2D" /> 
    513         <field id="dia_vor_mean-pvo"    long_name="vertical mean vorticity tendencies for pvo"  standard_name="vertical_mean_vorticity_tendencies_pvo"   unit="s-2"      grid_ref="grid_W_2D" /> 
    514         <field id="dia_vor_int-zad"    long_name="vertical integral vorticity tendencies for zad"  standard_name="vertically_integrated_vorticity_tendencies_zad"   unit="s-2"      grid_ref="grid_W_2D" /> 
    515         <field id="dia_vor_mean-zad"    long_name="vertical mean vorticity tendencies for zad"  standard_name="vertical_mean_vorticity_tendencies_zad"   unit="s-2"      grid_ref="grid_W_2D" /> 
    516         <field id="dia_vor_int-zdf"    long_name="vertical integral vorticity tendencies for zdf"  standard_name="vertically_integrated_vorticity_tendencies_zdf"   unit="s-2"      grid_ref="grid_W_2D" /> 
    517         <field id="dia_vor_mean-zdf"    long_name="vertical mean vorticity tendencies for zdf"  standard_name="vertical_mean_vorticity_tendencies_zdf"   unit="s-2"      grid_ref="grid_W_2D" /> 
    518         <field id="dia_vor_int-bfr"    long_name="vertical integral vorticity tendencies for bfr"  standard_name="vertically_integrated_vorticity_tendencies_bfr"   unit="s-2"      grid_ref="grid_W_2D" /> 
    519         <field id="dia_vor_mean-bfr"    long_name="vertical mean vorticity tendencies for bfr"  standard_name="vertical_mean_vorticity_tendencies_bfr"   unit="s-2"      grid_ref="grid_W_2D" /> 
    520         <field id="dia_vor_int-atf"    long_name="vertical integral vorticity tendencies for atf"  standard_name="vertically_integrated_vorticity_tendencies_atf"   unit="s-2"      grid_ref="grid_W_2D" /> 
    521         <field id="dia_vor_mean-atf"    long_name="vertical mean vorticity tendencies for atf"  standard_name="vertical_mean_vorticity_tendencies_atf"   unit="s-2"      grid_ref="grid_W_2D" /> 
     501        <!-- variables for trdvor output using ln_vor_trd --> 
     502        <!-- trdvor vertical average trends--> 
     503        <field id="sovortPh_avr"      long_name="averaged vertical average tendencies for pressure gradient trend"       standard_name="averaged_vor_grad_Ph"            unit="s-2"      grid_ref="grid_W_2D" /> 
     504        <field id="sovortEk_avr"      long_name="averaged vertical average tendencies for kinetic energy trend"          standard_name="averaged_vor_Energy"             unit="s-2"      grid_ref="grid_W_2D" /> 
     505        <field id="sovozeta_avr"      long_name="averaged vertical average tendencies for relative vorticity trend"      standard_name="averaged_vor_rel_vorticity"      unit="s-2"      grid_ref="grid_W_2D" /> 
     506        <field id="sovortif_avr"      long_name="averaged vertical average tendencies for coriolis trend"                standard_name="averaged_vor_coriolis"           unit="s-2"      grid_ref="grid_W_2D" /> 
     507        <field id="sovodifl_avr"      long_name="averaged vertical average tendencies for horizontal diffusion"          standard_name="averaged_vor_lat_diff"           unit="s-2"      grid_ref="grid_W_2D" /> 
     508        <field id="sovoadvv_avr"      long_name="averaged vertical average tendencies for vertical advection"            standard_name="averaged_vor_vert_adv"           unit="s-2"      grid_ref="grid_W_2D" /> 
     509        <field id="sovodifv_avr"      long_name="averaged vertical average tendencies for vertical diffusion"            standard_name="averaged_vor_vert_diff"          unit="s-2"      grid_ref="grid_W_2D" /> 
     510        <field id="sovortPs_avr"      long_name="averaged vertical average tendencies for surface pressure gradient"     standard_name="averaged_vor_grad_Ps"            unit="s-2"      grid_ref="grid_W_2D" /> 
     511        <field id="sovowind_avr"      long_name="averaged vertical average tendencies for wind stress"                   standard_name="averaged_vor_wind_stress"        unit="s-2"      grid_ref="grid_W_2D" /> 
     512        <field id="sovobfri_avr"      long_name="averaged vertical average tendencies for bottom friction"               standard_name="averaged_vor_bottom_friction"    unit="s-2"      grid_ref="grid_W_2D" /> 
     513        <field id="1st_mbre_avr"      long_name="averaged vertical average tendencies for 1st member of vrt eqn"         standard_name="averaged_vor_1st_mbre"           unit="s-2"      grid_ref="grid_W_2D" /> 
     514        <field id="sovorgap_avr"      long_name="averaged vertical average tendencies gap between 1st & 2nd member"      standard_name="averaged_vor_gap"                unit="s-2"      grid_ref="grid_W_2D" /> 
     515        <!-- trdvor vertical integral trends--> 
     516        <field id="sovortPh_int"      long_name="averaged vertical integral tendencies for pressure gradient trend"      standard_name="integrated_vor_grad_Ph"            unit="s-2"      grid_ref="grid_W_2D" /> 
     517        <field id="sovortEk_int"      long_name="averaged vertical integral tendencies for kinetic energy trend"         standard_name="integrated_vor_Energy"             unit="s-2"      grid_ref="grid_W_2D" /> 
     518        <field id="sovozeta_int"      long_name="averaged vertical integral tendencies for relative vorticity trend"     standard_name="integrated_vor_rel_vorticity"      unit="s-2"      grid_ref="grid_W_2D" /> 
     519        <field id="sovortif_int"      long_name="averaged vertical integral tendencies for coriolis trend"               standard_name="integrated_vor_coriolis"           unit="s-2"      grid_ref="grid_W_2D" /> 
     520        <field id="sovodifl_int"      long_name="averaged vertical integral tendencies for horizontal diffusion"         standard_name="integrated_vor_lat_diff"           unit="s-2"      grid_ref="grid_W_2D" /> 
     521        <field id="sovoadvv_int"      long_name="averaged vertical integral tendencies for vertical advection"           standard_name="integrated_vor_vert_adv"           unit="s-2"      grid_ref="grid_W_2D" /> 
     522        <field id="sovodifv_int"      long_name="averaged vertical integral tendencies for vertical diffusion"           standard_name="integrated_vor_vert_diff"          unit="s-2"      grid_ref="grid_W_2D" /> 
     523        <field id="sovortPs_int"      long_name="averaged vertical integral tendencies for surface pressure gradient"    standard_name="integrated_vor_grad_Ps"            unit="s-2"      grid_ref="grid_W_2D" /> 
     524        <field id="sovowind_int"      long_name="averaged vertical integral tendencies for wind stress"                  standard_name="integrated_vor_wind_stress"        unit="s-2"      grid_ref="grid_W_2D" /> 
     525        <field id="sovobfri_int"      long_name="averaged vertical integral tendencies for bottom friction"              standard_name="integrated_vor_bottom_friction"    unit="s-2"      grid_ref="grid_W_2D" /> 
     526        <field id="1st_mbre_int"      long_name="averaged vertical integral tendencies for 1st member of vrt eqn"        standard_name="integrated_vor_1st_mbre"           unit="s-2"      grid_ref="grid_W_2D" /> 
     527        <field id="sovorgap_int"      long_name="averaged vertical integral tendencies gap between 1st & 2nd member"     standard_name="integrated_vor_gap"                unit="s-2"      grid_ref="grid_W_2D" /> 
    522528      </field_group> 
    523529           
  • branches/UKMO/GO6_dyn_vrt_diag/NEMOGCM/NEMO/OPA_SRC/DYN/divcur.F90

    r8300 r9022  
    3535   USE wrk_nemo        ! Memory Allocation 
    3636   USE timing          ! Timing 
    37    USE iom             ! I/O Manager for dyn_vrt_dia_2d 
    3837 
    3938   IMPLICIT NONE 
    4039   PRIVATE 
    4140 
    42    PUBLIC   div_cur        ! routine called by step.F90 and istate.F90 
    43    PUBLIC   dyn_vrt_dia_3d ! routine called by various modules 
    44    PUBLIC   dyn_vrt_dia_2d ! routine called by various modules 
     41   PUBLIC   div_cur    ! routine called by step.F90 and istate.F90 
    4542 
    4643   !! * Substitutions 
     
    341338    
    342339#endif 
    343  
    344  
    345    SUBROUTINE dyn_vrt_dia_3d( utend, vtend, id_dia_vor) 
    346  
    347       !!---------------------------------------------------------------------- 
    348       !! 
    349       !! ** Purpose :  compute the vertical integrals of utend and vtend, and 
    350       !!               then pass to dyn_vrt_dia_2d to calculate vorticity 
    351       !!               tendencies. 
    352       !! 
    353       !! ** Action : a) Calculate the vertical integrals of utend & of vtend 
    354       !!                (u_int & v_int) 
    355       !!             b) Call dyn_vrt_dia_2d with vertical integrals 
    356       !! 
    357       !!---------------------------------------------------------------------- 
    358       REAL             :: utend(jpi,jpj,jpk) ! contribution to du/dt 
    359       REAL             :: vtend(jpi,jpj,jpk) ! contribution to dv/dt 
    360       CHARACTER(len=3) :: id_dia_vor         ! identifier for the diagnostic 
    361       ! 
    362       !!---------------------------------------------------------------------- 
    363       ! 
    364       INTEGER  ::   ji, jj, jk     ! dummy loop indices 
    365       ! 
    366       REAL(wp), POINTER, DIMENSION(:,:) :: u_int   ! u vertical integral 
    367       REAL(wp), POINTER, DIMENSION(:,:) :: v_int   ! v vertical integral 
    368  
    369       CALL wrk_alloc(jpi, jpj, u_int) 
    370       CALL wrk_alloc(jpi, jpj, v_int) 
    371  
    372       u_int(:,:) = 0.0_wp 
    373       v_int(:,:) = 0.0_wp 
    374  
    375       ! 
    376       ! Calculate the vertical integrals of utend & of vtend 
    377       ! 
    378       DO jk = 1, jpkm1 
    379           DO jj = 2, jpjm1 
    380               DO ji = fs_2, fs_jpim1 
    381                   u_int(ji,jj) = u_int(ji,jj) + ( utend(ji,jj,jk) * fse3u(ji,jj,jk) & 
    382                                  &                 * e1u(ji,jj) * umask(ji,jj,jk) ) 
    383                   v_int(ji,jj) = v_int(ji,jj) + ( vtend(ji,jj,jk) * fse3v(ji,jj,jk) & 
    384                                  &                 * e2v(ji,jj) * vmask(ji,jj,jk) ) 
    385               END DO 
    386           END DO 
    387       END DO 
    388  
    389       CALL dyn_vrt_dia_2d(u_int, v_int, id_dia_vor) 
    390  
    391       CALL wrk_dealloc(jpi, jpj, u_int) 
    392       CALL wrk_dealloc(jpi, jpj, v_int) 
    393  
    394    END SUBROUTINE dyn_vrt_dia_3d 
    395  
    396  
    397    SUBROUTINE dyn_vrt_dia_2d( u_int, v_int, id_dia_vor) 
    398  
    399       !!---------------------------------------------------------------------- 
    400       !! 
    401       !! ** Purpose :  compute the integral and mean vorticity tendencies. 
    402       !! 
    403       !! ** Action : a) Calculate the vorticity tendencies for the vertical 
    404       !!                integrals. 
    405       !!             b) Calculate the vertical means, u_mn, v_mn from u_int 
    406       !!                and v_int by dividing by the depth 
    407       !!             c) Calculate the vorticity tendencies for the vertical 
    408       !!                means 
    409       !!             d) Call iom_put for the vertical integral vorticity 
    410       !!                tendencies (using cid_dia_vor_int) 
    411       !!             e) Call iom_put for the vertical mean vorticity 
    412       !!                tendencies (using cid_dia_vor_mn) 
    413       !! 
    414       !!---------------------------------------------------------------------- 
    415       REAL             :: u_int(jpi,jpj)  ! u vertical integral 
    416       REAL             :: v_int(jpi,jpj)  ! v vertical integral 
    417       CHARACTER(len=3) :: id_dia_vor      ! identifier for the vorticity diagnostic 
    418       ! 
    419       !!---------------------------------------------------------------------- 
    420       ! 
    421       INTEGER  ::   ji, jj       ! dummy loop indices 
    422       INTEGER  ::   ikbu, ikbv   ! dummy loop indices 
    423       ! 
    424       REAL(wp), POINTER, DIMENSION(:,:) :: u_mn    ! u vertical means 
    425       REAL(wp), POINTER, DIMENSION(:,:) :: v_mn    ! u vertical means 
    426       REAL(wp), POINTER, DIMENSION(:,:) :: vor_int ! vort trend of vert integrals 
    427       REAL(wp), POINTER, DIMENSION(:,:) :: vor_mn  ! vort trend of vert means 
    428       CHARACTER(len=20) ::   cid_dia_vor_int       ! Identifier for iom_put 
    429       CHARACTER(len=20) ::   cid_dia_vor_mn        ! Identifier for iom_put 
    430  
    431       CALL wrk_alloc(jpi, jpj, u_mn) 
    432       CALL wrk_alloc(jpi, jpj, v_mn) 
    433       CALL wrk_alloc(jpi, jpj, vor_int) 
    434       CALL wrk_alloc(jpi, jpj, vor_mn) 
    435  
    436       CALL lbc_lnk( u_int, 'U', 1. ) 
    437       CALL lbc_lnk( v_int, 'V', 1. ) 
    438  
    439       WRITE ( cid_dia_vor_int, "(A16,A3)" ) "dia_vor_int-",  id_dia_vor 
    440       WRITE ( cid_dia_vor_mn,  "(A17,A3)" ) "dia_vor_mean-", id_dia_vor 
    441  
    442       ! 
    443       ! Calculate the vorticity tendencies for the vertical integrals. 
    444       ! 1/e1e2 * ((d(vtend)/dx) - (d(utend)/dy)) 
    445       ! 
    446  
    447       DO jj = 2, jpjm1 
    448           DO ji = fs_2, fs_jpim1 
    449               vor_int(ji,jj) = (   ( v_int(ji+1,jj) - v_int(ji,jj) )     & 
    450                   &              - ( u_int(ji,jj+1) - u_int(ji,jj) ) )   & 
    451                   &            / ( e1f(ji,jj)    * e2f(ji,jj) ) 
    452           END DO 
    453       END DO 
    454  
    455       ! Multiply by the surface mask 
    456       vor_int(:,:) = vor_int(:,:) * fmask(:,:,1) 
    457  
    458       ! 
    459       ! Calculate the vertical means, u_mn, v_mn from u_int & v_int by dividing 
    460       ! by the depth 
    461       ! mbku & mbkv - vertical index of the bottom last U- & W- ocean level 
    462       ! 
    463  
    464       DO jj = 2, jpjm1 
    465           DO ji = fs_2, fs_jpim1 
    466               ikbu = mbku(ji,jj) 
    467               ikbv = mbkv(ji,jj) 
    468  
    469               IF (gdepw_n(ji,jj,ikbu+1) .ne. 0.0_wp) THEN      ! Don't divide by 0! 
    470                   u_mn(ji,jj) = u_int(ji,jj) / gdepw_n(ji,jj,ikbu+1) 
    471               ELSE 
    472                   u_mn(ji,jj) = 0.0_wp 
    473               END IF 
    474  
    475               IF (gdepw_n(ji,jj,ikbv+1) .ne. 0.0_wp) THEN      ! Don't divide by 0! 
    476                   v_mn(ji,jj) = v_int(ji,jj) / gdepw_n(ji,jj,ikbv+1) 
    477               ELSE 
    478                   v_mn(ji,jj) = 0.0_wp 
    479               END IF 
    480           END DO 
    481       END DO 
    482  
    483       ! 
    484       ! Calculate the vorticity tendencies for the vertical means 
    485       ! 1/e1e2 * ((d(v_mn)/dx) - (d(u_mn)/dy)) 
    486       ! 
    487  
    488       DO jj = 2, jpjm1 
    489           DO ji = fs_2, fs_jpim1 
    490               vor_mn(ji,jj) = (   ( v_mn(ji+1,jj) - v_mn(ji,jj) )        & 
    491                   &             - ( u_mn(ji,jj+1) - u_mn(ji,jj) ) )      & 
    492                   &          / ( e1f(ji,jj)   * e2f(ji,jj) ) 
    493           END DO 
    494       END DO 
    495  
    496       ! Multiply by the surface mask 
    497       vor_mn(:,:) = vor_mn(:,:) * fmask(:,:,1) 
    498  
    499       ! Call iom_put for the vertical integral vorticity tendencies 
    500       CALL iom_put( cid_dia_vor_int, vor_int(:,:)) 
    501  
    502       ! Call iom_put for the vertical mean vorticity tendencies 
    503       CALL iom_put( cid_dia_vor_mn, vor_mn(:,:)) 
    504  
    505       CALL wrk_dealloc(jpi, jpj, u_mn) 
    506       CALL wrk_dealloc(jpi, jpj, v_mn) 
    507       CALL wrk_dealloc(jpi, jpj, vor_int) 
    508       CALL wrk_dealloc(jpi, jpj, vor_mn) 
    509  
    510    END SUBROUTINE dyn_vrt_dia_2d 
    511  
    512  
    513340   !!====================================================================== 
    514341END MODULE divcur 
  • branches/UKMO/GO6_dyn_vrt_diag/NEMOGCM/NEMO/OPA_SRC/DYN/dynbfr.F90

    r8300 r9022  
    2222   USE timing         ! Timing 
    2323   USE wrk_nemo       ! Memory Allocation 
    24    USE divcur         ! for dyn_vrt_dia 
    2524 
    2625   IMPLICIT NONE 
     
    5352      INTEGER  ::   ikbu, ikbv   ! local integers 
    5453      REAL(wp) ::   zm1_2dt      ! local scalar 
    55       CHARACTER(len=3) ::  id_dyn_vrt_bfr = "bfr"    ! TODO remove once flags set properly 
    5654      REAL(wp), POINTER, DIMENSION(:,:,:) ::  ztrdu, ztrdv 
    5755      !!--------------------------------------------------------------------- 
     
    6765        zm1_2dt = - 1._wp / ( 2._wp * rdt ) 
    6866 
    69         IF( l_trddyn .OR. ( id_dyn_vrt_bfr == "bfr" ) )   THEN                      ! temporary save of ua and va trends 
     67        IF( l_trddyn )   THEN                      ! temporary save of ua and va trends 
    7068           CALL wrk_alloc( jpi,jpj,jpk, ztrdu, ztrdv ) 
    7169           ztrdu(:,:,:) = ua(:,:,:) 
     
    103101 
    104102        ! 
    105         IF( l_trddyn .OR. ( id_dyn_vrt_bfr == "bfr" ) )   THEN                      ! save the vertical diffusive trends for further diagnostics 
     103        IF( l_trddyn )   THEN                      ! save the vertical diffusive trends for further diagnostics 
    106104           ztrdu(:,:,:) = ua(:,:,:) - ztrdu(:,:,:) 
    107105           ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:) 
    108            IF( l_trddyn )  THEN 
    109               CALL trd_dyn( ztrdu(:,:,:), ztrdv(:,:,:), jpdyn_bfr, kt ) 
    110            ENDIF 
    111            IF( id_dyn_vrt_bfr == "bfr" )  THEN 
    112               CALL dyn_vrt_dia_3d( ztrdu, ztrdv, id_dyn_vrt_bfr ) 
    113            ENDIF 
     106           CALL trd_dyn( ztrdu(:,:,:), ztrdv(:,:,:), jpdyn_bfr, kt ) 
    114107           CALL wrk_dealloc( jpi,jpj,jpk, ztrdu, ztrdv ) 
    115108        ENDIF 
  • branches/UKMO/GO6_dyn_vrt_diag/NEMOGCM/NEMO/OPA_SRC/DYN/dynhpg.F90

    r8300 r9022  
    3636   USE trd_oce         ! trends: ocean variables 
    3737   USE trddyn          ! trend manager: dynamics 
    38    USE divcur          ! for dyn_vrt_dia_3d 
    3938   ! 
    4039   USE in_out_manager  ! I/O manager 
     
    8382      !!             - send trends to trd_dyn for futher diagnostics (l_trddyn=T) 
    8483      !!---------------------------------------------------------------------- 
    85       INTEGER, INTENT(in) ::   kt                     ! ocean time-step index 
    86       CHARACTER(len=3)    ::   id_vrt_dia_hpg = "hpg" ! TODO remove once flags set properly 
     84      INTEGER, INTENT(in) ::   kt   ! ocean time-step index 
    8785      REAL(wp), POINTER, DIMENSION(:,:,:) ::  ztrdu, ztrdv 
    8886      !!---------------------------------------------------------------------- 
     
    9088      IF( nn_timing == 1 )  CALL timing_start('dyn_hpg') 
    9189      ! 
    92       IF( l_trddyn .or. ( id_vrt_dia_hpg == "hpg" ) ) THEN    ! Temporary saving of ua and va trends (l_trddyn) 
     90      IF( l_trddyn ) THEN                    ! Temporary saving of ua and va trends (l_trddyn) 
    9391         CALL wrk_alloc( jpi,jpj,jpk, ztrdu, ztrdv ) 
    9492         ztrdu(:,:,:) = ua(:,:,:) 
     
    105103      END SELECT 
    106104      ! 
    107       IF( l_trddyn .or. ( id_vrt_dia_hpg == "hpg" ) ) THEN      ! save the hydrostatic pressure gradient trends for momentum trend diagnostics 
     105      IF( l_trddyn ) THEN      ! save the hydrostatic pressure gradient trends for momentum trend diagnostics 
    108106         ztrdu(:,:,:) = ua(:,:,:) - ztrdu(:,:,:) 
    109107         ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:) 
    110          ! 
    111          IF( id_vrt_dia_hpg == "hpg" ) THEN 
    112              CALL dyn_vrt_dia_3d( ztrdu, ztrdv, id_vrt_dia_hpg ) 
    113          END IF 
    114          IF( l_trddyn ) THEN 
    115              CALL trd_dyn( ztrdu, ztrdv, jpdyn_hpg, kt ) 
    116          END IF 
     108         CALL trd_dyn( ztrdu, ztrdv, jpdyn_hpg, kt ) 
    117109         CALL wrk_dealloc( jpi,jpj,jpk, ztrdu, ztrdv ) 
    118110      ENDIF 
  • branches/UKMO/GO6_dyn_vrt_diag/NEMOGCM/NEMO/OPA_SRC/DYN/dynkeg.F90

    r8300 r9022  
    1717   USE trd_oce         ! trends: ocean variables 
    1818   USE trddyn          ! trend manager: dynamics 
    19    USE divcur          ! for dyn_vrt_dia_3d 
    2019   ! 
    2120   USE in_out_manager  ! I/O manager 
     
    7574      INTEGER, INTENT( in ) ::   kscheme   ! =0/1   type of KEG scheme  
    7675      ! 
    77       CHARACTER(len=3) :: id_vrt_dia_keg = "keg" ! TODO remove once flags set properly 
    78       INTEGER  ::   ji, jj, jk             ! dummy loop indices 
    79       REAL(wp) ::   zu, zv                 ! temporary scalars 
     76      INTEGER  ::   ji, jj, jk   ! dummy loop indices 
     77      REAL(wp) ::   zu, zv       ! temporary scalars 
    8078      REAL(wp), POINTER, DIMENSION(:,:,:) :: zhke 
    8179      REAL(wp), POINTER, DIMENSION(:,:,:) :: ztrdu, ztrdv  
     
    9290      ENDIF 
    9391 
    94       IF( l_trddyn .or. ( id_vrt_dia_keg == "keg" ) ) THEN           ! Save ua and va trends 
     92      IF( l_trddyn ) THEN           ! Save ua and va trends 
    9593         CALL wrk_alloc( jpi,jpj,jpk,   ztrdu, ztrdv ) 
    9694         ztrdu(:,:,:) = ua(:,:,:)  
     
    145143      END DO 
    146144      ! 
    147       IF( l_trddyn .or. ( id_vrt_dia_keg == "keg" ) ) THEN      ! save the Kinetic Energy trends for diagnostic 
     145      IF( l_trddyn ) THEN                 ! save the Kinetic Energy trends for diagnostic 
    148146         ztrdu(:,:,:) = ua(:,:,:) - ztrdu(:,:,:) 
    149147         ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:) 
    150          ! 
    151          IF( id_vrt_dia_keg == "keg" ) THEN 
    152              CALL dyn_vrt_dia_3d( ztrdu, ztrdv, id_vrt_dia_keg ) 
    153          END IF 
    154          ! 
    155          IF( l_trddyn ) THEN 
    156              CALL trd_dyn( ztrdu, ztrdv, jpdyn_keg, kt ) 
    157          END IF 
    158          ! 
     148         CALL trd_dyn( ztrdu, ztrdv, jpdyn_keg, kt ) 
    159149         CALL wrk_dealloc( jpi,jpj,jpk,   ztrdu, ztrdv ) 
    160150      ENDIF 
  • branches/UKMO/GO6_dyn_vrt_diag/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf.F90

    r8300 r9022  
    3030   USE wrk_nemo       ! Memory Allocation 
    3131   USE timing         ! Timing 
    32    USE divcur         ! Used by dyn_vrt_dia 
    3332 
    3433   IMPLICIT NONE 
     
    5857      INTEGER, INTENT(in) ::   kt   ! ocean time-step index 
    5958      ! 
    60       CHARACTER(len=3)    ::   id_vrt_dia_ldf = "ldf"  ! TODO Replace once proper flags in place 
    6159      REAL(wp), POINTER, DIMENSION(:,:,:) ::  ztrdu, ztrdv 
    6260      !!---------------------------------------------------------------------- 
     
    6462      IF( nn_timing == 1 )  CALL timing_start('dyn_ldf') 
    6563      ! 
    66       IF( l_trddyn .or. ( id_vrt_dia_ldf == "ldf" ) )   THEN               ! temporary save of ta and sa trends 
     64      IF( l_trddyn )   THEN                      ! temporary save of ta and sa trends 
    6765         CALL wrk_alloc( jpi, jpj, jpk, ztrdu, ztrdv ) 
    6866         ztrdu(:,:,:) = ua(:,:,:)  
     
    105103      END SELECT 
    106104 
    107       IF( l_trddyn .or. ( id_vrt_dia_ldf == "ldf" ) ) THEN    ! save the horizontal diffusive trends for further diagnostics 
     105      IF( l_trddyn ) THEN                        ! save the horizontal diffusive trends for further diagnostics 
    108106         ztrdu(:,:,:) = ua(:,:,:) - ztrdu(:,:,:) 
    109107         ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:) 
    110          IF( id_vrt_dia_ldf == "ldf" ) THEN 
    111              CALL dyn_vrt_dia_3d( ztrdu, ztrdv, id_vrt_dia_ldf ) 
    112          END IF 
    113          IF( l_trddyn ) THEN 
    114              CALL trd_dyn( ztrdu, ztrdv, jpdyn_ldf, kt ) 
    115          END IF 
     108         CALL trd_dyn( ztrdu, ztrdv, jpdyn_ldf, kt ) 
    116109         CALL wrk_dealloc( jpi, jpj, jpk, ztrdu, ztrdv ) 
    117110      ENDIF 
  • branches/UKMO/GO6_dyn_vrt_diag/NEMOGCM/NEMO/OPA_SRC/DYN/dynnxt.F90

    r8300 r9022  
    4646   USE prtctl          ! Print control 
    4747   USE timing          ! Timing 
    48    USE divcur          ! for dyn_vrt_dia 
    4948#if defined key_agrif 
    5049   USE agrif_opa_interp 
     
    107106      REAL(wp) ::   zue3a, zue3n, zue3b, zuf, zec      ! local scalars 
    108107      REAL(wp) ::   zve3a, zve3n, zve3b, zvf, z1_2dt   !   -      - 
    109       CHARACTER(len=3) ::  id_dyn_vrt_atf = "atf"      ! TODO remove once flags done 
    110108      REAL(wp), POINTER, DIMENSION(:,:)   ::  zue, zve 
    111109      REAL(wp), POINTER, DIMENSION(:,:,:) ::  ze3u_f, ze3v_f, zua, zva  
     
    205203#endif 
    206204 
    207       IF( l_trddyn .OR. ( id_dyn_vrt_atf == "atf" ) ) THEN             ! prepare the atf trend computation + some diagnostics 
     205      IF( l_trddyn ) THEN             ! prepare the atf trend computation + some diagnostics 
    208206         z1_2dt = 1._wp / (2. * rdt)        ! Euler or leap-frog time step  
    209          IF( l_trddyn )  THEN 
    210             IF( neuler == 0 .AND. kt == nit000 )   z1_2dt = 1._wp / rdt 
    211             ! 
    212             !                                  ! Kinetic energy and Conversion 
    213             IF( ln_KE_trd  )   CALL trd_dyn( ua, va, jpdyn_ken, kt ) 
    214             ! 
    215             IF( ln_dyn_trd ) THEN              ! 3D output: total momentum trends 
    216                zua(:,:,:) = ( ua(:,:,:) - ub(:,:,:) ) * z1_2dt 
    217                zva(:,:,:) = ( va(:,:,:) - vb(:,:,:) ) * z1_2dt 
    218                CALL iom_put( "utrd_tot", zua )        ! total momentum trends, except the asselin time filter 
    219                CALL iom_put( "vtrd_tot", zva ) 
    220             ENDIF 
     207         IF( neuler == 0 .AND. kt == nit000 )   z1_2dt = 1._wp / rdt 
     208         ! 
     209         !                                  ! Kinetic energy and Conversion 
     210         IF( ln_KE_trd  )   CALL trd_dyn( ua, va, jpdyn_ken, kt ) 
     211         ! 
     212         IF( ln_dyn_trd ) THEN              ! 3D output: total momentum trends 
     213            zua(:,:,:) = ( ua(:,:,:) - ub(:,:,:) ) * z1_2dt 
     214            zva(:,:,:) = ( va(:,:,:) - vb(:,:,:) ) * z1_2dt 
     215            CALL iom_put( "utrd_tot", zua )        ! total momentum trends, except the asselin time filter 
     216            CALL iom_put( "vtrd_tot", zva ) 
    221217         ENDIF 
    222218         ! 
     
    396392      ! 
    397393 
    398       IF( l_trddyn .OR. ( id_dyn_vrt_atf == "atf" ) ) THEN                ! 3D output: asselin filter trends on momentum 
     394      IF( l_trddyn ) THEN                ! 3D output: asselin filter trends on momentum 
    399395         zua(:,:,:) = ( ub(:,:,:) - zua(:,:,:) ) * z1_2dt 
    400396         zva(:,:,:) = ( vb(:,:,:) - zva(:,:,:) ) * z1_2dt 
    401          IF( id_dyn_vrt_atf == "atf" ) THEN                ! 3D output: asselin filter trends on momentum 
    402             CALL dyn_vrt_dia_3d( zua, zva, id_dyn_vrt_atf ) 
    403          ENDIF 
    404          IF( l_trddyn ) THEN                ! 3D output: asselin filter trends on momentum 
    405             CALL trd_dyn( zua, zva, jpdyn_atf, kt ) 
    406          ENDIF 
     397         CALL trd_dyn( zua, zva, jpdyn_atf, kt ) 
    407398      ENDIF 
    408399      ! 
  • branches/UKMO/GO6_dyn_vrt_diag/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg.F90

    r8300 r9022  
    3535   USE wrk_nemo       ! Memory Allocation 
    3636   USE timing         ! Timing 
    37    USE divcur         ! for dyn_vrt_dia_3d 
    3837 
    3938 
     
    9089      REAL(wp), POINTER, DIMENSION(:,:,:) ::  ztrdu, ztrdv 
    9190      REAL(wp), POINTER, DIMENSION(:,:)   ::  zpice 
    92       CHARACTER(len=3) ::  id_vrt_dia_spg = "spg"          ! TODO remove once flags set properly 
    9391      !!---------------------------------------------------------------------- 
    9492      ! 
     
    10199 
    102100 
    103       IF( l_trddyn .or. ( id_vrt_dia_spg == "spg" ) )   THEN      ! temporary save of ta and sa trends 
    104          CALL wrk_alloc( jpi, jpj, jpk, ztrdu, ztrdv )  
     101      IF( l_trddyn )   THEN                      ! temporary save of ta and sa trends 
     102         CALL wrk_alloc( jpi, jpj, jpk, ztrdu, ztrdv ) 
    105103         ztrdu(:,:,:) = ua(:,:,:) 
    106104         ztrdv(:,:,:) = va(:,:,:) 
     
    190188      END SELECT 
    191189      !                     
    192       IF( l_trddyn .or. ( id_vrt_dia_spg == "spg" ) )   THEN      ! save the surface pressure gradient trends for further diagnostics 
     190      IF( l_trddyn )   THEN                      ! save the surface pressure gradient trends for further diagnostics 
    193191         SELECT CASE ( nspg ) 
    194192         CASE ( 0, 1 ) 
     
    201199            ztrdv(:,:,:) = ( va(:,:,:) - vb(:,:,:) ) / z2dt - ztrdv(:,:,:) 
    202200         END SELECT 
    203          ! 
    204          IF( id_vrt_dia_spg == "spg" ) THEN 
    205              CALL dyn_vrt_dia_3d( ztrdu, ztrdv, id_vrt_dia_spg ) 
    206          END IF 
    207          ! 
    208          IF( l_trddyn ) THEN 
    209              CALL trd_dyn( ztrdu, ztrdv, jpdyn_spg, kt ) 
    210          END IF 
     201         CALL trd_dyn( ztrdu, ztrdv, jpdyn_spg, kt ) 
    211202         ! 
    212203         CALL wrk_dealloc( jpi, jpj, jpk, ztrdu, ztrdv )  
  • branches/UKMO/GO6_dyn_vrt_diag/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_flt.F90

    r8300 r9022  
    137137      z2dtg  = grav * z2dt 
    138138 
     139      IF( l_trddyn )   CALL wrk_alloc( jpi, jpj, jpk, ztrdu, ztrdv ) ! temporary save of spg trends 
     140 
    139141      ! Evaluate the masked next velocity (effect of the additional force not included) 
    140142      ! ---------------------------------   
     
    185187         ! 
    186188         IF( l_trddyn )   THEN                      ! temporary save of spg trends 
    187             CALL wrk_alloc( jpi, jpj, jpk, ztrdu, ztrdv ) 
    188189            DO jk = 1, jpkm1              ! unweighted time stepping  
    189190               DO jj = 2, jpjm1 
  • branches/UKMO/GO6_dyn_vrt_diag/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_ts.F90

    r7900 r9022  
    4141   USE timing          ! Timing     
    4242   USE sbcapr          ! surface boundary condition: atmospheric pressure 
    43    USE divcur          ! for dyn_vrt_dia_2d 
    4443   USE dynadv, ONLY: ln_dynadv_vec 
    4544#if defined key_agrif 
     
    145144      INTEGER  ::   ji, jj, jk, jn        ! dummy loop indices 
    146145      INTEGER  ::   ikbu, ikbv, noffset      ! local integers 
    147       INTEGER  ::   id_dia_vrt_spg_int  = 1  ! TODO remove once flags set properly 
    148       INTEGER  ::   id_dia_vrt_spg_mean = 1  ! TODO remove once flags set properly 
    149146      REAL(wp) ::   zraur, z1_2dt_b, z2dt_bf    ! local scalars 
    150147      REAL(wp) ::   zx1, zy1, zx2, zy2         !   -      - 
    151148      REAL(wp) ::   z1_12, z1_8, z1_4, z1_2    !   -      - 
     149      REAL(wp) ::   zu_spg, zv_spg             !   -      - 
    152150      REAL(wp) ::   zhura, zhvra               !   -      - 
    153151      REAL(wp) ::   za0, za1, za2, za3           !   -      - 
     
    158156      REAL(wp), POINTER, DIMENSION(:,:) :: zhup2_e, zhvp2_e, zhust_e, zhvst_e 
    159157      REAL(wp), POINTER, DIMENSION(:,:) :: zsshu_a, zsshv_a 
    160       REAL(wp), POINTER, DIMENSION(:,:) :: zu_dia, zv_dia, zu_wfrc, zv_wfrc 
    161       REAL(wp), POINTER, DIMENSION(:,:) :: zu_spg, zv_spg 
    162158      REAL(wp), POINTER, DIMENSION(:,:) :: zhf 
    163159      !!---------------------------------------------------------------------- 
     
    171167      CALL wrk_alloc( jpi, jpj, zhup2_e, zhvp2_e, zhust_e, zhvst_e) 
    172168      CALL wrk_alloc( jpi, jpj, zsshu_a, zsshv_a                                   ) 
    173       CALL wrk_alloc( jpi, jpj, zu_dia, zv_dia, zu_wfrc, zv_wfrc                   ) 
    174       CALL wrk_alloc( jpi, jpj, zu_spg, zv_spg                                     ) 
    175169      CALL wrk_alloc( jpi, jpj, zhf ) 
    176170      ! 
     
    400394      END DO  
    401395      ! 
    402       !                                         ! Set zu_dia for use in dyn_vrt_dia 
    403       zu_dia(:,:) = zu_frc(:,:) 
    404       zv_dia(:,:) = zv_frc(:,:) 
    405       ! 
    406396      !                 ! Add bottom stress contribution from baroclinic velocities:       
    407397      IF (ln_bt_fw) THEN 
     
    430420      !        
    431421      IF (ln_bt_fw) THEN                        ! Add wind forcing 
    432          zu_wfrc(:,:) = zraur * utau(:,:) * hur(:,:) 
    433          zv_wfrc(:,:) = zraur * vtau(:,:) * hvr(:,:) 
     422         zu_frc(:,:) =  zu_frc(:,:) + zraur * utau(:,:) * hur(:,:) 
     423         zv_frc(:,:) =  zv_frc(:,:) + zraur * vtau(:,:) * hvr(:,:) 
    434424      ELSE 
    435          zu_wfrc(:,:) = zraur * z1_2 * ( utau_b(:,:) + utau(:,:) ) * hur(:,:) 
    436          zv_wfrc(:,:) = zraur * z1_2 * ( vtau_b(:,:) + vtau(:,:) ) * hvr(:,:) 
     425         zu_frc(:,:) =  zu_frc(:,:) + zraur * z1_2 * ( utau_b(:,:) + utau(:,:) ) * hur(:,:) 
     426         zv_frc(:,:) =  zv_frc(:,:) + zraur * z1_2 * ( vtau_b(:,:) + vtau(:,:) ) * hvr(:,:) 
    437427      ENDIF   
    438       zu_frc(:,:) = zu_frc(:,:) + zu_wfrc(:,:) 
    439       zv_frc(:,:) = zv_frc(:,:) + zv_wfrc(:,:) 
    440  
    441                                         ! Calculate dia_vor_int & dia_vor_mn if required 
    442       IF ( ( id_dia_vrt_spg_int == 1 ) .or. (id_dia_vrt_spg_mean == 1 ) ) THEN 
    443          zu_dia(:,:) = zu_dia(:,:) + zu_wfrc(:,:) 
    444          zv_dia(:,:) = zv_dia(:,:) + zv_wfrc(:,:) 
    445          CALL dyn_vrt_dia_2d(zu_dia, zv_dia, id_dia_vrt_spg_int, id_dia_vrt_spg_mean) 
    446          !     Reset zu_dia for use again later. 
    447          zu_dia(:,:) = zu_dia(:,:) - zu_wfrc(:,:) 
    448          zv_dia(:,:) = zv_dia(:,:) - zv_wfrc(:,:) 
    449       END IF 
    450       ! 
    451428      ! 
    452429      IF ( ln_apr_dyn ) THEN                    ! Add atm pressure forcing 
     
    454431            DO jj = 2, jpjm1               
    455432               DO ji = fs_2, fs_jpim1   ! vector opt. 
    456                   zu_spg(ji,jj) =  grav * (  ssh_ib (ji+1,jj  ) - ssh_ib (ji,jj) ) /e1u(ji,jj) 
    457                   zv_spg(ji,jj) =  grav * (  ssh_ib (ji  ,jj+1) - ssh_ib (ji,jj) ) /e2v(ji,jj) 
     433                  zu_spg =  grav * (  ssh_ib (ji+1,jj  ) - ssh_ib (ji,jj) ) /e1u(ji,jj) 
     434                  zv_spg =  grav * (  ssh_ib (ji  ,jj+1) - ssh_ib (ji,jj) ) /e2v(ji,jj) 
     435                  zu_frc(ji,jj) = zu_frc(ji,jj) + zu_spg 
     436                  zv_frc(ji,jj) = zv_frc(ji,jj) + zv_spg 
    458437               END DO 
    459438            END DO 
     
    461440            DO jj = 2, jpjm1               
    462441               DO ji = fs_2, fs_jpim1   ! vector opt. 
    463                   zu_spg(ji,jj) =  grav * z1_2 * (  ssh_ib (ji+1,jj  ) - ssh_ib (ji,jj)    & 
     442                  zu_spg =  grav * z1_2 * (  ssh_ib (ji+1,jj  ) - ssh_ib (ji,jj)    & 
    464443                      &                    + ssh_ibb(ji+1,jj  ) - ssh_ibb(ji,jj)  ) /e1u(ji,jj) 
    465                   zv_spg(ji,jj) =  grav * z1_2 * (  ssh_ib (ji  ,jj+1) - ssh_ib (ji,jj)    & 
     444                  zv_spg =  grav * z1_2 * (  ssh_ib (ji  ,jj+1) - ssh_ib (ji,jj)    & 
    466445                      &                    + ssh_ibb(ji  ,jj+1) - ssh_ibb(ji,jj)  ) /e2v(ji,jj) 
     446                  zu_frc(ji,jj) = zu_frc(ji,jj) + zu_spg 
     447                  zv_frc(ji,jj) = zv_frc(ji,jj) + zv_spg 
    467448               END DO 
    468449            END DO 
    469450         ENDIF  
    470          zu_frc(:,:) = zu_frc(:,:) + zu_spg(:,:) 
    471          zv_frc(:,:) = zv_frc(:,:) + zv_spg(:,:) 
    472                                         ! Calculate dia_vor_int & dia_vor_mn if required 
    473          IF ( ( id_dia_vrt_spg_int == 1 ) .or. (id_dia_vrt_spg_mean == 1 ) ) THEN 
    474             zu_dia(:,:) = zu_dia(:,:) + zu_spg(:,:) 
    475             zv_dia(:,:) = zv_dia(:,:) + zv_spg(:,:) 
    476             CALL dyn_vrt_dia_2d(zu_dia, zv_dia, id_dia_vrt_spg_int, id_dia_vrt_spg_mean) 
    477          END IF 
    478       ENDIF 
    479       ! 
     451      ENDIF 
    480452      !                                   !* Right-Hand-Side of the barotropic ssh equation 
    481453      !                                   ! ----------------------------------------------- 
     
    952924      CALL wrk_dealloc( jpi, jpj, zhup2_e, zhvp2_e, zhust_e, zhvst_e ) 
    953925      CALL wrk_dealloc( jpi, jpj, zsshu_a, zsshv_a                                   ) 
    954       CALL wrk_dealloc( jpi, jpj, zu_dia, zv_dia, zu_wfrc, zv_wfrc                   ) 
    955       CALL wrk_dealloc( jpi, jpj, zu_spg, zv_spg                                     ) 
    956926      CALL wrk_dealloc( jpi, jpj, zhf ) 
    957927      ! 
  • branches/UKMO/GO6_dyn_vrt_diag/NEMOGCM/NEMO/OPA_SRC/DYN/dynvor.F90

    r8300 r9022  
    3838   USE wrk_nemo       ! Memory Allocation 
    3939   USE timing         ! Timing 
    40    USE divcur         ! For dyn_vrt_dia_3d 
    4140 
    4241 
     
    8180      INTEGER, INTENT( in ) ::   kt   ! ocean time-step index 
    8281      ! 
    83       CHARACTER(len=3)      ::   id_vrt_dia_rvo = "rvo"      ! TODO remove once flags set properly 
    84       CHARACTER(len=3)      ::   id_vrt_dia_pvo = "pvo"      ! TODO remove once flags set properly 
    8582      REAL(wp), POINTER, DIMENSION(:,:,:) ::  ztrdu, ztrdv 
    8683      !!---------------------------------------------------------------------- 
     
    8885      IF( nn_timing == 1 )  CALL timing_start('dyn_vor') 
    8986      ! 
    90       IF( l_trddyn .or. ( id_vrt_dia_rvo == "rvo" ) )   CALL wrk_alloc( jpi,jpj,jpk, ztrdu, ztrdv ) 
     87      IF( l_trddyn )   CALL wrk_alloc( jpi,jpj,jpk, ztrdu, ztrdv ) 
    9188      ! 
    9289      !                                          ! vorticity term  
     
    108105         ! 
    109106      CASE ( 0 )                                       ! energy conserving scheme 
    110          IF( l_trddyn .or. ( id_vrt_dia_rvo == "rvo" ) )   THEN 
     107         IF( l_trddyn )   THEN 
    111108            ztrdu(:,:,:) = ua(:,:,:) 
    112109            ztrdv(:,:,:) = va(:,:,:) 
     
    114111            ztrdu(:,:,:) = ua(:,:,:) - ztrdu(:,:,:) 
    115112            ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:) 
    116             IF( l_trddyn ) THEN 
    117                 CALL trd_dyn( ztrdu, ztrdv, jpdyn_rvo, kt ) 
    118             END IF 
    119             IF( id_vrt_dia_rvo == "rvo" ) THEN 
    120                 CALL dyn_vrt_dia_3d( ztrdu, ztrdv, id_vrt_dia_rvo ) 
    121             END IF 
     113            CALL trd_dyn( ztrdu, ztrdv, jpdyn_rvo, kt ) 
    122114            ztrdu(:,:,:) = ua(:,:,:) 
    123115            ztrdv(:,:,:) = va(:,:,:) 
     
    125117            ztrdu(:,:,:) = ua(:,:,:) - ztrdu(:,:,:) 
    126118            ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:) 
    127             IF( l_trddyn ) THEN 
    128                 CALL trd_dyn( ztrdu, ztrdv, jpdyn_pvo, kt ) 
    129             END IF 
    130             IF( id_vrt_dia_pvo == "pvo" ) THEN 
    131                 CALL dyn_vrt_dia_3d( ztrdu, ztrdv, id_vrt_dia_pvo ) 
    132             END IF 
     119            CALL trd_dyn( ztrdu, ztrdv, jpdyn_pvo, kt ) 
    133120         ELSE 
    134121            CALL vor_ene( kt, ntot, ua, va )                ! total vorticity 
     
    136123         ! 
    137124      CASE ( 1 )                                       ! enstrophy conserving scheme 
    138          IF( l_trddyn .or. ( id_vrt_dia_rvo == "rvo" ) )   THEN 
     125         IF( l_trddyn )   THEN     
    139126            ztrdu(:,:,:) = ua(:,:,:) 
    140127            ztrdv(:,:,:) = va(:,:,:) 
     
    142129            ztrdu(:,:,:) = ua(:,:,:) - ztrdu(:,:,:) 
    143130            ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:) 
    144             IF( l_trddyn ) THEN 
    145                 CALL trd_dyn( ztrdu, ztrdv, jpdyn_rvo, kt ) 
    146             END IF 
    147             IF( id_vrt_dia_rvo == "rvo" ) THEN 
    148                 CALL dyn_vrt_dia_3d( ztrdu, ztrdv, id_vrt_dia_rvo ) 
    149             END IF 
     131            CALL trd_dyn( ztrdu, ztrdv, jpdyn_rvo, kt ) 
    150132            ztrdu(:,:,:) = ua(:,:,:) 
    151133            ztrdv(:,:,:) = va(:,:,:) 
     
    153135            ztrdu(:,:,:) = ua(:,:,:) - ztrdu(:,:,:) 
    154136            ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:) 
    155             IF( l_trddyn ) THEN 
    156                 CALL trd_dyn( ztrdu, ztrdv, jpdyn_pvo, kt ) 
    157             END IF 
    158             IF( id_vrt_dia_pvo == "pvo" ) THEN 
    159                 CALL dyn_vrt_dia_3d( ztrdu, ztrdv, id_vrt_dia_pvo ) 
    160             END IF 
     137            CALL trd_dyn( ztrdu, ztrdv, jpdyn_pvo, kt ) 
    161138         ELSE 
    162139            CALL vor_ens( kt, ntot, ua, va )                ! total vorticity 
     
    164141         ! 
    165142      CASE ( 2 )                                       ! mixed ene-ens scheme 
    166          IF( l_trddyn .or. ( id_vrt_dia_rvo == "rvo" ) )   THEN 
     143         IF( l_trddyn )   THEN 
    167144            ztrdu(:,:,:) = ua(:,:,:) 
    168145            ztrdv(:,:,:) = va(:,:,:) 
     
    170147            ztrdu(:,:,:) = ua(:,:,:) - ztrdu(:,:,:) 
    171148            ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:) 
    172             IF( l_trddyn ) THEN 
    173                 CALL trd_dyn( ztrdu, ztrdv, jpdyn_rvo, kt ) 
    174             END IF 
    175             IF( id_vrt_dia_rvo == "rvo" ) THEN 
    176                 CALL dyn_vrt_dia_3d( ztrdu, ztrdv, id_vrt_dia_rvo ) 
    177             END IF 
     149            CALL trd_dyn( ztrdu, ztrdv, jpdyn_rvo, kt ) 
    178150            ztrdu(:,:,:) = ua(:,:,:) 
    179151            ztrdv(:,:,:) = va(:,:,:) 
     
    181153            ztrdu(:,:,:) = ua(:,:,:) - ztrdu(:,:,:) 
    182154            ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:) 
    183             IF( l_trddyn ) THEN 
    184                 CALL trd_dyn( ztrdu, ztrdv, jpdyn_pvo, kt ) 
    185             END IF 
    186             IF( id_vrt_dia_pvo == "pvo" ) THEN 
    187                 CALL dyn_vrt_dia_3d( ztrdu, ztrdv, id_vrt_dia_pvo ) 
    188             END IF 
     155            CALL trd_dyn( ztrdu, ztrdv, jpdyn_pvo, kt ) 
    189156         ELSE 
    190157            CALL vor_mix( kt )                               ! total vorticity (mix=ens-ene) 
     
    192159         ! 
    193160      CASE ( 3 )                                       ! energy and enstrophy conserving scheme 
    194          IF( l_trddyn .or. ( id_vrt_dia_rvo == "rvo" ) )   THEN 
     161         IF( l_trddyn )   THEN 
    195162            ztrdu(:,:,:) = ua(:,:,:) 
    196163            ztrdv(:,:,:) = va(:,:,:) 
     
    198165            ztrdu(:,:,:) = ua(:,:,:) - ztrdu(:,:,:) 
    199166            ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:) 
    200             IF( l_trddyn ) THEN 
    201                 CALL trd_dyn( ztrdu, ztrdv, jpdyn_rvo, kt ) 
    202             END IF 
    203             IF( id_vrt_dia_rvo == "rvo" ) THEN 
    204                 CALL dyn_vrt_dia_3d( ztrdu, ztrdv, id_vrt_dia_rvo ) 
    205             END IF 
     167            CALL trd_dyn( ztrdu, ztrdv, jpdyn_rvo, kt ) 
    206168            ztrdu(:,:,:) = ua(:,:,:) 
    207169            ztrdv(:,:,:) = va(:,:,:) 
     
    209171            ztrdu(:,:,:) = ua(:,:,:) - ztrdu(:,:,:) 
    210172            ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:) 
    211             IF( l_trddyn ) THEN 
    212                 CALL trd_dyn( ztrdu, ztrdv, jpdyn_pvo, kt ) 
    213             END IF 
    214             IF( id_vrt_dia_pvo == "pvo" ) THEN 
    215                 CALL dyn_vrt_dia_3d( ztrdu, ztrdv, id_vrt_dia_pvo ) 
    216             END IF 
     173            CALL trd_dyn( ztrdu, ztrdv, jpdyn_pvo, kt ) 
    217174         ELSE 
    218175            CALL vor_een( kt, ntot, ua, va )                ! total vorticity 
     
    225182         &                     tab3d_2=va, clinfo2=       ' Va: ', mask2=vmask, clinfo3='dyn' ) 
    226183      ! 
    227       IF( l_trddyn .or. (id_vrt_dia_rvo == "rvo" ) )   CALL wrk_dealloc( jpi,jpj,jpk, ztrdu, ztrdv ) 
     184      IF( l_trddyn )   CALL wrk_dealloc( jpi,jpj,jpk, ztrdu, ztrdv ) 
    228185      ! 
    229186      IF( nn_timing == 1 )  CALL timing_stop('dyn_vor') 
  • branches/UKMO/GO6_dyn_vrt_diag/NEMOGCM/NEMO/OPA_SRC/DYN/dynzad.F90

    r8300 r9022  
    1818   USE trd_oce        ! trends: ocean variables 
    1919   USE trddyn         ! trend manager: dynamics 
    20    USE divcur         ! for dyn_vrt_dia_3d 
    2120   ! 
    2221   USE in_out_manager ! I/O manager 
     
    6059      INTEGER, INTENT(in) ::   kt   ! ocean time-step inedx 
    6160      ! 
    62       CHARACTER(len=3) :: id_vrt_dia_zad = "zad" ! TODO remove once flags set properly 
    63       INTEGER  ::   ji, jj, jk                   ! dummy loop indices 
    64       INTEGER  ::   zua, zva                     ! temporary scalars 
     61      INTEGER  ::   ji, jj, jk      ! dummy loop indices 
     62      REAL(wp) ::   zua, zva        ! temporary scalars 
    6563      REAL(wp), POINTER, DIMENSION(:,:,:) ::  zwuw , zwvw 
    6664      REAL(wp), POINTER, DIMENSION(:,:  ) ::  zww 
     
    7876      ENDIF 
    7977 
    80       IF( l_trddyn .OR. ( id_vrt_dia_zad == "zad" ) )   THEN         ! Save ua and va trends 
     78      IF( l_trddyn )   THEN         ! Save ua and va trends 
    8179         CALL wrk_alloc( jpi, jpj, jpk, ztrdu, ztrdv )  
    8280         ztrdu(:,:,:) = ua(:,:,:)  
     
    132130      END DO 
    133131 
    134       IF( l_trddyn .OR. ( id_vrt_dia_zad == "zad" ) ) THEN      ! save the vertical advection trends for diagnostic 
     132      IF( l_trddyn ) THEN           ! save the vertical advection trends for diagnostic 
    135133         ztrdu(:,:,:) = ua(:,:,:) - ztrdu(:,:,:) 
    136134         ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:) 
    137          ! 
    138          IF( id_vrt_dia_zad == "zad" ) THEN 
    139              CALL dyn_vrt_dia_3d( ztrdu, ztrdv, id_vrt_dia_zad ) 
    140          ENDIF 
    141          ! 
    142          IF( l_trddyn ) THEN 
    143              CALL trd_dyn( ztrdu, ztrdv, jpdyn_zad, kt ) 
    144          ENDIF 
    145          ! 
     135         CALL trd_dyn( ztrdu, ztrdv, jpdyn_zad, kt ) 
    146136         CALL wrk_dealloc( jpi, jpj, jpk, ztrdu, ztrdv )  
    147137      ENDIF 
  • branches/UKMO/GO6_dyn_vrt_diag/NEMOGCM/NEMO/OPA_SRC/DYN/dynzdf.F90

    r8300 r9022  
    2727   USE wrk_nemo        ! Memory Allocation 
    2828   USE timing          ! Timing 
    29    USE divcur          ! for dyn_vrt_dia 
    3029 
    3130   IMPLICIT NONE 
     
    6059      ! 
    6160      REAL(wp), POINTER, DIMENSION(:,:,:) ::  ztrdu, ztrdv 
    62       CHARACTER(len=3) :: id_vrt_dia_zdf = "zdf"   ! TODO remove once flags properly set 
    6361      !!--------------------------------------------------------------------- 
    6462      ! 
     
    7068      ENDIF 
    7169 
    72       IF( l_trddyn .or. ( id_vrt_dia_zdf == "zdf" ) )   THEN      ! temporary save of ta and sa trends 
     70      IF( l_trddyn )   THEN                      ! temporary save of ta and sa trends 
    7371         CALL wrk_alloc( jpi, jpj, jpk, ztrdu, ztrdv )  
    7472         ztrdu(:,:,:) = ua(:,:,:) 
     
    9088      END SELECT 
    9189 
    92       IF( l_trddyn .or. ( id_vrt_dia_zdf == "zdf" ) )   THEN      ! save the vertical diffusive trends for further diagnostics 
     90      IF( l_trddyn )   THEN                      ! save the vertical diffusive trends for further diagnostics 
    9391         ztrdu(:,:,:) = ua(:,:,:) - ztrdu(:,:,:) 
    9492         ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:) 
    95          ! 
    96          IF( id_vrt_dia_zdf == "zdf" ) THEN 
    97              CALL dyn_vrt_dia_3d( ztrdu, ztrdv, id_vrt_dia_zdf ) 
    98          END IF 
    99          ! 
    100          IF( l_trddyn ) THEN 
    101              CALL trd_dyn( ztrdu, ztrdv, jpdyn_zdf, kt ) 
    102          END IF 
    103          ! 
     93         CALL trd_dyn( ztrdu, ztrdv, jpdyn_zdf, kt ) 
    10494         CALL wrk_dealloc( jpi, jpj, jpk, ztrdu, ztrdv )  
    10595      ENDIF 
  • branches/UKMO/GO6_dyn_vrt_diag/NEMOGCM/NEMO/OPA_SRC/TRD/trdini.F90

    r6486 r9022  
    8888 
    8989!!gm end 
    90       IF( ln_tra_mxl .OR. ln_vor_trd )   CALL ctl_stop( 'ML tracer and Barotropic vorticity diags are still using old IOIPSL' ) 
     90      IF( ln_tra_mxl )   CALL ctl_stop( 'ML tracer and Barotropic vorticity diags are still using old IOIPSL' ) 
    9191!!gm end 
    9292      ! 
    93       IF( lk_vvl .AND. ( l_trdtra .OR. l_trddyn ) )  CALL ctl_stop( 'trend diagnostics with variable volume not validated' ) 
     93      IF( lk_vvl .AND. l_trdtra )  CALL ctl_stop( 'trend diagnostics with variable volume not validated' ) 
    9494       
    9595!!gm  : Potential BUG : 3D output only for vector invariant form!  add a ctl_stop or code the flux form case 
  • branches/UKMO/GO6_dyn_vrt_diag/NEMOGCM/NEMO/OPA_SRC/TRD/trdvor.F90

    r6486 r9022  
    2828   USE lib_mpp         ! MPP library 
    2929   USE wrk_nemo        ! Memory allocation 
     30   USE iom             ! I/O manager library 
    3031 
    3132   IMPLICIT NONE 
     
    4041   PUBLIC   trd_vor_alloc  ! routine called by nemogcm.F90 
    4142 
    42    INTEGER ::   nh_t, nmoydpvor, nidvor, nhoridvor, ndimvor1, icount   ! needs for IOIPSL output 
    43    INTEGER, SAVE, ALLOCATABLE, DIMENSION(:) ::   ndexvor1   ! needed for IOIPSL output 
     43   INTEGER ::   nmoydpvor  ! needs for IOIPSL output 
    4444   INTEGER ::   ndebug     ! (0/1) set it to 1 in case of problem to have more print 
    4545 
    46    REAL(wp), SAVE, ALLOCATABLE, DIMENSION(:,:)   ::   vor_avr      ! average 
    47    REAL(wp), SAVE, ALLOCATABLE, DIMENSION(:,:)   ::   vor_avrb     ! before vorticity (kt-1) 
    48    REAL(wp), SAVE, ALLOCATABLE, DIMENSION(:,:)   ::   vor_avrbb    ! vorticity at begining of the nwrite-1 timestep averaging period 
    49    REAL(wp), SAVE, ALLOCATABLE, DIMENSION(:,:)   ::   vor_avrbn    ! after vorticity at time step after the 
    50    REAL(wp), SAVE, ALLOCATABLE, DIMENSION(:,:)   ::   rotot        ! begining of the NWRITE-1 timesteps 
    51    REAL(wp), SAVE, ALLOCATABLE, DIMENSION(:,:)   ::   vor_avrtot   ! 
    52    REAL(wp), SAVE, ALLOCATABLE, DIMENSION(:,:)   ::   vor_avrres   ! 
    53    REAL(wp), SAVE, ALLOCATABLE, DIMENSION(:,:,:) ::   vortrd       ! curl of trends 
     46   REAL(wp), SAVE, ALLOCATABLE, DIMENSION(:,:,:)    ::   vor        ! time average 
     47   REAL(wp), SAVE, ALLOCATABLE, DIMENSION(:,:,:)    ::   vor_b      ! before vorticity (kt-1) 
     48   REAL(wp), SAVE, ALLOCATABLE, DIMENSION(:,:,:)    ::   vor_bb     ! vorticity at begining of the nwrite-1 timestep averaging period 
     49   REAL(wp), SAVE, ALLOCATABLE, DIMENSION(:,:,:)    ::   vor_bn     ! after vorticity at time step after 
     50   REAL(wp), SAVE, ALLOCATABLE, DIMENSION(:,:,:)    ::   rotot      ! begining of the NWRITE-1 timesteps 
     51   REAL(wp), SAVE, ALLOCATABLE, DIMENSION(:,:)      ::   intvor_tot ! 
     52   REAL(wp), SAVE, ALLOCATABLE, DIMENSION(:,:)      ::   avrvor_tot ! 
     53   REAL(wp), SAVE, ALLOCATABLE, DIMENSION(:,:)      ::   intvor_res ! 
     54   REAL(wp), SAVE, ALLOCATABLE, DIMENSION(:,:)      ::   avrvor_res ! 
     55   REAL(wp), SAVE, ALLOCATABLE, DIMENSION(:,:,:)    ::   intvortrd  ! vertically integrated curl of trends 
     56   REAL(wp), SAVE, ALLOCATABLE, DIMENSION(:,:,:)    ::   avrvortrd  ! vertically averaged curl of trends 
    5457          
    5558   CHARACTER(len=12) ::   cvort 
     
    7073      !!                  ***  ROUTINE trd_vor_alloc  *** 
    7174      !!---------------------------------------------------------------------------- 
    72       ALLOCATE( vor_avr   (jpi,jpj) , vor_avrb(jpi,jpj) , vor_avrbb (jpi,jpj) ,   & 
    73          &      vor_avrbn (jpi,jpj) , rotot   (jpi,jpj) , vor_avrtot(jpi,jpj) ,   & 
    74          &      vor_avrres(jpi,jpj) , vortrd  (jpi,jpj,jpltot_vor) ,              & 
    75          &      ndexvor1  (jpi*jpj)                                ,   STAT= trd_vor_alloc ) 
     75      ALLOCATE( vor   (jpi,jpj,jpvor_types) , vor_b  (jpi,jpj,jpvor_types) ,   & 
     76         &      vor_bb(jpi,jpj,jpvor_types) , vor_bn (jpi,jpj,jpvor_types) ,   & 
     77         &      rotot (jpi,jpj,jpvor_types) , intvor_tot(jpi,jpj)          ,   & 
     78         &      avrvor_tot(jpi,jpj)         , intvor_res(jpi,jpj)          ,   & 
     79         &      avrvor_res(jpi,jpj)         , intvortrd(jpi,jpj,jpltot_vor),   & 
     80         &      avrvortrd(jpi,jpj,jpltot_vor),              STAT= trd_vor_alloc ) 
    7681         ! 
    7782      IF( lk_mpp             )   CALL mpp_sum ( trd_vor_alloc ) 
     
    122127         CALL trd_vor_iom( kt ) 
    123128      END SELECT 
    124       ! 
    125129      CALL wrk_dealloc( jpi, jpj, ztswu, ztswv ) 
    126130      ! 
     
    138142      !! 
    139143      !! ** Action :   trends : 
    140       !!                  vortrd (,, 1) = Pressure Gradient Trend 
    141       !!                  vortrd (,, 2) = KE Gradient Trend 
    142       !!                  vortrd (,, 3) = Relative Vorticity Trend 
    143       !!                  vortrd (,, 4) = Coriolis Term Trend 
    144       !!                  vortrd (,, 5) = Horizontal Diffusion Trend 
    145       !!                  vortrd (,, 6) = Vertical Advection Trend 
    146       !!                  vortrd (,, 7) = Vertical Diffusion Trend 
    147       !!                  vortrd (,, 8) = Surface Pressure Grad. Trend 
    148       !!                  vortrd (,, 9) = Beta V 
    149       !!                  vortrd (,,10) = forcing term 
    150       !!                  vortrd (,,11) = bottom friction term 
    151       !!                  rotot(,) : total cumulative trends over nwrite-1 time steps 
    152       !!                  vor_avrtot(,) : first membre of vrticity equation 
    153       !!                  vor_avrres(,) : residual = dh/dt entrainment 
     144      !!                  intvortrd (,, 1) = Vertical integral of Pressure Gradient Trend 
     145      !!                  intvortrd (,, 2) = Vertical integral of KE Gradient Trend 
     146      !!                  intvortrd (,, 3) = Vertical integral of Relative Vorticity Trend 
     147      !!                  intvortrd (,, 4) = Vertical integral of Coriolis Term Trend 
     148      !!                  intvortrd (,, 5) = Vertical integral of Horizontal Diffusion Trend 
     149      !!                  intvortrd (,, 6) = Vertical integral of Vertical Advection Trend 
     150      !!                  intvortrd (,, 7) = Vertical integral of Vertical Diffusion Trend 
     151      !!                  intvortrd (,, 8) = Vertical integral of Surface Pressure Grad. Trend 
     152      !!                  intvortrd (,,10) = Vertical integral of forcing term 
     153      !!      intvortrd (,,11) = Vertical integral of bottom friction term 
     154      !!                  avrvortrd (,, 1) = Vertical average of Pressure Gradient Trend 
     155      !!                  avrvortrd (,, 2) = Vertical average of KE Gradient Trend 
     156      !!                  avrvortrd (,, 3) = Vertical average of Relative Vorticity Trend 
     157      !!                  avrvortrd (,, 4) = Vertical average of Coriolis Term Trend 
     158      !!                  avrvortrd (,, 5) = Vertical average of Horizontal Diffusion Trend 
     159      !!                  avrvortrd (,, 6) = Vertical average of Vertical Advection Trend 
     160      !!                  avrvortrd (,, 7) = Vertical average of Vertical Diffusion Trend 
     161      !!                  avrvortrd (,, 8) = Vertical average of Surface Pressure Grad. Trend 
     162      !!                  avrvortrd (,,10) = Vertical average of forcing term 
     163      !!      avrvortrd (,,11) = Vertical average of bottom friction term 
     164      !!                  rotot  (,,1) = total cumulative vertical integral trends over nwrite-1 time steps 
     165      !!                  rotot  (,,2) = total cumulative vertical average trends over nwrite-1 time steps 
     166      !!                  intvor_tot(,,) = first membre of vrticity equation for vertical integral 
     167      !!                  avrvor_tot(,,) = first membre of vrticity equation for vertical average 
     168      !!                  intvor_res(,,) = residual = dh/dt entrainment for vertical integral 
     169      !!                  avrvor_res(,,) = residual = dh/dt entrainment for vertical average 
    154170      !! 
    155171      !!      trends output in netCDF format using ioipsl 
     
    194210      END SELECT 
    195211 
    196       ! Average except for Beta.V 
     212      ! Curl of the vertical integral 
     213      DO ji = 1, jpim1 
     214         DO jj = 1, jpjm1 
     215            intvortrd(ji,jj,ktrd) = (    zvdpvor(ji+1,jj) - zvdpvor(ji,jj)       & 
     216                    &                - ( zudpvor(ji,jj+1) - zudpvor(ji,jj) )   ) / ( e1f(ji,jj) * e2f(ji,jj) ) 
     217         END DO 
     218      END DO 
     219      intvortrd(:,:,ktrd) = intvortrd(:,:,ktrd) * fmask(:,:,1)      ! Surface mask 
     220 
    197221      zudpvor(:,:) = zudpvor(:,:) * hur(:,:) 
    198222      zvdpvor(:,:) = zvdpvor(:,:) * hvr(:,:) 
    199223    
    200       ! Curl 
     224      ! Curl of the vertical average 
    201225      DO ji = 1, jpim1 
    202226         DO jj = 1, jpjm1 
    203             vortrd(ji,jj,ktrd) = (    zvdpvor(ji+1,jj) - zvdpvor(ji,jj)       & 
     227            avrvortrd(ji,jj,ktrd) = (    zvdpvor(ji+1,jj) - zvdpvor(ji,jj)       & 
    204228                 &                - ( zudpvor(ji,jj+1) - zudpvor(ji,jj) )   ) / ( e1f(ji,jj) * e2f(ji,jj) ) 
    205229         END DO 
    206230      END DO 
    207       vortrd(:,:,ktrd) = vortrd(:,:,ktrd) * fmask(:,:,1)      ! Surface mask 
     231      avrvortrd(:,:,ktrd) = avrvortrd(:,:,ktrd) * fmask(:,:,1)      ! Surface mask 
    208232 
    209233      IF( ndebug /= 0 ) THEN 
     
    227251      !! 
    228252      !! ** Action :     trends : 
    229       !!                  vortrd (,,1) = Pressure Gradient Trend 
    230       !!                  vortrd (,,2) = KE Gradient Trend 
    231       !!                  vortrd (,,3) = Relative Vorticity Trend 
    232       !!                  vortrd (,,4) = Coriolis Term Trend 
    233       !!                  vortrd (,,5) = Horizontal Diffusion Trend 
    234       !!                  vortrd (,,6) = Vertical Advection Trend 
    235       !!                  vortrd (,,7) = Vertical Diffusion Trend 
    236       !!                  vortrd (,,8) = Surface Pressure Grad. Trend 
    237       !!                  vortrd (,,9) = Beta V 
    238       !!                  vortrd (,,10) = forcing term 
    239       !!      vortrd (,,11) = bottom friction term 
    240       !!                  rotot(,) : total cumulative trends over nwrite-1 time steps 
    241       !!                  vor_avrtot(,) : first membre of vrticity equation 
    242       !!                  vor_avrres(,) : residual = dh/dt entrainment 
    243       !! 
    244       !!      trends output in netCDF format using ioipsl 
     253      !!                  intvortrd (,, 1) = Vertical integral of Pressure Gradient Trend 
     254      !!                  intvortrd (,, 2) = Vertical integral of KE Gradient Trend 
     255      !!                  intvortrd (,, 3) = Vertical integral of Relative Vorticity Trend 
     256      !!                  intvortrd (,, 4) = Vertical integral of Coriolis Term Trend 
     257      !!                  intvortrd (,, 5) = Vertical integral of Horizontal Diffusion Trend 
     258      !!                  intvortrd (,, 6) = Vertical integral of Vertical Advection Trend 
     259      !!                  intvortrd (,, 7) = Vertical integral of Vertical Diffusion Trend 
     260      !!                  intvortrd (,, 8) = Vertical integral of Surface Pressure Grad. Trend 
     261      !!                  intvortrd (,,10) = Vertical integral of forcing term 
     262      !!      intvortrd (,,11) = Vertical integral of bottom friction term 
     263      !!                  avrvortrd (,, 1) = Vertical average of Pressure Gradient Trend 
     264      !!                  avrvortrd (,, 2) = Vertical average of KE Gradient Trend 
     265      !!                  avrvortrd (,, 3) = Vertical average of Relative Vorticity Trend 
     266      !!                  avrvortrd (,, 4) = Vertical average of Coriolis Term Trend 
     267      !!                  avrvortrd (,, 5) = Vertical average of Horizontal Diffusion Trend 
     268      !!                  avrvortrd (,, 6) = Vertical average of Vertical Advection Trend 
     269      !!                  avrvortrd (,, 7) = Vertical average of Vertical Diffusion Trend 
     270      !!                  avrvortrd (,, 8) = Vertical average of Surface Pressure Grad. Trend 
     271      !!                  avrvortrd (,,10) = Vertical average of forcing term 
     272      !!      avrvortrd (,,11) = Vertical average of bottom friction term 
     273      !!                  rotot  (,,1) = total cumulative vertical integral trends over nwrite-1 time steps 
     274      !!                  rotot  (,,2) = total cumulative vertical average trends over nwrite-1 time steps 
     275      !!                  intvor_tot(,,) = first membre of vrticity equation for vertical integral 
     276      !!                  avrvor_tot(,,) = first membre of vrticity equation for vertical average 
     277      !!                  intvor_res(,,) = residual = dh/dt entrainment for vertical integral 
     278      !!                  avrvor_res(,,) = residual = dh/dt entrainment for vertical average 
     279      !! 
     280      !!      trends output in netCDF format using iom 
    245281      !!---------------------------------------------------------------------- 
    246282      ! 
     
    250286      ! 
    251287      INTEGER ::   ji, jj, jk   ! dummy loop indices 
    252       REAL(wp), POINTER, DIMENSION(:,:) :: zubet  , zvbet    ! Beta.V    
    253288      REAL(wp), POINTER, DIMENSION(:,:) :: zudpvor, zvdpvor  ! total cmulative trends 
    254289      !!---------------------------------------------------------------------- 
    255290      
    256       CALL wrk_alloc( jpi,jpj, zubet, zvbet, zudpvor, zvdpvor )                                    
     291      CALL wrk_alloc( jpi,jpj, zudpvor, zvdpvor )                                    
    257292 
    258293      ! Initialization 
    259       zubet  (:,:) = 0._wp 
    260       zvbet  (:,:) = 0._wp 
    261294      zudpvor(:,:) = 0._wp 
    262295      zvdpvor(:,:) = 0._wp 
     
    273306        zvdpvor(:,:) = zvdpvor(:,:) + pvtrdvor(:,:,jk) * fse3v(:,:,jk) * e2v(:,:) * vmask(:,:,jk) 
    274307      END DO 
    275  
    276       ! Planetary vorticity: 2nd computation (Beta.V term) store the vertical sum 
    277       ! as Beta.V term need intergration, not average 
    278       IF( ktrd == jpvor_pvo ) THEN  
    279          zubet(:,:) = zudpvor(:,:) 
    280          zvbet(:,:) = zvdpvor(:,:) 
    281          DO ji = 1, jpim1 
    282             DO jj = 1, jpjm1 
    283                vortrd(ji,jj,jpvor_bev) = (    zvbet(ji+1,jj) - zvbet(ji,jj)     & 
    284                   &                       - ( zubet(ji,jj+1) - zubet(ji,jj) ) ) / ( e1f(ji,jj) * e2f(ji,jj) ) 
    285             END DO 
    286          END DO 
    287          ! Average of the Curl and Surface mask 
    288          vortrd(:,:,jpvor_bev) = vortrd(:,:,jpvor_bev) * hur(:,:) * fmask(:,:,1) 
    289       ENDIF 
     308      ! 
     309      ! Curl of the integral 
     310      DO ji=1,jpim1 
     311         DO jj=1,jpjm1 
     312            intvortrd(ji,jj,ktrd) = (    zvdpvor(ji+1,jj) - zvdpvor(ji,jj)     & 
     313                  &                  - ( zudpvor(ji,jj+1) - zudpvor(ji,jj) ) ) / ( e1f(ji,jj) * e2f(ji,jj) ) 
     314         END DO 
     315      END DO 
     316      ! Surface mask 
     317      intvortrd(:,:,ktrd) = intvortrd(:,:,ktrd) * fmask(:,:,1) 
    290318      ! 
    291319      ! Average  
     
    293321      zvdpvor(:,:) = zvdpvor(:,:) * hvr(:,:) 
    294322      ! 
    295       ! Curl 
     323      ! Curl of the average 
    296324      DO ji=1,jpim1 
    297325         DO jj=1,jpjm1 
    298             vortrd(ji,jj,ktrd) = (    zvdpvor(ji+1,jj) - zvdpvor(ji,jj)     & 
    299                &                  - ( zudpvor(ji,jj+1) - zudpvor(ji,jj) ) ) / ( e1f(ji,jj) * e2f(ji,jj) ) 
     326            avrvortrd(ji,jj,ktrd) = (    zvdpvor(ji+1,jj) - zvdpvor(ji,jj)     & 
     327                 &                  - ( zudpvor(ji,jj+1) - zudpvor(ji,jj) ) ) / ( e1f(ji,jj) * e2f(ji,jj) ) 
    300328         END DO 
    301329      END DO 
    302330      ! Surface mask 
    303       vortrd(:,:,ktrd) = vortrd(:,:,ktrd) * fmask(:,:,1) 
     331      avrvortrd(:,:,ktrd) = avrvortrd(:,:,ktrd) * fmask(:,:,1) 
    304332    
    305333      IF( ndebug /= 0 ) THEN 
     
    308336      ENDIF 
    309337      ! 
    310       CALL wrk_dealloc( jpi,jpj, zubet, zvbet, zudpvor, zvdpvor )                                    
     338      CALL wrk_dealloc( jpi,jpj, zudpvor, zvdpvor )                                    
    311339      ! 
    312340   END SUBROUTINE trd_vor_zint_3d 
     
    338366      ! --------------------------------------------------- 
    339367 
    340       IF( kt > nit000 )   vor_avrb(:,:) = vor_avr(:,:) 
     368      IF( kt > nit000 ) THEN 
     369         vor_b(:,:,:) = vor(:,:,:) 
     370      ENDIF 
    341371 
    342372      ! I.2 vertically integrated vorticity 
    343373      !  ---------------------------------- 
    344374 
    345       vor_avr   (:,:) = 0._wp 
     375      vor     (:,:,:) = 0._wp 
    346376      zun       (:,:) = 0._wp 
    347377      zvn       (:,:) = 0._wp 
    348       vor_avrtot(:,:) = 0._wp 
    349       vor_avrres(:,:) = 0._wp 
     378      intvor_tot(:,:) = 0._wp 
     379      avrvor_tot(:,:) = 0._wp 
     380      intvor_res(:,:) = 0._wp 
     381      avrvor_res(:,:) = 0._wp 
    350382       
    351       ! Vertically averaged velocity 
     383      ! Vertically integrated velocity 
    352384      DO jk = 1, jpk - 1 
    353385         zun(:,:) = zun(:,:) + e1u(:,:) * un(:,:,jk) * fse3u(:,:,jk) 
     
    355387      END DO 
    356388  
     389      ! Curl of the vertical integral 
     390      DO ji = 1, jpim1 
     391         DO jj = 1, jpjm1 
     392            vor(ji,jj,jpvor_int) = (  ( zvn(ji+1,jj) - zvn(ji,jj) )    & 
     393               &              - ( zun(ji,jj+1) - zun(ji,jj) ) ) / ( e1f(ji,jj) * e2f(ji,jj) ) * fmask(ji,jj,1) 
     394         END DO 
     395      END DO 
     396 
     397      ! Vertically averaged velocity 
    357398      zun(:,:) = zun(:,:) * hur(:,:) 
    358399      zvn(:,:) = zvn(:,:) * hvr(:,:) 
    359400 
    360       ! Curl 
     401      ! Curl of the vertical average 
    361402      DO ji = 1, jpim1 
    362403         DO jj = 1, jpjm1 
    363             vor_avr(ji,jj) = (  ( zvn(ji+1,jj) - zvn(ji,jj) )    & 
     404            vor(ji,jj,jpvor_avr) = (  ( zvn(ji+1,jj) - zvn(ji,jj) )    & 
    364405               &              - ( zun(ji,jj+1) - zun(ji,jj) ) ) / ( e1f(ji,jj) * e2f(ji,jj) ) * fmask(ji,jj,1) 
    365406         END DO 
     
    373414      ! ------------------------------------------------------ 
    374415      IF( kt == nit000+1 ) THEN 
    375          vor_avrbb(:,:) = vor_avrb(:,:) 
    376          vor_avrbn(:,:) = vor_avr (:,:) 
     416         vor_bb(:,:,:) = vor_b(:,:,:) 
     417         vor_bn(:,:,:) = vor  (:,:,:) 
    377418      ENDIF 
    378419 
     
    384425         nmoydpvor = nmoydpvor + 1 
    385426         DO jl = 1, jpltot_vor 
    386             IF( jl /= 9 ) THEN 
    387                rotot(:,:) = rotot(:,:) + vortrd(:,:,jl) 
    388             ENDIF 
     427            rotot(:,:,jpvor_int) = rotot(:,:,jpvor_int) + intvortrd(:,:,jl) 
     428            rotot(:,:,jpvor_avr) = rotot(:,:,jpvor_avr) + avrvortrd(:,:,jl) 
    389429         END DO 
    390430      ENDIF 
     
    398438      itmod = kt - nit000 + 1 
    399439 
    400       IF( MOD( it, nn_trd ) == 0 ) THEN 
     440      IF( MOD( it, nwrite ) == 0 ) THEN 
    401441 
    402442         ! III.1 compute total trend 
    403443         ! ------------------------ 
    404444         zmean = 1._wp / (  REAL( nmoydpvor, wp ) * 2._wp * rdt  ) 
    405          vor_avrtot(:,:) = (  vor_avr(:,:) - vor_avrbn(:,:) + vor_avrb(:,:) - vor_avrbb(:,:) ) * zmean 
     445         intvor_tot(:,:) = (  vor(:,:,jpvor_int) - vor_bn(:,:,jpvor_int)      & 
     446                    &        + vor_b(:,:,jpvor_int) - vor_bb(:,:,jpvor_int) ) * zmean 
     447         avrvor_tot(:,:) = (  vor(:,:,jpvor_avr) - vor_bn(:,:,jpvor_avr)      & 
     448                    &        + vor_b(:,:,jpvor_avr) - vor_bb(:,:,jpvor_avr) ) * zmean 
    406449 
    407450 
     
    409452         ! --------------------- 
    410453         zmean = 1._wp / REAL( nmoydpvor, wp ) 
    411          vor_avrres(:,:) = vor_avrtot(:,:) - rotot(:,:) / zmean 
     454         intvor_res(:,:) = intvor_tot(:,:) - rotot(:,:,jpvor_int) * zmean 
     455         avrvor_res(:,:) = avrvor_tot(:,:) - rotot(:,:,jpvor_avr) * zmean 
    412456 
    413457         ! Boundary conditions 
    414          CALL lbc_lnk( vor_avrtot, 'F', 1. ) 
    415          CALL lbc_lnk( vor_avrres, 'F', 1. ) 
     458         CALL lbc_lnk( intvor_tot, 'F', 1. ) 
     459         CALL lbc_lnk( avrvor_tot, 'F', 1. ) 
     460         CALL lbc_lnk( intvor_res, 'F', 1. ) 
     461         CALL lbc_lnk( avrvor_res, 'F', 1. ) 
    416462 
    417463 
    418464         ! III.3 time evolution array swap 
    419465         ! ------------------------------ 
    420          vor_avrbb(:,:) = vor_avrb(:,:) 
    421          vor_avrbn(:,:) = vor_avr (:,:) 
     466         vor_bb(:,:,:) = vor_b(:,:,:) 
     467         vor_bn(:,:,:) = vor  (:,:,:) 
    422468         ! 
    423469         nmoydpvor = 0 
     
    430476      IF( kt >=  nit000+1 ) THEN 
    431477 
    432          IF( lwp .AND. MOD( itmod, nn_trd ) == 0 ) THEN 
     478         IF( lwp .AND. MOD( itmod, nwrite ) == 0 ) THEN 
    433479            WRITE(numout,*) '' 
    434480            WRITE(numout,*) 'trd_vor : write trends in the NetCDF file at kt = ', kt 
    435481            WRITE(numout,*) '~~~~~~~  ' 
    436482         ENDIF 
    437   
    438          CALL histwrite( nidvor,"sovortPh",it,vortrd(:,:,jpvor_prg),ndimvor1,ndexvor1)  ! grad Ph 
    439          CALL histwrite( nidvor,"sovortEk",it,vortrd(:,:,jpvor_keg),ndimvor1,ndexvor1)  ! Energy 
    440          CALL histwrite( nidvor,"sovozeta",it,vortrd(:,:,jpvor_rvo),ndimvor1,ndexvor1)  ! rel vorticity 
    441          CALL histwrite( nidvor,"sovortif",it,vortrd(:,:,jpvor_pvo),ndimvor1,ndexvor1)  ! coriolis 
    442          CALL histwrite( nidvor,"sovodifl",it,vortrd(:,:,jpvor_ldf),ndimvor1,ndexvor1)  ! lat diff 
    443          CALL histwrite( nidvor,"sovoadvv",it,vortrd(:,:,jpvor_zad),ndimvor1,ndexvor1)  ! vert adv 
    444          CALL histwrite( nidvor,"sovodifv",it,vortrd(:,:,jpvor_zdf),ndimvor1,ndexvor1)  ! vert diff 
    445          CALL histwrite( nidvor,"sovortPs",it,vortrd(:,:,jpvor_spg),ndimvor1,ndexvor1)  ! grad Ps 
    446          CALL histwrite( nidvor,"sovortbv",it,vortrd(:,:,jpvor_bev),ndimvor1,ndexvor1)  ! beta.V 
    447          CALL histwrite( nidvor,"sovowind",it,vortrd(:,:,jpvor_swf),ndimvor1,ndexvor1) ! wind stress 
    448          CALL histwrite( nidvor,"sovobfri",it,vortrd(:,:,jpvor_bfr),ndimvor1,ndexvor1) ! bottom friction 
    449          CALL histwrite( nidvor,"1st_mbre",it,vor_avrtot    ,ndimvor1,ndexvor1) ! First membre 
    450          CALL histwrite( nidvor,"sovorgap",it,vor_avrres    ,ndimvor1,ndexvor1) ! gap between 1st and 2 nd mbre 
    451          ! 
     483 
     484         ! Output the values for the vertical integral 
     485         CALL iom_put( "sovortPh_int", intvortrd (:,:,jpvor_prg) )  ! grad Ph 
     486         CALL iom_put( "sovortEk_int", intvortrd (:,:,jpvor_keg) )  ! Energy 
     487         CALL iom_put( "sovozeta_int", intvortrd (:,:,jpvor_rvo) )  ! rel vorticity 
     488         CALL iom_put( "sovortif_int", intvortrd (:,:,jpvor_pvo) )  ! coriolis 
     489         CALL iom_put( "sovodifl_int", intvortrd (:,:,jpvor_ldf) )  ! lat diff 
     490         CALL iom_put( "sovoadvv_int", intvortrd (:,:,jpvor_zad) )  ! vert adv 
     491         CALL iom_put( "sovodifv_int", intvortrd (:,:,jpvor_zdf) )  ! vert diff 
     492         CALL iom_put( "sovortPs_int", intvortrd (:,:,jpvor_spg) )  ! grad Ps 
     493         CALL iom_put( "sovowind_int", intvortrd (:,:,jpvor_swf) )  ! wind stress 
     494         CALL iom_put( "sovobfri_int", intvortrd (:,:,jpvor_bfr) )  ! bottom friction 
     495         CALL iom_put( "1st_mbre_int", intvor_tot(:,:)           )  ! First membre 
     496         CALL iom_put( "sovorgap_int", intvor_res(:,:)           )  ! gap between 1st and 2 nd mbre 
     497 
     498         ! Output the values for the vertical average 
     499         CALL iom_put( "sovortPh_avr", avrvortrd (:,:,jpvor_prg) )  ! grad Ph 
     500         CALL iom_put( "sovortEk_avr", avrvortrd (:,:,jpvor_keg) )  ! Energy 
     501         CALL iom_put( "sovozeta_avr", avrvortrd (:,:,jpvor_rvo) )  ! rel vorticity 
     502         CALL iom_put( "sovortif_avr", avrvortrd (:,:,jpvor_pvo) )  ! coriolis 
     503         CALL iom_put( "sovodifl_avr", avrvortrd (:,:,jpvor_ldf) )  ! lat diff 
     504         CALL iom_put( "sovoadvv_avr", avrvortrd (:,:,jpvor_zad) )  ! vert adv 
     505         CALL iom_put( "sovodifv_avr", avrvortrd (:,:,jpvor_zdf) )  ! vert diff 
     506         CALL iom_put( "sovortPs_avr", avrvortrd (:,:,jpvor_spg) )  ! grad Ps 
     507         CALL iom_put( "sovowind_avr", avrvortrd (:,:,jpvor_swf) )  ! wind stress 
     508         CALL iom_put( "sovobfri_avr", avrvortrd (:,:,jpvor_bfr) )  ! bottom friction 
     509         CALL iom_put( "1st_mbre_avr", avrvor_tot(:,:)           )  ! First membre 
     510         CALL iom_put( "sovorgap_avr", avrvor_res(:,:)           )  ! gap between 1st and 2 nd mbre 
     511 
    452512         IF( ndebug /= 0 ) THEN 
    453513            WRITE(numout,*) ' debuging trd_vor: III.4 done' 
     
    457517      ENDIF 
    458518      ! 
    459       IF( MOD( it, nn_trd ) == 0 ) rotot(:,:)=0 
    460       ! 
    461       IF( kt == nitend )   CALL histclo( nidvor ) 
     519      IF( MOD( it, nwrite ) == 0 ) THEN 
     520         rotot(:,:,:)=0 
     521      ENDIF 
    462522      ! 
    463523      CALL wrk_dealloc( jpi, jpj, zun, zvn )                                    
     
    473533      !!      from ocean surface down to control surface (NetCDF output) 
    474534      !!---------------------------------------------------------------------- 
    475       REAL(wp) ::   zjulian, zsto, zout 
    476       CHARACTER (len=40) ::   clhstnam 
    477       CHARACTER (len=40) ::   clop 
     535      !REAL(wp) ::   zjulian, zsto, zout 
     536      !CHARACTER (len=40) ::   clhstnam 
     537      !CHARACTER (len=40) ::   clop 
    478538      !!---------------------------------------------------------------------- 
    479539 
     
    504564      ! cumulated trends array init 
    505565      nmoydpvor = 0 
    506       rotot(:,:)=0 
    507       vor_avrtot(:,:)=0 
    508       vor_avrres(:,:)=0 
     566      rotot(:,:,:)=0 
     567      intvor_tot(:,:)=0 
     568      avrvor_tot(:,:)=0 
     569      intvor_res(:,:)=0 
     570      avrvor_res(:,:)=0 
    509571 
    510572      IF( ndebug /= 0 ) THEN 
     
    512574         CALL FLUSH(numout) 
    513575      ENDIF 
    514  
    515       !  ================================= 
    516       !   II. netCDF output initialization 
    517       !  ================================= 
    518  
    519       !----------------------------------------- 
    520       ! II.1 Define frequency of output and means 
    521       ! ----------------------------------------- 
    522       IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!) 
    523       ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time) 
    524       ENDIF 
    525 #if defined key_diainstant 
    526       zsto = nwrite*rdt 
    527       clop = "inst("//TRIM(clop)//")" 
    528 #else 
    529       zsto = rdt 
    530       clop = "ave("//TRIM(clop)//")" 
    531 #endif 
    532       zout = nn_trd*rdt 
    533  
    534       IF(lwp) WRITE(numout,*) '               netCDF initialization' 
    535  
    536       ! II.2 Compute julian date from starting date of the run 
    537       ! ------------------------ 
    538       CALL ymds2ju( nyear, nmonth, nday, rdt, zjulian ) 
    539       zjulian = zjulian - adatrj   !   set calendar origin to the beginning of the experiment 
    540       IF(lwp) WRITE(numout,*)' '   
    541       IF(lwp) WRITE(numout,*)'               Date 0 used :',nit000,    & 
    542          &                   ' YEAR ', nyear,' MONTH '      , nmonth,   & 
    543          &                   ' DAY ' , nday, 'Julian day : ', zjulian 
    544  
    545       ! II.3 Define the T grid trend file (nidvor) 
    546       ! --------------------------------- 
    547       CALL dia_nam( clhstnam, nn_trd, 'vort' )                  ! filename 
    548       IF(lwp) WRITE(numout,*) ' Name of NETCDF file ', clhstnam 
    549       CALL histbeg( clhstnam, jpi, glamf, jpj, gphif,1, jpi,   &  ! Horizontal grid : glamt and gphit 
    550          &          1, jpj, nit000-1, zjulian, rdt, nh_t, nidvor, domain_id=nidom, snc4chunks=snc4set ) 
    551       CALL wheneq( jpi*jpj, fmask, 1, 1., ndexvor1, ndimvor1 )    ! surface 
    552  
    553       ! Declare output fields as netCDF variables 
    554       CALL histdef( nidvor, "sovortPh", cvort//"grad Ph" , "s-2",        & ! grad Ph 
    555          &          jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 
    556       CALL histdef( nidvor, "sovortEk", cvort//"Energy", "s-2",          & ! Energy 
    557          &          jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 
    558       CALL histdef( nidvor, "sovozeta", cvort//"rel vorticity", "s-2",   & ! rel vorticity 
    559          &          jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 
    560       CALL histdef( nidvor, "sovortif", cvort//"coriolis", "s-2",        & ! coriolis 
    561          &          jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 
    562       CALL histdef( nidvor, "sovodifl", cvort//"lat diff ", "s-2",       & ! lat diff 
    563          &          jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 
    564       CALL histdef( nidvor, "sovoadvv", cvort//"vert adv", "s-2",        & ! vert adv 
    565          &          jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 
    566       CALL histdef( nidvor, "sovodifv", cvort//"vert diff" , "s-2",      & ! vert diff 
    567          &          jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 
    568       CALL histdef( nidvor, "sovortPs", cvort//"grad Ps", "s-2",         & ! grad Ps 
    569          &          jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 
    570       CALL histdef( nidvor, "sovortbv", cvort//"Beta V", "s-2",          & ! beta.V 
    571          &          jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 
    572       CALL histdef( nidvor, "sovowind", cvort//"wind stress", "s-2",     & ! wind stress 
    573          &          jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 
    574       CALL histdef( nidvor, "sovobfri", cvort//"bottom friction", "s-2", & ! bottom friction 
    575          &          jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 
    576       CALL histdef( nidvor, "1st_mbre", cvort//"1st mbre", "s-2",        & ! First membre 
    577          &          jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 
    578       CALL histdef( nidvor, "sovorgap", cvort//"gap", "s-2",             & ! gap between 1st and 2 nd mbre 
    579          &          jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 
    580       CALL histend( nidvor, snc4set ) 
    581  
    582       IF( ndebug /= 0 ) THEN 
    583          WRITE(numout,*) ' debuging trd_vor_init: II. done' 
    584          CALL FLUSH(numout) 
    585       ENDIF 
    586576      ! 
    587577   END SUBROUTINE trd_vor_init 
  • branches/UKMO/GO6_dyn_vrt_diag/NEMOGCM/NEMO/OPA_SRC/TRD/trdvor_oce.F90

    r6486 r9022  
    1313 
    1414   !                                               !!* vorticity trends index 
    15    INTEGER, PUBLIC, PARAMETER ::   jpltot_vor = 11  !: Number of vorticity trend terms 
     15   INTEGER, PUBLIC, PARAMETER ::   jpltot_vor = 10  !: Number of vorticity trend terms 
    1616   ! 
    1717   INTEGER, PUBLIC, PARAMETER ::   jpvor_prg =  1   !: Pressure Gradient Trend 
     
    2323   INTEGER, PUBLIC, PARAMETER ::   jpvor_zdf =  7   !: Vertical Diffusion Trend 
    2424   INTEGER, PUBLIC, PARAMETER ::   jpvor_spg =  8   !: Surface Pressure Grad. Trend 
    25    INTEGER, PUBLIC, PARAMETER ::   jpvor_bev =  9   !: Beta V 
    26    INTEGER, PUBLIC, PARAMETER ::   jpvor_swf = 10   !: wind stress forcing term 
    27    INTEGER, PUBLIC, PARAMETER ::   jpvor_bfr = 11   !: bottom friction term 
     25   INTEGER, PUBLIC, PARAMETER ::   jpvor_swf =  9   !: wind stress forcing term 
     26   INTEGER, PUBLIC, PARAMETER ::   jpvor_bfr = 10   !: bottom friction term 
     27   ! 
     28   INTEGER, PUBLIC, PARAMETER ::   jpvor_types = 2  !: bottom friction term 
     29   INTEGER, PUBLIC, PARAMETER ::   jpvor_int = 1    !: bottom friction term 
     30   INTEGER, PUBLIC, PARAMETER ::   jpvor_avr = 2    !: bottom friction term 
    2831 
    2932   !!---------------------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.