- Timestamp:
- 2017-09-27T11:11:01+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icedyn_adv.F90
r8562 r8564 89 89 ! !-----------------------! 90 90 CALL ice_dyn_adv_umx( nn_UMx, kt, u_ice, v_ice, & 91 & ato_i, v_i, v_s, s mv_i, oa_i, a_i, a_ip, v_ip, e_s, e_i )91 & ato_i, v_i, v_s, sv_i, oa_i, a_i, a_ip, v_ip, e_s, e_i ) 92 92 ! !-----------------------! 93 93 CASE( np_advPRA ) ! PRATHER scheme ! 94 94 ! !-----------------------! 95 95 CALL ice_dyn_adv_pra( kt, u_ice, v_ice, & 96 & ato_i, v_i, v_s, s mv_i, oa_i, a_i, a_ip, v_ip, e_s, e_i )96 & ato_i, v_i, v_s, sv_i, oa_i, a_i, a_ip, v_ip, e_s, e_i ) 97 97 END SELECT 98 98 … … 100 100 ! diagnostics 101 101 !------------ 102 diag_trp_ei (:,:) = SUM(SUM( e_i (:,:,1:nlay_i,:) - e_i_b(:,:,1:nlay_i,:), dim=4 ), dim=3 ) * r1_rdtice103 diag_trp_es (:,:) = SUM(SUM( e_s (:,:,1:nlay_s,:) - e_s_b(:,:,1:nlay_s,:), dim=4 ), dim=3 ) * r1_rdtice104 diag_trp_s mv(:,:) = SUM( smv_i(:,:,:) - smv_i_b(:,:,:) , dim=3 ) * r1_rdtice105 diag_trp_vi (:,:) = SUM( v_i (:,:,:) - v_i_b(:,:,:) , dim=3 ) * r1_rdtice106 diag_trp_vs (:,:) = SUM( v_s (:,:,:) - v_s_b(:,:,:) , dim=3 ) * r1_rdtice102 diag_trp_ei(:,:) = SUM(SUM( e_i (:,:,1:nlay_i,:) - e_i_b (:,:,1:nlay_i,:), dim=4 ), dim=3 ) * r1_rdtice 103 diag_trp_es(:,:) = SUM(SUM( e_s (:,:,1:nlay_s,:) - e_s_b (:,:,1:nlay_s,:), dim=4 ), dim=3 ) * r1_rdtice 104 diag_trp_sv(:,:) = SUM( sv_i(:,:,:) - sv_i_b(:,:,:) , dim=3 ) * r1_rdtice 105 diag_trp_vi(:,:) = SUM( v_i (:,:,:) - v_i_b (:,:,:) , dim=3 ) * r1_rdtice 106 diag_trp_vs(:,:) = SUM( v_s (:,:,:) - v_s_b (:,:,:) , dim=3 ) * r1_rdtice 107 107 IF( iom_use('icetrp') ) CALL iom_put( "icetrp" , diag_trp_vi ) ! ice volume transport 108 108 IF( iom_use('snwtrp') ) CALL iom_put( "snwtrp" , diag_trp_vs ) ! snw volume transport 109 IF( iom_use('saltrp') ) CALL iom_put( "saltrp" , diag_trp_s mv * rhoic )! salt content transport109 IF( iom_use('saltrp') ) CALL iom_put( "saltrp" , diag_trp_sv * rhoic ) ! salt content transport 110 110 IF( iom_use('deitrp') ) CALL iom_put( "deitrp" , diag_trp_ei ) ! advected ice enthalpy (W/m2) 111 111 IF( iom_use('destrp') ) CALL iom_put( "destrp" , diag_trp_es ) ! advected snw enthalpy (W/m2)
Note: See TracChangeset
for help on using the changeset viewer.