Changeset 4602
- Timestamp:
- 2014-03-28T11:32:43+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/OPA_SRC/SBC/sbcana.F90
r3294 r4602 89 89 nn_tau000 = MAX( nn_tau000, 1 ) ! must be >= 1 90 90 ! 91 utau(:,:) = rn_utau0 92 vtau(:,:) = rn_vtau0 93 taum(:,:) = SQRT ( rn_utau0 * rn_utau0 + rn_vtau0 * rn_vtau0 ) 94 wndm(:,:) = SQRT ( taum(1,1) / ( zrhoa * zcdrag ) ) 95 ! 91 96 qns (:,:) = rn_qns0 92 97 qsr (:,:) = rn_qsr0 … … 94 99 emps(:,:) = rn_emp0 95 100 ! 96 utau(:,:) = rn_utau097 vtau(:,:) = rn_vtau098 taum(:,:) = SQRT ( rn_utau0 * rn_utau0 + rn_vtau0 * rn_vtau0 )99 wndm(:,:) = SQRT ( taum(1,1) / ( zrhoa * zcdrag ) )100 !101 101 ENDIF 102 102 103 104 ! Increase the surface stress to its nominal value during the first nn_tau000 time-steps 105 IF( kt <= nn_tau000 ) THEN 106 zfacto = 0.5 * ( 1. - COS( rpi * FLOAT( kt ) / FLOAT( nn_tau000 ) ) ) 107 zcoef = 1. / ( zrhoa * zcdrag ) 108 ztx = zfacto * rn_utau0 109 zty = zfacto * rn_vtau0 110 zmod = SQRT( ztx * ztx + zty * zty ) 111 utau(:,:) = ztx 112 vtau(:,:) = zty 113 taum(:,:) = zmod 114 zmod = SQRT( zmod * zcoef ) 115 wndm(:,:) = zmod 103 IF( MOD( kt - 1, nn_fsbc ) == 0 ) THEN 104 ! 105 ! Increase the surface stress to its nominal value during the first nn_tau000 time-steps 106 IF( kt <= nn_tau000 ) THEN 107 zfacto = 0.5 * ( 1. - COS( rpi * FLOAT( kt ) / FLOAT( nn_tau000 ) ) ) 108 zcoef = 1. / ( zrhoa * zcdrag ) 109 ztx = zfacto * rn_utau0 110 zty = zfacto * rn_vtau0 111 zmod = SQRT( ztx * ztx + zty * zty ) 112 utau(:,:) = ztx 113 vtau(:,:) = zty 114 taum(:,:) = zmod 115 zmod = SQRT( zmod * zcoef ) 116 wndm(:,:) = zmod 117 ENDIF 118 ! 119 ! retore heat and fresh water fluxes as they may have been changed by sbcssr module 120 qns (:,:) = rn_qns0 121 qsr (:,:) = rn_qsr0 122 emp (:,:) = rn_emp0 123 emps(:,:) = rn_emp0 124 ! 116 125 ENDIF 117 126 !
Note: See TracChangeset
for help on using the changeset viewer.