- Timestamp:
- 2012-10-26T12:13:21+02:00 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/LIM_SRC_3/limthd.F90
r3419 r3517 110 110 e_i(ji,jj,jk,jl) = e_i(ji,jj,jk,jl) / ( area(ji,jj) * MAX( v_i(ji,jj,jl) , epsi06 ) ) * nlay_i 111 111 !0 if no ice and 1 if yes 112 zindb = 1.0 - MAX ( 0.0 , SIGN ( 1.0 , - ht_i(ji,jj,jl) ))112 zindb = 1.0 - MAX( 0.0 , SIGN( 1.0 , - ht_i(ji,jj,jl) ) ) 113 113 !convert units ! very important that this line is here 114 114 e_i(ji,jj,jk,jl) = e_i(ji,jj,jk,jl) * unit_fac * zindb … … 122 122 e_s(ji,jj,jk,jl) = e_s(ji,jj,jk,jl) / ( area(ji,jj) * MAX( v_s(ji,jj,jl) , epsi06 ) ) * nlay_s 123 123 !0 if no ice and 1 if yes 124 zindb = 1.0 - MAX ( 0.0 , SIGN ( 1.0 , - ht_s(ji,jj,jl) ))124 zindb = 1.0 - MAX( 0.0 , SIGN( 1.0 , - ht_s(ji,jj,jl) ) ) 125 125 !convert units ! very important that this line is here 126 126 e_s(ji,jj,jk,jl) = e_s(ji,jj,jk,jl) * unit_fac * zindb … … 284 284 CALL tab_2d_1d( nbpb, fbif_1d (1:nbpb), fbif , jpi, jpj, npb(1:nbpb) ) 285 285 CALL tab_2d_1d( nbpb, qldif_1d (1:nbpb), qldif , jpi, jpj, npb(1:nbpb) ) 286 CALL tab_2d_1d( nbpb, rdm icif_1d (1:nbpb), rdm_ice , jpi, jpj, npb(1:nbpb) )287 CALL tab_2d_1d( nbpb, rdm snif_1d (1:nbpb), rdm_snw , jpi, jpj, npb(1:nbpb) )286 CALL tab_2d_1d( nbpb, rdm_ice_1d (1:nbpb), rdm_ice , jpi, jpj, npb(1:nbpb) ) 287 CALL tab_2d_1d( nbpb, rdm_snw_1d (1:nbpb), rdm_snw , jpi, jpj, npb(1:nbpb) ) 288 288 CALL tab_2d_1d( nbpb, dmgwi_1d (1:nbpb), dmgwi , jpi, jpj, npb(1:nbpb) ) 289 289 CALL tab_2d_1d( nbpb, qlbbq_1d (1:nbpb), zqlbsbq , jpi, jpj, npb(1:nbpb) ) … … 352 352 CALL tab_1d_2d( nbpb, qldif , npb, qldif_1d (1:nbpb), jpi, jpj ) 353 353 CALL tab_1d_2d( nbpb, qfvbq , npb, qfvbq_1d (1:nbpb), jpi, jpj ) 354 CALL tab_1d_2d( nbpb, rdm_ice, npb, rdm icif_1d(1:nbpb), jpi, jpj )355 CALL tab_1d_2d( nbpb, rdm_snw, npb, rdm snif_1d(1:nbpb), jpi, jpj )354 CALL tab_1d_2d( nbpb, rdm_ice, npb, rdm_ice_1d(1:nbpb), jpi, jpj ) 355 CALL tab_1d_2d( nbpb, rdm_snw, npb, rdm_snw_1d(1:nbpb), jpi, jpj ) 356 356 CALL tab_1d_2d( nbpb, dmgwi , npb, dmgwi_1d (1:nbpb), jpi, jpj ) 357 357 CALL tab_1d_2d( nbpb, rdvosif, npb, dvsbq_1d (1:nbpb), jpi, jpj ) … … 419 419 !-------------------------------------------- 420 420 d_v_i_thd(:,:,:) = v_i (:,:,:) - old_v_i(:,:,:) ! ice volumes 421 dv_dt_thd(:,:,:) = d_v_i_thd(:,:,:) / rdt_ice * 86400.0421 dv_dt_thd(:,:,:) = d_v_i_thd(:,:,:) * r1_rdtice * rday 422 422 423 423 IF( con_i ) fbif(:,:) = fbif(:,:) + zqlbsbq(:,:) … … 488 488 ! 489 489 IF(lwp) WRITE(numout,*) ' lim_thd_glohec ' 490 IF(lwp) WRITE(numout,*) ' qt_i_in : ', eti(jiindex_1d,jl) / rdt_ice491 IF(lwp) WRITE(numout,*) ' qt_s_in : ', ets(jiindex_1d,jl) / rdt_ice492 IF(lwp) WRITE(numout,*) ' qt_in : ', ( eti(jiindex_1d,jl) + ets(jiindex_1d,jl) ) / rdt_ice490 IF(lwp) WRITE(numout,*) ' qt_i_in : ', eti(jiindex_1d,jl) * r1_rdtice 491 IF(lwp) WRITE(numout,*) ' qt_s_in : ', ets(jiindex_1d,jl) * r1_rdtice 492 IF(lwp) WRITE(numout,*) ' qt_in : ', ( eti(jiindex_1d,jl) + ets(jiindex_1d,jl) ) * r1_rdtice 493 493 ! 494 494 END SUBROUTINE lim_thd_glohec … … 538 538 !-------------------- 539 539 DO ji = kideb, kiut 540 cons_error(ji,jl) = ABS( dq_i(ji,jl) / rdt_ice + sum_fluxq(ji,jl) )540 cons_error(ji,jl) = ABS( dq_i(ji,jl) * r1_rdtice + sum_fluxq(ji,jl) ) 541 541 END DO 542 542 … … 597 597 WRITE(numout,*) ' cons_error : ', cons_error(ji,jl) 598 598 WRITE(numout,*) ' surf_error : ', surf_error(ji,jl) 599 WRITE(numout,*) ' dq_i : ', - dq_i(ji,jl) / rdt_ice599 WRITE(numout,*) ' dq_i : ', - dq_i(ji,jl) * r1_rdtice 600 600 WRITE(numout,*) ' Fdt : ', sum_fluxq(ji,jl) 601 601 WRITE(numout,*) … … 631 631 WRITE(numout,*) 632 632 WRITE(numout,*) ' Layer by layer ... ' 633 WRITE(numout,*) ' dq_snow : ', ( qt_s_fin(ji,jl) - qt_s_in(ji,jl) ) / rdt_ice633 WRITE(numout,*) ' dq_snow : ', ( qt_s_fin(ji,jl) - qt_s_in(ji,jl) ) * r1_rdtice 634 634 WRITE(numout,*) ' dfc_snow : ', fc_s(ji,1) - fc_s(ji,0) 635 635 DO jk = 1, nlay_i 636 636 WRITE(numout,*) ' layer : ', jk 637 WRITE(numout,*) ' dq_ice : ', dq_i_layer(ji,jk) / rdt_ice637 WRITE(numout,*) ' dq_ice : ', dq_i_layer(ji,jk) * r1_rdtice 638 638 WRITE(numout,*) ' radab : ', radab(ji,jk) 639 639 WRITE(numout,*) ' dfc_i : ', fc_i(ji,jk) - fc_i(ji,jk-1) … … 681 681 fatm (ji,jl) = qnsr_ice_1d(ji) + qsr_ice_1d(ji) ! total heat flux 682 682 sum_fluxq (ji,jl) = fatm(ji,jl) + fbif_1d(ji) - ftotal_fin(ji) - fstroc(zji,zjj,jl) 683 cons_error(ji,jl) = ABS( dq_i(ji,jl) / rdt_ice + sum_fluxq(ji,jl) )683 cons_error(ji,jl) = ABS( dq_i(ji,jl) * r1_rdtice + sum_fluxq(ji,jl) ) 684 684 END DO 685 685 … … 688 688 !-------------------- 689 689 DO ji = kideb, kiut 690 cons_error(ji,jl) = ABS( dq_i(ji,jl) / rdt_ice + sum_fluxq(ji,jl) )690 cons_error(ji,jl) = ABS( dq_i(ji,jl) * r1_rdtice + sum_fluxq(ji,jl) ) 691 691 END DO 692 692 … … 722 722 WRITE(numout,*) ' * ' 723 723 WRITE(numout,*) ' Ftotal : ', sum_fluxq(ji,jl) 724 WRITE(numout,*) ' dq_t : ', - dq_i(ji,jl) / rdt_ice725 WRITE(numout,*) ' dq_i : ', - ( qt_i_fin(ji,jl) - qt_i_in(ji,jl) ) / rdt_ice726 WRITE(numout,*) ' dq_s : ', - ( qt_s_fin(ji,jl) - qt_s_in(ji,jl) ) / rdt_ice724 WRITE(numout,*) ' dq_t : ', - dq_i(ji,jl) * r1_rdtice 725 WRITE(numout,*) ' dq_i : ', - ( qt_i_fin(ji,jl) - qt_i_in(ji,jl) ) * r1_rdtice 726 WRITE(numout,*) ' dq_s : ', - ( qt_s_fin(ji,jl) - qt_s_in(ji,jl) ) * r1_rdtice 727 727 WRITE(numout,*) ' cons_error : ', cons_error(ji,jl) 728 728 WRITE(numout,*) ' * ' … … 734 734 WRITE(numout,*) ' * ' 735 735 WRITE(numout,*) ' Heat contents --- : ' 736 WRITE(numout,*) ' qt_s_in : ', qt_s_in(ji,jl) / rdt_ice737 WRITE(numout,*) ' qt_i_in : ', qt_i_in(ji,jl) / rdt_ice738 WRITE(numout,*) ' qt_in : ', ( qt_i_in(ji,jl) + qt_s_in(ji,jl) ) / rdt_ice739 WRITE(numout,*) ' qt_s_fin : ', qt_s_fin(ji,jl) / rdt_ice740 WRITE(numout,*) ' qt_i_fin : ', qt_i_fin(ji,jl) / rdt_ice741 WRITE(numout,*) ' qt_fin : ', ( qt_i_fin(ji,jl) + qt_s_fin(ji,jl) ) / rdt_ice736 WRITE(numout,*) ' qt_s_in : ', qt_s_in(ji,jl) * r1_rdtice 737 WRITE(numout,*) ' qt_i_in : ', qt_i_in(ji,jl) * r1_rdtice 738 WRITE(numout,*) ' qt_in : ', ( qt_i_in(ji,jl) + qt_s_in(ji,jl) ) * r1_rdtice 739 WRITE(numout,*) ' qt_s_fin : ', qt_s_fin(ji,jl) * r1_rdtice 740 WRITE(numout,*) ' qt_i_fin : ', qt_i_fin(ji,jl) * r1_rdtice 741 WRITE(numout,*) ' qt_fin : ', ( qt_i_fin(ji,jl) + qt_s_fin(ji,jl) ) * r1_rdtice 742 742 WRITE(numout,*) ' * ' 743 743 WRITE(numout,*) ' Ice variables --- : '
Note: See TracChangeset
for help on using the changeset viewer.