Changeset 886 for branches/dev_001_SBC/NEMO/LIM_SRC_3/limthd_dh.F90
- Timestamp:
- 2008-04-11T11:24:17+02:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/dev_001_SBC/NEMO/LIM_SRC_3/limthd_dh.F90
r884 r886 16 16 USE phycst ! physical constants (OCE directory) 17 17 USE ice_oce ! ice variables 18 USE sbc_oce ! Surface boundary condition: ocean fields 18 19 USE thd_ice 19 20 USE iceini … … 338 339 zjj = ( npb(ji) - 1 ) / jpi + 1 339 340 zfsalt_melt(ji) = zfsalt_melt(ji) + & 340 ( sss_ io(zji,zjj) - sm_i_b(ji) ) *&341 ( sss_m(zji,zjj) - sm_i_b(ji) ) * & 341 342 a_i_b(ji) * & 342 343 MIN( zdeltah(ji,jk) , 0.0 ) * rhoic / rdt_ice … … 368 369 WRITE(numout,*) ' qlbbq_1d: ', qlbbq_1d(ji) 369 370 WRITE(numout,*) ' s_i_new : ', s_i_new(ji) 370 WRITE(numout,*) ' sss_ io : ', sss_io(zji,zjj)371 WRITE(numout,*) ' sss_m : ', sss_m(zji,zjj) 371 372 ENDIF 372 373 … … 494 495 zswi2 * 0.26 / & 495 496 ( 0.26 + 0.74 * EXP ( - 724300.0 * zgrr ) ) 496 zds = zfracs*sss_ io(zji,zjj) - s_i_new(ji)497 s_i_new(ji) = zfracs * sss_ io(zji,zjj)497 zds = zfracs*sss_m(zji,zjj) - s_i_new(ji) 498 s_i_new(ji) = zfracs * sss_m(zji,zjj) 498 499 ENDIF ! fc_bo_i 499 500 END DO ! ji … … 567 568 zjj = ( npb(ji) - 1 ) / jpi + 1 568 569 zfsalt_melt(ji) = zfsalt_melt(ji) + & 569 ( sss_ io(zji,zjj) - sm_i_b(ji) ) *&570 ( sss_m(zji,zjj) - sm_i_b(ji) ) * & 570 571 a_i_b(ji) * & 571 572 MIN( zdeltah(ji,jk) , 0.0 ) * rhoic / rdt_ice … … 596 597 WRITE(numout,*) ' qlbbq_1d: ', qlbbq_1d(ji) 597 598 WRITE(numout,*) ' s_i_new : ', s_i_new(ji) 598 WRITE(numout,*) ' sss_ io : ', sss_io(zji,zjj)599 WRITE(numout,*) ' sss_m : ', sss_m(zji,zjj) 599 600 WRITE(numout,*) ' dh_i_bott : ', dh_i_bott(ji) 600 601 WRITE(numout,*) ' innermelt : ', innermelt(ji) … … 701 702 fseqv_1d(ji) = fseqv_1d(ji) + zihgnew * zfsalt_melt(ji) + & 702 703 (1.0 - zihgnew) * rdmicif_1d(ji) * & 703 ( sss_ io(zji,zjj) - sm_i_b(ji) ) / rdt_ice704 ( sss_m(zji,zjj) - sm_i_b(ji) ) / rdt_ice 704 705 ! new lines 705 706 IF ( num_sal .EQ. 4 ) & 706 707 fseqv_1d(ji) = fseqv_1d(ji) + zihgnew * zfsalt_melt(ji) + & 707 708 (1.0 - zihgnew) * rdmicif_1d(ji) * & 708 ( sss_ io(zji,zjj) - bulk_sal ) / rdt_ice709 ( sss_m(zji,zjj) - bulk_sal ) / rdt_ice 709 710 ! Heat flux 710 711 ! excessive bottom ablation energy (fsup) - 0 except if jpl = 1 … … 775 776 776 777 zsm_snowice = ( rhoic - rhosn ) / rhoic * & 777 sss_ io(zji,zjj)778 sss_m(zji,zjj) 778 779 779 780 IF ( num_sal .NE. 2 ) zsm_snowice = sm_i_b(ji) … … 781 782 IF ( num_sal .NE. 4 ) & 782 783 fseqv_1d(ji) = fseqv_1d(ji) + & 783 ( sss_ io(zji,zjj) - zsm_snowice ) * &784 ( sss_m(zji,zjj) - zsm_snowice ) * & 784 785 a_i_b(ji) * & 785 786 ( zhgnew(ji) - ht_i_b(ji) ) * rhoic / rdt_ice … … 787 788 IF ( num_sal .EQ. 4 ) & 788 789 fseqv_1d(ji) = fseqv_1d(ji) + & 789 ( sss_ io(zji,zjj) - bulk_sal ) * &790 ( sss_m(zji,zjj) - bulk_sal ) * & 790 791 a_i_b(ji) * & 791 792 ( zhgnew(ji) - ht_i_b(ji) ) * rhoic / rdt_ice … … 804 805 rdmicif_1d(ji) = rdmicif_1d(ji) + a_i_b(ji) & 805 806 * ( zhgnew(ji) - ht_i_b(ji) ) * rhoic & 806 + ( zhnnew - ht_s_b(ji) ) * rhosn )807 + ( zhnnew - ht_s_b(ji) ) * rhosn 807 808 #endif 808 809 ! Actualize new snow and ice thickness.
Note: See TracChangeset
for help on using the changeset viewer.