Changeset 7646 for trunk/NEMOGCM/TOOLS/SIREN/src/vgrid.f90
- Timestamp:
- 2017-02-06T10:25:03+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/TOOLS/SIREN/src/vgrid.f90
r6393 r7646 14 14 !> CALL vgrid_zgr_z(dd_gdepw(:), dd_gdept(:), dd_e3w(:), dd_e3t(:), 15 15 !> dd_ppkth, dd_ppkth2, dd_ppacr, dd_ppacr2, 16 !> dd_ppdzmin, dd_pphmax, dd_pp_to_be_computed,16 !> dd_ppdzmin, dd_pphmax, 17 17 !> dd_ppa0, dd_ppa1, dd_ppa2, dd_ppsur) 18 18 !> @endcode … … 26 26 !> - dd_ppdzmin see NEMO documentation 27 27 !> - dd_pphmax see NEMO documentation 28 !> - dd_pp_to_be_computed see NEMO documentation29 28 !> - dd_ppa1 see NEMO documentation 30 29 !> - dd_ppa2 see NEMO documentation … … 133 132 !> @param[in] dd_ppdzmin 134 133 !> @param[in] dd_pphmax 135 !> @param[in] dd_pp_to_be_computed136 134 !> @param[in] dd_ppa1 137 135 !> @param[in] dd_ppa2 … … 140 138 !------------------------------------------------------------------- 141 139 SUBROUTINE vgrid_zgr_z( dd_gdepw, dd_gdept, dd_e3w, dd_e3t, & 142 & dd_e3w_1d, dd_e3t_1d, &140 & dd_e3w_1d, dd_e3t_1d, & 143 141 & dd_ppkth, dd_ppkth2, dd_ppacr, dd_ppacr2, & 144 & dd_ppdzmin, dd_pphmax, dd_pp_to_be_computed,&142 & dd_ppdzmin, dd_pphmax, & 145 143 & dd_ppa0, dd_ppa1, dd_ppa2, dd_ppsur ) 146 144 IMPLICIT NONE … … 160 158 REAL(dp) , INTENT(IN ) :: dd_ppdzmin 161 159 REAL(dp) , INTENT(IN ) :: dd_pphmax 162 REAL(dp) , INTENT(IN ) :: dd_pp_to_be_computed160 REAL(dp), PARAMETER :: dp_pp_to_be_computed = NF90_FILL_DOUBLE 163 161 164 162 REAL(dp) , INTENT(IN ) :: dd_ppa0 … … 207 205 ! za0, za1, zsur are computed from ppdzmin , pphmax, ppkth, ppacr 208 206 ! 209 IF( dd_ppa1 == d d_pp_to_be_computed .AND. &210 & dd_ppa0 == d d_pp_to_be_computed .AND. &211 & dd_ppsur == d d_pp_to_be_computed ) THEN207 IF( dd_ppa1 == dp_pp_to_be_computed .AND. & 208 & dd_ppa0 == dp_pp_to_be_computed .AND. & 209 & dd_ppsur == dp_pp_to_be_computed ) THEN 212 210 dl_za1 = ( dl_zdzmin - dl_zhmax / REAL((il_jpk-1),dp) ) & 213 211 & / ( TANH((1-dl_zkth)/dl_zacr) - dl_zacr/REAL((il_jpk-1),dp) & … … 667 665 668 666 !namelist (intialise with GLORYS 75 levels parameters) 669 REAL(dp) :: dn_pp_to_be_computed = 0._dp670 667 REAL(dp) :: dn_ppsur = -3958.951371276829_dp 671 668 REAL(dp) :: dn_ppa0 = 103.9530096000000_dp … … 684 681 !---------------------------------------------------------------- 685 682 NAMELIST /namzgr/ & 686 & dn_pp_to_be_computed, &687 683 & dn_ppsur, & 688 684 & dn_ppa0, & … … 788 784 & dl_e3w_1d, dl_e3t_1d, & 789 785 & dn_ppkth, dn_ppkth2, dn_ppacr, dn_ppacr2, & 790 & dn_ppdzmin, dn_pphmax, dn_pp_to_be_computed,&786 & dn_ppdzmin, dn_pphmax, & 791 787 & dn_ppa0, dn_ppa1, dn_ppa2, dn_ppsur ) 792 788
Note: See TracChangeset
for help on using the changeset viewer.