Changeset 2715 for trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcana.F90
 Timestamp:
 20110330T17:58:35+02:00 (13 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcana.F90
r2548 r2715 27 27 PUBLIC sbc_gyre ! routine called in sbcmod module 28 28 29 ! !!* Namelist namsbc_ana *30 INTEGER :: nn_tau000 = 1 ! nb of timestep during which the surface stress31 ! ! increase from 0 to its nominal value32 REAL(wp) :: rn_utau0 = 0. e0! constant wind stress value in idirection33 REAL(wp) :: rn_vtau0 = 0. e0! constant wind stress value in jdirection34 REAL(wp) :: rn_qns0 = 0. e0! non solar heat flux35 REAL(wp) :: rn_qsr0 = 0. e0! solar heat flux36 REAL(wp) :: rn_emp0 = 0. e0! net freshwater flux29 ! !!* Namelist namsbc_ana * 30 INTEGER :: nn_tau000 = 1 ! nb of timestep during which the surface stress 31 ! ! increase from 0 to its nominal value 32 REAL(wp) :: rn_utau0 = 0._wp ! constant wind stress value in idirection 33 REAL(wp) :: rn_vtau0 = 0._wp ! constant wind stress value in jdirection 34 REAL(wp) :: rn_qns0 = 0._wp ! non solar heat flux 35 REAL(wp) :: rn_qsr0 = 0._wp ! solar heat flux 36 REAL(wp) :: rn_emp0 = 0._wp ! net freshwater flux 37 37 38 38 !! * Substitutions … … 42 42 !! NEMO/OPA 3.3 , NEMO Consortium (2010) 43 43 !! $Id$ 44 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 45 !! 46 44 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 45 !! 47 46 CONTAINS 48 47 … … 64 63 !! 65 64 INTEGER, INTENT(in) :: kt ! ocean time step 66 ! !67 REAL(wp) 68 REAL(wp) :: zrhoa = 1.22! Air density kg/m369 REAL(wp) :: zcdrag = 1.5e3! drag coefficient70 REAL(wp) 65 ! 66 REAL(wp) :: zfacto ! local scalar 67 REAL(wp) :: zrhoa = 1.22_wp ! Air density kg/m3 68 REAL(wp) :: zcdrag = 1.5e3_wp ! drag coefficient 69 REAL(wp) :: ztx, zty, zmod, zcoef ! temporary variables 71 70 !! 72 71 NAMELIST/namsbc_ana/ nn_tau000, rn_utau0, rn_vtau0, rn_qns0, rn_qsr0, rn_emp0 … … 75 74 IF( kt == nit000 ) THEN 76 75 ! 77 REWIND ( numnam )! Read Namelist namsbc : surface fluxes78 READ 76 REWIND( numnam ) ! Read Namelist namsbc : surface fluxes 77 READ ( numnam, namsbc_ana ) 79 78 ! 80 79 IF(lwp) WRITE(numout,*)' ' … … 88 87 IF(lwp) WRITE(numout,*)' net heat flux rn_emp0 = ', rn_emp0 , ' Kg/m2/s' 89 88 ! 90 nn_tau000 = MAX( nn_tau000, 1 ) ! must be >= 189 nn_tau000 = MAX( nn_tau000, 1 ) ! must be >= 1 91 90 ! 92 91 ENDIF 93 92 94 qns 95 qsr 96 emp 97 emps 93 qns (:,:) = rn_qns0 94 qsr (:,:) = rn_qsr0 95 emp (:,:) = rn_emp0 96 emps(:,:) = rn_emp0 98 97 99 98 ! Increase the surface stress to its nominal value during the first nn_tau000 timesteps … … 305 304 WRITE(numout,*)' ndastp = ', ndastp 306 305 WRITE(numout,*)' adatrj = ', adatrj 307 308 306 ENDIF 309 307 !
Note: See TracChangeset
for help on using the changeset viewer.