Changeset 3545 for branches/2012
- Timestamp:
- 2012-11-14T10:04:08+01:00 (12 years ago)
- Location:
- branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/LIM_SRC_3/limdia.F90
r3525 r3545 238 238 ij0 = 136 ; ij1 = 136 ! Fram strait : 83 = area export 239 239 ii0 = 134 ; ii1 = 138 ! 84 = volume export 240 DO jj = mj0(ij0) :mj1(ij1)241 DO ji = mi0(ii0) :mi1(ii1)240 DO jj = mj0(ij0),mj1(ij1) 241 DO ji = mi0(ii0),mi1(ii1) 242 242 vinfor(83) = vinfor(83) - v_ice(ji,jj) * e1t(ji,jj)*at_i(ji,jj)*rdt_ice * 1.e-12_wp 243 243 vinfor(84) = vinfor(84) - v_ice(ji,jj) * e1t(ji,jj)*vt_i(ji,jj)*rdt_ice * 1.e-12_wp 244 244 END DO 245 245 END DO 246 END SELECT 246 247 !!gm just above, this is NOT the correct way of evaluating the transport ! 247 248 !!gm mass of snow is missing and v_ice should be the mean between jj and jj+1 -
branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90
r3524 r3545 103 103 INTEGER :: iflt, ial , iadv , ifral, ifrdv ! - - 104 104 REAL(wp) :: zinda, zemp, zemp_snow, zfmm ! local scalars 105 REAL(wp) :: zemp_snw ! - - 105 106 REAL(wp) :: zfcm1 , zfcm2 ! - - 106 107 REAL(wp), POINTER, DIMENSION(:,:,:) :: zalb, zalbp ! 2D/3D workspace -
branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/LIM_SRC_3/limthd_dh.F90
r3524 r3545 593 593 ! 594 594 ! Salt flux 595 sfx_thd_1d(ji) = sfx_thd_1d(ji) + zihgnew * zsfx_melt(ji) 596 & - (1.0 - zihgnew) * zfmass_i (ji) * sm_i_b(ji) )* r1_rdtice595 sfx_thd_1d(ji) = sfx_thd_1d(ji) + zihgnew * zsfx_melt(ji) & 596 & - (1.0 - zihgnew) * zfmass_i (ji) * sm_i_b(ji) * r1_rdtice 597 597 ! 598 598 ! Heat flux -
branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim.F90
r3517 r3545 178 178 d_oa_i_thd (:,:,:) = 0._wp ; d_oa_i_trp (:,:,:) = 0._wp 179 179 ! 180 fseqv(:,:) = 0._wp181 fsbri (:,:) = 0._wp ; fsalt_res(:,:) = 0._wp182 f salt_rpo(:,:) = 0._wp183 fhmec (:,:) = 0._wp ; fhbri (:,:) = 0._wp184 fmmec (:,:) = 0._wp ; fheat_res(:,:) = 0._wp185 f heat_rpo(:,:) = 0._wp ; focea2D(:,:) = 0._wp186 fsup2D 180 sfx (:,:) = 0._wp 181 sfx_bri(:,:) = 0._wp ; sfx_mec (:,:) = 0._wp ; sfx_res (:,:) = 0._wp 182 fhbri (:,:) = 0._wp ; fheat_mec(:,:) = 0._wp ; fheat_res(:,:) = 0._wp 183 fhmec (:,:) = 0._wp ; 184 fmmec (:,:) = 0._wp 185 focea2D(:,:) = 0._wp 186 fsup2D (:,:) = 0._wp 187 187 ! 188 188 diag_sni_gr(:,:) = 0._wp ; diag_lat_gr(:,:) = 0._wp … … 428 428 WRITE(numout,*) ' fdtcn : ', fdtcn(ji,jj) 429 429 WRITE(numout,*) ' fhmec : ', fhmec(ji,jj) 430 WRITE(numout,*) ' fheat_ rpo : ', fheat_rpo(ji,jj)430 WRITE(numout,*) ' fheat_mec : ', fheat_mec(ji,jj) 431 431 WRITE(numout,*) ' fheat_res : ', fheat_res(ji,jj) 432 432 WRITE(numout,*) ' fhbri : ', fhbri(ji,jj) … … 605 605 WRITE(numout,*) ' - Heat / FW fluxes ' 606 606 WRITE(numout,*) ' ~~~~~~~~~~~~~~~~ ' 607 ! WRITE(numout,*) ' fsbri : ', fsbri(ki,kj)608 ! WRITE(numout,*) ' fseqv : ', fseqv(ki,kj)607 ! WRITE(numout,*) ' sfx_bri : ', sfx_bri (ki,kj) 608 ! WRITE(numout,*) ' sfx : ', sfx (ki,kj) 609 609 ! WRITE(numout,*) ' fsalt_res : ', fsalt_res(ki,kj) 610 WRITE(numout,*) ' fmmec : ', fmmec (ki,kj)611 WRITE(numout,*) ' fhmec : ', fhmec (ki,kj)612 WRITE(numout,*) ' fhbri : ', fhbri (ki,kj)613 WRITE(numout,*) ' fheat_ rpo : ', fheat_rpo(ki,kj)610 WRITE(numout,*) ' fmmec : ', fmmec (ki,kj) 611 WRITE(numout,*) ' fhmec : ', fhmec (ki,kj) 612 WRITE(numout,*) ' fhbri : ', fhbri (ki,kj) 613 WRITE(numout,*) ' fheat_mec : ', fheat_mec(ki,kj) 614 614 WRITE(numout,*) 615 615 WRITE(numout,*) ' sst : ', sst_m(ki,kj) … … 620 620 WRITE(numout,*) ' utau_ice : ', utau_ice(ki,kj) 621 621 WRITE(numout,*) ' vtau_ice : ', vtau_ice(ki,kj) 622 WRITE(numout,*) ' utau : ', utau (ki,kj)623 WRITE(numout,*) ' vtau : ', vtau (ki,kj)624 WRITE(numout,*) ' oc. vel. u : ', u_oce (ki,kj)625 WRITE(numout,*) ' oc. vel. v : ', v_oce (ki,kj)622 WRITE(numout,*) ' utau : ', utau (ki,kj) 623 WRITE(numout,*) ' vtau : ', vtau (ki,kj) 624 WRITE(numout,*) ' oc. vel. u : ', u_oce (ki,kj) 625 WRITE(numout,*) ' oc. vel. v : ', v_oce (ki,kj) 626 626 ENDIF 627 627 … … 639 639 WRITE(numout,*) 640 640 WRITE(numout,*) ' - Heat fluxes at bottom interface ***' 641 WRITE(numout,*) ' qsr 642 WRITE(numout,*) ' qns 641 WRITE(numout,*) ' qsr : ', qsr(ki,kj) 642 WRITE(numout,*) ' qns : ', qns(ki,kj) 643 643 WRITE(numout,*) 644 644 WRITE(numout,*) ' - Salt fluxes at bottom interface ***' 645 WRITE(numout,*) ' sfx : ', sfx (ki,kj) 646 WRITE(numout,*) ' emp : ', emp(ki,kj) 647 WRITE(numout,*) ' fsbri : ', fsbri(ki,kj) 648 WRITE(numout,*) ' fseqv : ', fseqv(ki,kj) 649 WRITE(numout,*) ' fsalt_res : ', fsalt_res(ki,kj) 650 WRITE(numout,*) ' fsalt_rpo : ', fsalt_rpo(ki,kj) 645 WRITE(numout,*) ' emp : ', emp (ki,kj) 646 WRITE(numout,*) ' sfx_bri : ', sfx_bri(ki,kj) 647 WRITE(numout,*) ' sfx : ', sfx (ki,kj) 648 WRITE(numout,*) ' sfx_res : ', sfx_res(ki,kj) 649 WRITE(numout,*) ' sfx_mec : ', sfx_mec(ki,kj) 651 650 WRITE(numout,*) ' - Heat fluxes at bottom interface ***' 652 WRITE(numout,*) ' fheat_res 651 WRITE(numout,*) ' fheat_res : ', fheat_res(ki,kj) 653 652 WRITE(numout,*) 654 653 WRITE(numout,*) ' - Momentum fluxes '
Note: See TracChangeset
for help on using the changeset viewer.