- Timestamp:
- 2012-05-17T18:33:12+02:00 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/OPA_SRC/TRA/trasbc.F90
r3294 r3396 60 60 !! at the surface by evaporation, precipitations and runoff (E-P-R); 61 61 !! (3) Fwe, tracer carried with the water that is exchanged. 62 !! - salinity : salt flux only due to freezing/melting 63 !! sa = sa + fsalt / rau0 / e3t for k=1 62 64 !! 63 65 !! Fext, flux through the air-sea interface for temperature and salt: … … 84 86 !! (Tp P - Te E) + SST (P-E) = 0 when Tp=Te=SST 85 87 !! - salinity : evaporation, precipitation and runoff 86 !! water has a zero salinity (Fwe=0), thus only Fwi remains: 87 !! sa = sa + emp * sn / e3t for k=1 88 !! water has a zero salinity but there is a salt flux due to 89 !! freezing/melting, thus: 90 !! sa = sa + emp * sn / rau0 / e3t for k=1 91 !! + fsalt / rau0 / e3t 88 92 !! where emp, the surface freshwater budget (evaporation minus 89 93 !! precipitation minus runoff) given in kg/m2/s is divided … … 109 113 !! 110 114 INTEGER :: ji, jj, jk, jn ! dummy loop indices 111 REAL(wp) :: zfact, z1_e3t, z srau, zdep115 REAL(wp) :: zfact, z1_e3t, zdep 112 116 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztrdt, ztrds 113 117 !!---------------------------------------------------------------------- … … 120 124 IF(lwp) WRITE(numout,*) '~~~~~~~ ' 121 125 ENDIF 122 123 zsrau = 1. / rau0 ! initialization124 126 125 127 IF( l_trdtra ) THEN !* Save ta and sa trends … … 163 165 ! evaporation, precipitation and qns, but not river runoff 164 166 165 IF( lk_vvl ) THEN ! Variable Volume case 167 IF( lk_vvl ) THEN ! Variable Volume case ==>> heat content of mass flux is in qns 166 168 DO jj = 1, jpj 167 169 DO ji = 1, jpi 168 ! temperature : heat flux + cooling/heating effet of EMP flux 169 sbc_tsc(ji,jj,jp_tem) = ro0cpr * qns(ji,jj) - zsrau * emp(ji,jj) * tsn(ji,jj,1,jp_tem) 170 ! concent./dilut. effect due to sea-ice melt/formation and (possibly) SSS restoration 171 sbc_tsc(ji,jj,jp_sal) = ( emps(ji,jj) - emp(ji,jj) ) * zsrau * tsn(ji,jj,1,jp_sal) 170 sbc_tsc(ji,jj,jp_tem) = r1_rau0_rcp * qns(ji,jj) ! non solar heat flux 171 sbc_tsc(ji,jj,jp_sal) = r1_rau0 * emps(ji,jj) ! salt flux (freezing/melting) 172 172 END DO 173 173 END DO … … 176 176 DO ji = fs_2, fs_jpim1 ! vector opt. 177 177 ! temperature : heat flux 178 sbc_tsc(ji,jj,jp_tem) = ro0cpr * qns(ji,jj) 178 sbc_tsc(ji,jj,jp_tem) = r1_rau0_rcp * qns(ji,jj) & ! non solar heat flux 179 & + r1_rau0 * emp(ji,jj) * tsn(ji,jj,1,jp_tem) ! concent./dilut. effect 179 180 ! salinity : salt flux + concent./dilut. effect (both in emps) 180 sbc_tsc(ji,jj,jp_sal) = zsrau * emps(ji,jj) * tsn(ji,jj,1,jp_sal) 181 sbc_tsc(ji,jj,jp_sal) = r1_rau0 * ( emps(ji,jj) & ! salt flux (freezing/melting) 182 & + emp (ji,jj) * tsn(ji,jj,1,jp_sal) ) ! concent./dilut. effect 181 183 END DO 182 184 END DO
Note: See TracChangeset
for help on using the changeset viewer.