- Timestamp:
- 2017-12-23T13:27:17+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/DOM/domvvl.F90
r9124 r9168 966 966 REWIND( numnam_ref ) ! Namelist nam_vvl in reference namelist : 967 967 READ ( numnam_ref, nam_vvl, IOSTAT = ios, ERR = 901) 968 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_vvl in reference namelist', lwp ) 969 ! 968 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_vvl in reference namelist', lwp ) 970 969 REWIND( numnam_cfg ) ! Namelist nam_vvl in configuration namelist : Parameters of the run 971 970 READ ( numnam_cfg, nam_vvl, IOSTAT = ios, ERR = 902 ) 972 902 IF( ios /=0 ) CALL ctl_nam ( ios , 'nam_vvl in configuration namelist', lwp )971 902 IF( ios > 0 ) CALL ctl_nam ( ios , 'nam_vvl in configuration namelist', lwp ) 973 972 IF(lwm) WRITE ( numond, nam_vvl ) 974 973 ! … … 977 976 WRITE(numout,*) 'dom_vvl_ctl : choice/control of the variable vertical coordinate' 978 977 WRITE(numout,*) '~~~~~~~~~~~' 979 WRITE(numout,*) ' 980 WRITE(numout,*) ' 981 WRITE(numout,*) ' 982 WRITE(numout,*) ' 983 WRITE(numout,*) ' 978 WRITE(numout,*) ' Namelist nam_vvl : chose a vertical coordinate' 979 WRITE(numout,*) ' zstar ln_vvl_zstar = ', ln_vvl_zstar 980 WRITE(numout,*) ' ztilde ln_vvl_ztilde = ', ln_vvl_ztilde 981 WRITE(numout,*) ' layer ln_vvl_layer = ', ln_vvl_layer 982 WRITE(numout,*) ' ztilde as zstar ln_vvl_ztilde_as_zstar = ', ln_vvl_ztilde_as_zstar 984 983 WRITE(numout,*) ' ztilde near the equator ln_vvl_zstar_at_eqtor = ', ln_vvl_zstar_at_eqtor 985 ! WRITE(numout,*) ' Namelist nam_vvl : chose kinetic-to-potential energy conservation' 986 ! WRITE(numout,*) ' ln_vvl_kepe = ', ln_vvl_kepe 987 WRITE(numout,*) ' Namelist nam_vvl : thickness diffusion coefficient' 988 WRITE(numout,*) ' rn_ahe3 = ', rn_ahe3 989 WRITE(numout,*) ' Namelist nam_vvl : maximum e3t deformation fractional change' 990 WRITE(numout,*) ' rn_zdef_max = ', rn_zdef_max 984 WRITE(numout,*) ' !' 985 WRITE(numout,*) ' thickness diffusion coefficient rn_ahe3 = ', rn_ahe3 986 WRITE(numout,*) ' maximum e3t deformation fractional change rn_zdef_max = ', rn_zdef_max 991 987 IF( ln_vvl_ztilde_as_zstar ) THEN 992 WRITE(numout,*) ' ztilde running in zstar emulation mode;'993 WRITE(numout,*) ' 994 WRITE(numout,*) ' 995 WRITE(numout,*) ' rn_rst_e3t = 0.0'996 WRITE(numout,*) ' 997 WRITE(numout,*) ' rn_lf_cutoff =1.0/rdt'988 WRITE(numout,*) ' ztilde running in zstar emulation mode (ln_vvl_ztilde_as_zstar=T) ' 989 WRITE(numout,*) ' ignoring namelist timescale parameters and using:' 990 WRITE(numout,*) ' hard-wired : z-tilde to zstar restoration timescale (days)' 991 WRITE(numout,*) ' rn_rst_e3t = 0.e0' 992 WRITE(numout,*) ' hard-wired : z-tilde cutoff frequency of low-pass filter (days)' 993 WRITE(numout,*) ' rn_lf_cutoff = 1.0/rdt' 998 994 ELSE 999 WRITE(numout,*) ' Namelist nam_vvl : z-tilde to zstar restoration timescale (days)' 1000 WRITE(numout,*) ' rn_rst_e3t = ', rn_rst_e3t 1001 WRITE(numout,*) ' Namelist nam_vvl : z-tilde cutoff frequency of low-pass filter (days)' 1002 WRITE(numout,*) ' rn_lf_cutoff = ', rn_lf_cutoff 995 WRITE(numout,*) ' z-tilde to zstar restoration timescale (days) rn_rst_e3t = ', rn_rst_e3t 996 WRITE(numout,*) ' z-tilde cutoff frequency of low-pass filter (days) rn_lf_cutoff = ', rn_lf_cutoff 1003 997 ENDIF 1004 WRITE(numout,*) ' Namelist nam_vvl : debug prints' 1005 WRITE(numout,*) ' ln_vvl_dbg = ', ln_vvl_dbg 998 WRITE(numout,*) ' debug prints flag ln_vvl_dbg = ', ln_vvl_dbg 1006 999 ENDIF 1007 1000 ! … … 1017 1010 IF(lwp) THEN ! Print the choice 1018 1011 WRITE(numout,*) 1019 IF( ln_vvl_zstar ) WRITE(numout,*) ' zstar vertical coordinate is used' 1020 IF( ln_vvl_ztilde ) WRITE(numout,*) ' ztilde vertical coordinate is used' 1021 IF( ln_vvl_layer ) WRITE(numout,*) ' layer vertical coordinate is used' 1022 IF( ln_vvl_ztilde_as_zstar ) WRITE(numout,*) ' to emulate a zstar coordinate' 1023 ! - ML - Option not developed yet 1024 ! IF( ln_vvl_kepe ) WRITE(numout,*) ' kinetic to potential energy transfer : option used' 1025 ! IF( .NOT. ln_vvl_kepe ) WRITE(numout,*) ' kinetic to potential energy transfer : option not used' 1012 IF( ln_vvl_zstar ) WRITE(numout,*) ' ==>>> zstar vertical coordinate is used' 1013 IF( ln_vvl_ztilde ) WRITE(numout,*) ' ==>>> ztilde vertical coordinate is used' 1014 IF( ln_vvl_layer ) WRITE(numout,*) ' ==>>> layer vertical coordinate is used' 1015 IF( ln_vvl_ztilde_as_zstar ) WRITE(numout,*) ' ==>>> to emulate a zstar coordinate' 1026 1016 ENDIF 1027 1017 !
Note: See TracChangeset
for help on using the changeset viewer.