Changeset 12475


Ignore:
Timestamp:
2020-02-27T12:27:50+01:00 (8 months ago)
Author:
dancopsey
Message:
  • Add more print statements.
  • Move away from using snow to ice diagnostics and use a new snow to pond one instead.
Location:
NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints
Files:
17 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/cfgs/SHARED/field_def_nemo-ice.xml

    r12449 r12475  
    221221          <field id="dmibog"       long_name="sea-ice mass change through basal growth"                standard_name="tendency_of_sea_ice_amount_due_to_congelation_ice_accumulation"          unit="kg/m2/s" /> 
    222222          <field id="dmisni"       long_name="sea-ice mass change through snow-to-ice conversion"      standard_name="tendency_of_sea_ice_amount_due_to_snow_conversion"                       unit="kg/m2/s" /> 
     223          <field id="dmisnp"       long_name="sea-ice mass change through snow melting into ponds"     standard_name="tendency_of_sea_ice_amount_due_to_snow_melt_into_ponds"                  unit="kg/m2/s" /> 
    223224          <field id="dmtsub"       long_name="snow and sea ice mass change through sublimation"        standard_name="tendency_of_snow_and_ice_amount_due_to_sublimation"                      unit="kg/m2/s" /> 
    224225          <field id="dmssub"       long_name="sea-ice mass change through evaporation and sublimation" standard_name="water_evaporation_flux"                                                  unit="kg/m2/s" /> 
  • NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/ice.F90

    r12449 r12475  
    245245   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   wfx_ice     !: mass flux from ice-ocean mass exchange                   [kg.m-2.s-1] 
    246246   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   wfx_sni     !: mass flux from snow ice growth component of wfx_ice      [kg.m-2.s-1] 
     247   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   wfx_snp     !: mass flux from snow melting into melt ponds              [kg.m-2.s-1] 
    247248   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   wfx_opw     !: mass flux from lateral ice growth component of wfx_ice   [kg.m-2.s-1] 
    248249   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   wfx_bog     !: mass flux from bottom ice growth component of wfx_ice    [kg.m-2.s-1] 
     
    415416         &      wfx_bog    (jpi,jpj) , wfx_dyn   (jpi,jpj) , wfx_bom(jpi,jpj) , wfx_sum(jpi,jpj) ,           & 
    416417         &      wfx_res    (jpi,jpj) , wfx_sni   (jpi,jpj) , wfx_opw(jpi,jpj) , wfx_spr(jpi,jpj) ,           & 
    417          &      afx_tot    (jpi,jpj) , rn_amax_2d(jpi,jpj) ,                                                 & 
     418         &      afx_tot    (jpi,jpj) , rn_amax_2d(jpi,jpj) , wfx_snp(jpi,jpj) ,                              & 
    418419         &      qsb_ice_bot(jpi,jpj) , qlead     (jpi,jpj) ,                                                 & 
    419420         &      sfx_res    (jpi,jpj) , sfx_bri   (jpi,jpj) , sfx_dyn(jpi,jpj) , sfx_sub(jpi,jpj) , sfx_lam(jpi,jpj) ,  & 
  • NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/ice1d.F90

    r12449 r12475  
    7777   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   wfx_sum_1d   
    7878   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   wfx_sni_1d  
     79   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   wfx_snp_1d  
    7980   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   wfx_opw_1d 
    8081   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   wfx_res_1d  
     
    203204         &      fhld_1d       (jpij) , wfx_sub_1d    (jpij) , wfx_bog_1d    (jpij) , wfx_bom_1d    (jpij) ,  & 
    204205         &      wfx_sum_1d    (jpij) , wfx_sni_1d    (jpij) , wfx_opw_1d    (jpij) , wfx_res_1d    (jpij) ,  & 
     206         &      wfx_snp_1d    (jpij) ,                                                                       & 
    205207         &      wfx_snw_sub_1d(jpij) , wfx_snw_dyn_1d(jpij) , wfx_ice_sub_1d(jpij) , wfx_err_sub_1d(jpij) ,  & 
    206208         &      wfx_lam_1d    (jpij) , wfx_dyn_1d    (jpij) , wfx_pnd_1d    (jpij) , dqns_ice_1d   (jpij) , evap_ice_1d (jpij) , & 
  • NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/icedyn.F90

    r12379 r12475  
    105105      ! 
    106106 
    107       write(numout,*)'ice_dyn 0: e_i = ',e_i (42,26,1,1) 
     107      write(numout,*)'ice_dyn 0:sv_i, e_i, t_i = ',sv_i(3,4,1), ' ', e_i(3,4,1,1), ' ', t_i(3,4,1,1) 
     108      write(numout,*)'ice_dyn 0:a_i, h_i, v_i = ',a_i(3,4,1), ' ', h_i(3,4,1), ' ', v_i(3,4,1) 
     109      write(numout,*)'ice_dyn 0:s_i = ',s_i(3,4,1) 
     110      write(numout,*)'ice_dyn 0: MEAN(sv_i) = ',SUM(sv_i)/SIZE(sv_i) 
     111      write(numout,*)'ice_dyn 0: MEAN(e_i) = ',SUM(e_i)/SIZE(e_i) 
     112      write(numout,*)'ice_dyn 0: MEAN(t_i) = ',SUM(t_i)/SIZE(t_i) 
     113      write(numout,*)'ice_dyn 0: MEAN(a_i) = ',SUM(a_i)/SIZE(a_i) 
     114      write(numout,*)'ice_dyn 0: MEAN(h_i) = ',SUM(h_i)/SIZE(h_i) 
     115      write(numout,*)'ice_dyn 0: MEAN(v_i) = ',SUM(v_i)/SIZE(v_i) 
     116      write(numout,*)'ice_dyn 0: MEAN(s_i) = ',SUM(s_i)/SIZE(s_i) 
     117      write(numout,*)'ice_dyn 0: MEAN(u_ice) = ',SUM(u_ice)/SIZE(u_ice) 
     118      write(numout,*)'ice_dyn 0: MEAN(e2u) = ',SUM(e2u)/SIZE(e2u) 
     119      write(numout,*)'ice_dyn 0: u_ice, e2u = ',u_ice(3,4), e2u(3,4) 
     120       
    108121      ! 
    109122      SELECT CASE( nice_dyn )          !-- Set which dynamics is running 
     
    111124      CASE ( np_dynALL )           !==  all dynamical processes  ==! 
    112125         ! 
    113          CALL ice_dyn_rhg   ( kt )                                          ! -- rheology   
     126         CALL ice_dyn_rhg   ( kt )                                          ! -- rheology  
     127          write(numout,*)'ice_dyn 1:sv_i, e_s = ',sv_i(3,4,1), ' ', e_i (3,4,1,1) 
     128          write(numout,*)'ice_dyn 1: u_ice, e2u = ',u_ice(3,4), e2u(3,4) 
    114129         CALL ice_dyn_adv   ( kt )                                          ! -- advection of ice 
     130          write(numout,*)'ice_dyn 2:sv_i, e_s = ',sv_i(3,4,1), ' ', e_i (3,4,1,1) 
    115131         CALL ice_dyn_rdgrft( kt )                                          ! -- ridging/rafting  
     132          write(numout,*)'ice_dyn 3:sv_i, e_s = ',sv_i(3,4,1), ' ', e_i (3,4,1,1) 
    116133         CALL ice_cor       ( kt , 1 )                                      ! -- Corrections 
    117134         ! 
     
    150167      END SELECT 
    151168 
    152       write(numout,*)'ice_dyn 1: e_i = ',e_i (42,26,1,1) 
     169      write(numout,*)'ice_dyn 4: e_i = ',e_i (42,26,1,1) 
    153170      ! 
    154171      ! 
     
    174191         ! 
    175192      ENDIF 
    176       write(numout,*)'ice_dyn 2: e_i = ',e_i (42,26,1,1) 
     193      write(numout,*)'ice_dyn 5: e_i = ',e_i (42,26,1,1) 
    177194      ! 
    178195      ! controls 
  • NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/icedyn_adv_umx.F90

    r12379 r12475  
    105105      ! 
    106106      IF( kt == nit000 .AND. lwp )   WRITE(numout,*) '-- ice_dyn_adv_umx: Ultimate-Macho advection scheme' 
     107 
     108      write(numout,*)'ice_dyn_adv_umx 1:psv_i, pe_s = ',psv_i(3,4,1), ' ', pe_i (3,4,1,1) 
    107109      ! 
    108110      ! --- Record max of the surrounding 9-pts ice thick. (for call Hbig) --- ! 
     
    142144      zdt = rdt_ice / REAL(icycle) 
    143145 
     146      write(numout,*)'ice_dyn_adv_umx 2: pu_ice, e2u = ',pu_ice(3,4), e2u(3,4) 
     147 
    144148      ! --- transport --- ! 
    145149      zudy(:,:) = pu_ice(:,:) * e2u(:,:) 
    146150      zvdx(:,:) = pv_ice(:,:) * e1v(:,:) 
     151 
     152      write(numout,*)'ice_dyn_adv_umx 2: zudy, zvdx, zu_cat, zv_cat, zcu_box, zcv_box = ',zudy(3,4), zvdx(3,4), zu_cat(3,4,1), zv_cat(3,4,1), zcu_box(3,4), zcv_box(3,4) 
     153 
    147154      ! 
    148155      ! setup transport for each ice cat  
     
    151158         zv_cat(:,:,jl) = zvdx(:,:) 
    152159      END DO 
     160 
     161      write(numout,*)'ice_dyn_adv_umx 3: zudy, zvdx, zu_cat, zv_cat, zcu_box, zcv_box = ',zudy(3,4), zvdx(3,4), zu_cat(3,4,1), zv_cat(3,4,1), zcu_box(3,4), zcv_box(3,4) 
     162 
    153163      ! 
    154164      ! --- define velocity for advection: u*grad(H) --- ! 
     
    166176         END DO 
    167177      END DO 
     178 
     179      write(numout,*)'ice_dyn_adv_umx 4: zudy, zvdx, zu_cat, zv_cat, zcu_box, zcv_box = ',zudy(3,4), zvdx(3,4), zu_cat(3,4,1), zv_cat(3,4,1), zcu_box(3,4), zcv_box(3,4) 
    168180 
    169181      !---------------! 
     
    209221         CALL adv_umx( zamsk, kn_umx, jt, kt, zdt, zudy, zvdx, zu_cat , zv_cat , zcu_box, zcv_box, & 
    210222            &                                      pa_i, pa_i, zua_ups, zva_ups, zua_ho , zva_ho ) 
     223 
     224         write(numout,*)'ice_dyn_adv_umx 5: zudy, zvdx, zu_cat, zv_cat, zcu_box, zcv_box = ',zudy(3,4), zvdx(3,4), zu_cat(3,4,1), zv_cat(3,4,1), zcu_box(3,4), zcv_box(3,4) 
     225 
    211226         ! 
    212227         !                             ! --------------------------------- ! 
     
    227242            CALL adv_umx( zamsk, kn_umx, jt, kt, zdt, zudy , zvdx , zu_cat, zv_cat, zcu_box, zcv_box, & 
    228243               &                                      psv_i, psv_i ) 
     244            write(numout,*)'ice_dyn_adv_umx 6:psv_i, pe_s, jt = ',psv_i(3,4,1), ' ', pe_i (3,4,1,1), ' ', jt 
    229245            !== Ice heat content ==! 
    230246            DO jk = 1, nlay_i 
     
    254270            CALL adv_umx( zamsk, kn_umx, jt, kt, zdt, zudy , zvdx , zua_ho , zva_ho , zcu_box, zcv_box, & 
    255271               &                                      zhvar, psv_i, zua_ups, zva_ups ) 
     272            write(numout,*)'ice_dyn_adv_umx 7:psv_i, pe_s, jt = ',psv_i(3,4,1), ' ', pe_i (3,4,1,1), ' ', jt 
    256273            !== Ice heat content ==! 
    257274            DO jk = 1, nlay_i 
     
    296313            CALL adv_umx( zamsk, kn_umx, jt, kt, zdt, zudy , zvdx , zuv_ho , zvv_ho , zcu_box, zcv_box, & 
    297314               &                                      zhvar, psv_i, zuv_ups, zvv_ups ) 
     315 
     316            write(numout,*)'ice_dyn_adv_umx 8:psv_i, pe_s, jt = ',psv_i(3,4,1), ' ', pe_i (3,4,1,1), ' ', jt 
    298317            !== Ice heat content ==! 
    299318            DO jk = 1, nlay_i 
     
    360379         !    (because advected fields are not perfectly bounded and tiny negative values can occur, e.g. -1.e-20) 
    361380         CALL ice_var_zapneg( zdt, pato_i, pv_i, pv_s, psv_i, poa_i, pa_i, pa_ip, pv_ip, plh_ip, pe_s, pe_i ) 
     381 
     382         write(numout,*)'ice_dyn_adv_umx 9:psv_i, pe_s, jt = ',psv_i(3,4,1), ' ', pe_i (3,4,1,1), ' ', jt 
    362383         ! 
    363384         ! Make sure ice thickness is not too big 
     
    365386         CALL Hbig( zdt, zhi_max, zhs_max, zhip_max, pv_i, pv_s, psv_i, poa_i, pa_i, pa_ip, pv_ip, pe_s, pe_i ) 
    366387 
    367       END DO 
     388         write(numout,*)'ice_dyn_adv_umx 10:psv_i, pe_s, jt = ',psv_i(3,4,1), ' ', pe_i (3,4,1,1), ' ', jt 
     389 
     390      END DO 
     391 
     392      write(numout,*)'ice_dyn_adv_umx 11:psv_i, pe_s = ',psv_i(3,4,1), ' ', pe_i (3,4,1,1) 
    368393      ! 
    369394   END SUBROUTINE ice_dyn_adv_umx 
  • NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/icedyn_rhg.F90

    r11081 r12475  
    6363      ! controls 
    6464      IF( ln_timing    )   CALL timing_start('icedyn_rhg')                                                             ! timing 
     65      write(numout,*)'ice_dyn_rhg 1: u_ice = ',u_ice(3,4) 
    6566      IF( ln_icediachk )   CALL ice_cons_hsm(0, 'icedyn_rhg', rdiag_v, rdiag_s, rdiag_t, rdiag_fv, rdiag_fs, rdiag_ft) ! conservation 
    6667      ! 
     
    7778         END DO 
    7879      ENDIF 
     80      write(numout,*)'ice_dyn_rhg 2: u_ice = ',u_ice(3,4) 
    7981      ! 
    8082      !--------------! 
     
    8688         !                             !------------------------! 
    8789         CALL ice_dyn_rhg_evp( kt, stress1_i, stress2_i, stress12_i, shear_i, divu_i, delta_i ) 
     90         write(numout,*)'ice_dyn_rhg 3: u_ice = ',u_ice(3,4) 
    8891         !          
    8992      END SELECT 
     
    9194      IF( lrst_ice ) THEN                       !* write EVP fields in the restart file 
    9295         IF( ln_rhg_EVP )   CALL rhg_evp_rst( 'WRITE', kt ) 
     96         write(numout,*)'ice_dyn_rhg 4: u_ice = ',u_ice(3,4) 
    9397      ENDIF 
    9498      ! 
    9599      ! controls 
    96100      IF( ln_icediachk )   CALL ice_cons_hsm(1, 'icedyn_rhg', rdiag_v, rdiag_s, rdiag_t, rdiag_fv, rdiag_fs, rdiag_ft) ! conservation 
     101      write(numout,*)'ice_dyn_rhg 5: u_ice = ',u_ice(3,4) 
    97102      IF( ln_ctl       )   CALL ice_prt3D   ('icedyn_rhg')                                                             ! prints 
    98       IF( ln_timing    )   CALL timing_stop ('icedyn_rhg')                                                             ! timing 
     103      write(numout,*)'ice_dyn_rhg 6: u_ice = ',u_ice(3,4) 
     104      IF( ln_timing    )   CALL timing_stop ('icedyn_rhg')   
     105      write(numout,*)'ice_dyn_rhg 7: u_ice = ',u_ice(3,4)                                                           ! timing 
    99106      ! 
    100107   END SUBROUTINE ice_dyn_rhg 
  • NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/icedyn_rhg_evp.F90

    r11081 r12475  
    185185 
    186186      IF( kt == nit000 .AND. lwp )   WRITE(numout,*) '-- ice_dyn_rhg_evp: EVP sea-ice rheology' 
     187 
     188      write(numout,*)'ice_dyn_rhg_evp 1: u_ice = ',u_ice(3,4) 
    187189      ! 
    188190!!gm for Clem:  OPTIMIZATION:  I think zfmask can be computed one for all at the initialization.... 
     
    225227      CALL lbc_lnk( 'icedyn_rhg_evp', zfmask, 'F', 1._wp ) 
    226228 
     229      write(numout,*)'ice_dyn_rhg_evp 3: u_ice = ',u_ice(3,4) 
     230 
    227231      !------------------------------------------------------------------------------! 
    228232      ! 1) define some variables and initialize arrays 
     
    267271      ELSE                                               ;   zkt = 0._wp 
    268272      ENDIF 
     273 
     274      write(numout,*)'ice_dyn_rhg_evp 3: u_ice = ',u_ice(3,4) 
    269275      ! 
    270276      !------------------------------------------------------------------------------! 
     
    328334      END DO 
    329335      CALL lbc_lnk_multi( 'icedyn_rhg_evp', zmf, 'T', 1., zdt_m, 'T', 1. ) 
     336 
     337      write(numout,*)'ice_dyn_rhg_evp 3: u_ice = ',u_ice(3,4) 
    330338      ! 
    331339      !                                  !== Landfast ice parameterization ==! 
     
    368376      IF( iom_use('tau_icebfr') )   CALL iom_put( 'tau_icebfr', tau_icebfr(:,:) ) 
    369377 
     378      write(numout,*)'ice_dyn_rhg_evp 3: u_ice = ',u_ice(3,4) 
     379 
    370380      !------------------------------------------------------------------------------! 
    371381      ! 3) Solution of the momentum equation, iterative procedure 
     
    397407         CALL lbc_lnk( 'icedyn_rhg_evp', zds, 'F', 1. ) 
    398408 
     409         write(numout,*)'ice_dyn_rhg_evp 3: u_ice = ',u_ice(3,4) 
     410 
    399411         DO jj = 2, jpj    ! loop to jpi,jpj to avoid making a communication for zs1,zs2,zs12 
    400412            DO ji = 2, jpi ! no vector loop 
     
    440452         END DO 
    441453         CALL lbc_lnk( 'icedyn_rhg_evp', zp_delt, 'T', 1. ) 
     454 
     455         write(numout,*)'ice_dyn_rhg_evp 3: u_ice = ',u_ice(3,4) 
    442456 
    443457         DO jj = 1, jpjm1 
     
    488502            END DO 
    489503         END DO 
     504         write(numout,*)'ice_dyn_rhg_evp 3: u_ice = ',u_ice(3,4) 
     505 
    490506         ! 
    491507         ! --- Computation of ice velocity --- ! 
     
    558574                     &    ( zmf(ji  ,jj) * ( e1v(ji  ,jj) * v_ice(ji  ,jj) + e1v(ji  ,jj-1) * v_ice(ji  ,jj-1) )  & 
    559575                     &    + zmf(ji+1,jj) * ( e1v(ji+1,jj) * v_ice(ji+1,jj) + e1v(ji+1,jj-1) * v_ice(ji+1,jj-1) ) ) 
     576 
     577                  IF (ji == 3 .AND. jj == 4) THEN 
     578                    write(numout,*)'ice_dyn_rhg_evp 6: u_ice(ji,jj), u_oce(ji,jj), v_iceU(ji,jj), v_oceU(ji,jj) = ',u_ice(ji,jj), u_oce(ji,jj), v_iceU(ji,jj), v_oceU(ji,jj) 
     579                    write(numout,*)'ice_dyn_rhg_evp 6: zTauO, zTauU_ib(ji,jj), zvel =',zTauO, zTauU_ib(ji,jj), zvel 
     580                    write(numout,*)'ice_dyn_rhg_evp 6: r1_e1u(ji,jj), zmf(ji,jj), e1v(ji,jj) = ',r1_e1u(ji,jj), zmf(ji,jj), e1v(ji,jj) 
     581                  ENDIF 
     582 
    560583                  ! 
    561584                  !                 !--- Sum of external forces (explicit solution) = F + tau_ia + Coriolis + spg + tau_io 
     
    580603                     &            ) * zmaskU(ji,jj) 
    581604                  ENDIF 
     605 
     606                  IF (ji == 3 .AND. jj == 4) THEN 
     607                    write(numout,*)'ice_dyn_rhg_evp 7: u_ice(ji,jj), zmU_t(ji,jj), zTauE, zTauO = ',u_ice(ji,jj), zmU_t(ji,jj), zTauE, zTauO 
     608                    write(numout,*)'ice_dyn_rhg_evp 7: zepsi, zTauO, zTauB, rswitch = ', zepsi, zTauO, zTauB, rswitch 
     609                    write(numout,*)'ice_dyn_rhg_evp 7: zdtevp, rn_lfrelax, zswitchU(ji,jj), zmaskU(ji,jj) = ', zdtevp, rn_lfrelax, zswitchU(ji,jj), zmaskU(ji,jj) 
     610                  ENDIF 
     611 
    582612               END DO 
    583613            END DO 
     
    608638                     &    ( zmf(ji  ,jj) * ( e1v(ji  ,jj) * v_ice(ji  ,jj) + e1v(ji  ,jj-1) * v_ice(ji  ,jj-1) )  & 
    609639                     &    + zmf(ji+1,jj) * ( e1v(ji+1,jj) * v_ice(ji+1,jj) + e1v(ji+1,jj-1) * v_ice(ji+1,jj-1) ) ) 
     640 
     641                  IF (ji == 3 .AND. jj == 4) THEN 
     642                    write(numout,*)'ice_dyn_rhg_evp 8: u_ice(ji,jj), u_oce(ji,jj), v_iceU(ji,jj), v_oceU(ji,jj) = ',u_ice(ji,jj), u_oce(ji,jj), v_iceU(ji,jj), v_oceU(ji,jj) 
     643                    write(numout,*)'ice_dyn_rhg_evp 8: zTauO, zTauU_ib(ji,jj), zvel =',zTauO, zTauU_ib(ji,jj), zvel 
     644                    write(numout,*)'ice_dyn_rhg_evp 8: r1_e1u(ji,jj), zmf(ji,jj), e1v(ji,jj) = ',r1_e1u(ji,jj), zmf(ji,jj), e1v(ji,jj) 
     645                  ENDIF 
     646 
    610647                  ! 
    611648                  !                 !--- Sum of external forces (explicit solution) = F + tau_ia + Coriolis + spg + tau_io 
     
    630667                     &            ) * zmaskU(ji,jj) 
    631668                  ENDIF 
     669 
     670                  IF (ji == 3 .AND. jj == 4) THEN 
     671                    write(numout,*)'ice_dyn_rhg_evp 9: u_ice(ji,jj), zmU_t(ji,jj), zTauE, zTauO = ',u_ice(ji,jj), zmU_t(ji,jj), zTauE, zTauO 
     672                    write(numout,*)'ice_dyn_rhg_evp 9: zepsi, zTauO, zTauB, rswitch = ', zepsi, zTauO, zTauB, rswitch 
     673                    write(numout,*)'ice_dyn_rhg_evp 9: zdtevp, rn_lfrelax, zswitchU(ji,jj), zmaskU(ji,jj) = ', zdtevp, rn_lfrelax, zswitchU(ji,jj), zmaskU(ji,jj) 
     674                  ENDIF 
    632675               END DO 
    633676            END DO 
    634677            CALL lbc_lnk( 'icedyn_rhg_evp', u_ice, 'U', -1. ) 
     678 
     679            write(numout,*)'ice_dyn_rhg_evp 10: u_ice = ',u_ice(3,4) 
    635680            ! 
    636681#if defined key_agrif 
     
    706751         CALL bdy_ice_dyn( 'V' ) 
    707752      ENDIF 
     753      write(numout,*)'ice_dyn_rhg_evp 5: u_ice = ',u_ice(3,4) 
    708754      ! 
    709755      !------------------------------------------------------------------------------! 
     
    757803      pstress2_i (:,:) = zs2 (:,:) 
    758804      pstress12_i(:,:) = zs12(:,:) 
     805      write(numout,*)'ice_dyn_rhg_evp 4: u_ice = ',u_ice(3,4) 
    759806      ! 
    760807 
     
    853900            END DO 
    854901         END DO 
     902 
     903         write(numout,*)'ice_dyn_rhg_evp 11: u_ice = ',u_ice(3,4) 
    855904          
    856905         CALL lbc_lnk_multi( 'icedyn_rhg_evp', zdiag_sig1   , 'T',  1., zdiag_sig2   , 'T',  1.,   & 
     
    887936 
    888937      ENDIF 
     938 
     939      write(numout,*)'ice_dyn_rhg_evp 2: u_ice = ',u_ice(3,4) 
    889940      ! 
    890941   END SUBROUTINE ice_dyn_rhg_evp 
  • NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/icestp.F90

    r12449 r12475  
    188188         kt_ice = kt                              ! -- Ice model time step 
    189189 
    190          write(numout,*)'icestp 0: e_s = ',e_s (42,26,1,1) 
     190         write(numout,*)'icestp 0: sv_i, e_i, u_ice, e2u, u_oce = ',sv_i(3,4,1), ' ', e_i (3,4,1,1), ' ', u_ice(3,4), ' ', e2u(3,4), ' ', u_oce(3,4) 
    191191         ! 
    192192         u_oce(:,:) = ssu_m(:,:)                  ! -- mean surface ocean current 
     
    215215                                        CALL ice_sbc_tau( kt, ksbc, utau_ice, vtau_ice )       
    216216 
    217          write(numout,*)'icestp 1: e_s = ',e_s (42,26,1,1) 
     217         write(numout,*)'icestp 1: sv_i, e_i, u_ice, e2u, u_oce = ',sv_i(3,4,1), ' ', e_i (3,4,1,1), ' ', u_ice(3,4), ' ', e2u(3,4), ' ', u_oce(3,4) 
    218218     
    219219         !-------------------------------------! 
     
    226226            &                           CALL ice_dyn( kt )            ! -- Ice dynamics 
    227227 
    228          write(numout,*)'icestp 2: e_s = ',e_s (42,26,1,1) 
     228         write(numout,*)'icestp 2: sv_i, e_i, u_ice, e2u, u_oce = ',sv_i(3,4,1), ' ', e_i (3,4,1,1), ' ', u_ice(3,4), ' ', e2u(3,4), ' ', u_oce(3,4) 
    229229         ! 
    230230         !                          !==  lateral boundary conditions  ==! 
    231231         IF( ln_icethd .AND. ln_bdy )   CALL bdy_ice( kt )            ! -- bdy ice thermo 
    232232 
    233          write(numout,*)'icestp 3: e_s = ',e_s (42,26,1,1) 
     233         write(numout,*)'icestp 3: sv_i, e_i, u_ice, e2u, u_oce = ',sv_i(3,4,1), ' ', e_i (3,4,1,1), ' ', u_ice(3,4), ' ', e2u(3,4), ' ', u_oce(3,4) 
    234234         ! 
    235235         !                          !==  previous lead fraction and ice volume for flux calculations 
    236236                                        CALL ice_var_glo2eqv          ! h_i and h_s for ice albedo calculation 
    237237 
    238          write(numout,*)'icestp 4: t_s, e_s = ',t_s(ji,jj,jk,jl) - rt0, ' ', e_s (42,26,1,1) 
     238         write(numout,*)'icestp 4: sv_i, e_i, u_ice, e2u, u_oce = ',sv_i(3,4,1), ' ', e_i (3,4,1,1), ' ', u_ice(3,4), ' ', e2u(3,4), ' ', u_oce(3,4) 
    239239 
    240240                                        CALL ice_var_agg(1)           ! at_i for coupling  
    241241 
    242          write(numout,*)'icestp 5: t_s, e_s = ',t_s(ji,jj,jk,jl) - rt0, ' ', e_s (42,26,1,1) 
     242         write(numout,*)'icestp 5: sv_i, e_i, u_ice, e2u, u_oce = ',sv_i(3,4,1), ' ', e_i (3,4,1,1), ' ', u_ice(3,4), ' ', e2u(3,4), ' ', u_oce(3,4) 
    243243 
    244244                                        CALL store_fields             ! Store now ice values 
    245245 
    246          write(numout,*)'icestp 6: t_s, e_s = ',t_s(ji,jj,jk,jl) - rt0, ' ', e_s (42,26,1,1) 
     246         write(numout,*)'icestp 6: sv_i, e_i, u_ice, e2u, u_oce = ',sv_i(3,4,1), ' ', e_i (3,4,1,1), ' ', u_ice(3,4), ' ', e2u(3,4), ' ', u_oce(3,4) 
    247247         ! 
    248248         !------------------------------------------------------! 
     
    260260         !------------------------------------------------------! 
    261261                                        CALL ice_sbc_flx( kt, ksbc ) 
    262          write(numout,*)'icestp 7: t_s, e_s = ',t_s(ji,jj,jk,jl) - rt0, ' ', e_s (42,26,1,1) 
     262         write(numout,*)'icestp 7: sv_i, e_s = ',sv_i(3,4,1), ' ', e_s (3,4,1,1) 
    263263         !----------------------------! 
    264264         ! --- ice thermodynamics --- ! 
    265265         !----------------------------! 
    266266         IF( ln_icethd )                CALL ice_thd( kt )            ! -- Ice thermodynamics    
    267          write(numout,*)'icestp 8: t_s, e_s = ',t_s(ji,jj,jk,jl) - rt0, ' ', e_s (42,26,1,1)    
     267         write(numout,*)'icestp 8: sv_i, e_s = ',sv_i(3,4,1), ' ', e_s (3,4,1,1)    
    268268         ! 
    269269                                        CALL ice_cor( kt , 2 )        ! -- Corrections 
    270          write(numout,*)'icestp 9: t_s, e_s = ',t_s(ji,jj,jk,jl) - rt0, ' ', e_s (42,26,1,1) 
     270         write(numout,*)'icestp 9: sv_i, e_s = ',sv_i(3,4,1), ' ', e_s (3,4,1,1) 
    271271         ! 
    272272                                        CALL ice_var_glo2eqv          ! necessary calls (at least for coupling) 
     
    490490      wfx_snw(:,:) = 0._wp   ;   wfx_ice(:,:) = 0._wp 
    491491      wfx_sni(:,:) = 0._wp   ;   wfx_opw(:,:) = 0._wp 
     492      wfx_snp(:,:) = 0._wp 
    492493      wfx_bog(:,:) = 0._wp   ;   wfx_dyn(:,:) = 0._wp 
    493494      wfx_bom(:,:) = 0._wp   ;   wfx_sum(:,:) = 0._wp 
  • NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/icethd.F90

    r12449 r12475  
    442442         CALL tab_2d_1d( npti, nptidx(1:npti), wfx_sum_1d (1:npti), wfx_sum          ) 
    443443         CALL tab_2d_1d( npti, nptidx(1:npti), wfx_sni_1d (1:npti), wfx_sni          ) 
     444         CALL tab_2d_1d( npti, nptidx(1:npti), wfx_snp_1d (1:npti), wfx_snp          ) 
    444445         CALL tab_2d_1d( npti, nptidx(1:npti), wfx_res_1d (1:npti), wfx_res          ) 
    445446         CALL tab_2d_1d( npti, nptidx(1:npti), wfx_spr_1d (1:npti), wfx_spr          ) 
     
    536537         CALL tab_1d_2d( npti, nptidx(1:npti), wfx_sum_1d (1:npti), wfx_sum        ) 
    537538         CALL tab_1d_2d( npti, nptidx(1:npti), wfx_sni_1d (1:npti), wfx_sni        ) 
     539         CALL tab_1d_2d( npti, nptidx(1:npti), wfx_snp_1d (1:npti), wfx_snp        ) 
    538540         CALL tab_1d_2d( npti, nptidx(1:npti), wfx_res_1d (1:npti), wfx_res        ) 
    539541         CALL tab_1d_2d( npti, nptidx(1:npti), wfx_spr_1d (1:npti), wfx_spr        ) 
  • NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/icethd_dh.F90

    r12449 r12475  
    296296               ! Snow melt trapped in ponds contrbutes to ice thickness and snow to ice conversion 
    297297               h_i_1d  (ji)    = MAX( 0._wp , h_i_1d(ji) - zdeltah_pnd(ji,jk) * rhos * r1_rhoi )  
    298                wfx_sni_1d(ji)     = wfx_sni_1d(ji)     - a_i_1d(ji) * zdeltah_pnd(ji,jk) * rhos * r1_rdtice 
    299                wfx_snw_sni_1d(ji) = wfx_snw_sni_1d(ji) + a_i_1d(ji) * zdeltah_pnd(ji,jk) * rhos * r1_rdtice 
     298               wfx_snp_1d(ji)     = wfx_snp_1d(ji)     + a_i_1d(ji) * zdeltah_pnd(ji,jk) * rhos * r1_rdtice  ! As zdeltah_pnd is negative this comes out as negative 
    300299               ! 
    301300            ENDIF 
  • NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/icethd_zdf_bl99.F90

    r12374 r12475  
    102102      ! 
    103103      REAL(wp), DIMENSION(jpij) ::   isnow        ! fraction of sea ice that is snow covered (for thermodynamic use only)  
    104       LOGICAL,  DIMENSION(jpij) ::   l_snow       ! is snow above the threshold 
    105104      REAL(wp), DIMENSION(jpij) ::   ztsub        ! surface temperature at previous iteration 
    106105      REAL(wp), DIMENSION(jpij) ::   zh_i, z1_h_i ! ice layer thickness 
     
    157156         ! If the snow thickness drops below zhs_min then reduce the snow fraction instead  
    158157         IF( h_s_1d(ji) < zhs_min ) THEN  
    159            ! isnow(ji) = h_s_1d(ji) / zhs_min  
    160            isnow(ji) = 0.0_wp 
     158           isnow(ji) = h_s_1d(ji) / zhs_min  
    161159           zh_s(ji) = zhs_min * r1_nlay_s  
    162            l_snow(ji) = .false. 
    163160         ELSE  
    164161           isnow(ji) = 1.0_wp  
    165162           zh_s(ji) = h_s_1d(ji) * r1_nlay_s  
    166            l_snow(ji) = .true. 
    167163         END IF  
    168164 
     
    322318               IF ( .NOT. l_T_converged(ji) ) & 
    323319                  zkappa_s(ji,jk) = zghe(ji) * rn_cnd_s * z1_h_s(ji) 
    324                IF ( .NOT. l_snow(ji) ) zkappa_s(ji,jk) = 0.0_wp 
    325320            END DO 
    326321         END DO 
     
    370365         DO jk = 1, nlay_s 
    371366            DO ji = 1, npti 
    372                IF (l_snow(ji)) then 
    373                   zeta_s(ji,jk) = rdt_ice * r1_rhos * r1_rcpi * z1_h_s(ji) 
    374                ELSE 
    375                   zeta_s(ji,jk) = 0.0_wp 
    376                END IF 
     367               zeta_s(ji,jk) = rdt_ice * r1_rhos * r1_rcpi * z1_h_s(ji) 
    377368            END DO 
    378369         END DO 
     
    722713            DO ji = 1, npti 
    723714               !                               !---------------------! 
    724                IF( l_snow(ji) ) THEN   !  snow-covered cells ! 
     715               IF( h_s_1d(ji) > 0._wp ) THEN   !  snow-covered cells ! 
    725716                  !                            !---------------------! 
    726717                  ! snow interior terms (bottom equation has the same form as the others) 
     
    779770                     zindterm(ji,jm_min(ji))   = ztiold(ji,1) + zeta_i(ji,1) *  & 
    780771                        &                                     ( zradab_i(ji,1) + zkappa_i(ji,1) * t_bo_1d(ji) + qcn_ice_1d(ji) ) 
    781                   ENDIF 
     772                  ENDIF                
    782773                   
    783774               ENDIF 
     775 
     776                  IF (to_print(ji) == 10) THEN 
     777                     write(numout,*)'after minmax: h_s_1d(ji), jm_min(ji), nlay_s = ',h_s_1d(ji), ' ', jm_min(ji), ' ', nlay_s 
     778                  ENDIF    
     779 
    784780               ! 
    785781               zindtbis(ji,jm_min(ji)) = zindterm(ji,jm_min(ji)) 
     
    802798            DO jk = jm_mint+1, jm_maxt 
    803799               DO ji = 1, npti 
     800 
     801  
    804802                  jm = MIN(MAX(jm_min(ji)+1,jk),jm_max(ji)) 
     803                  IF (to_print(ji) == 10) THEN 
     804                     write(numout,*)'before calc zdiagbis: jk, jm, zdiagbis(ji,jm-1) = ',jk, ' ', jm, ' ', zdiagbis(ji,jm-1) 
     805                     write(numout,*)'before calc zdiagbis:ztrid(ji,jm,2), ztrid(ji,jm,1), ztrid(ji,jm-1,3) = ', ztrid(ji,jm,2), ' ', ztrid(ji,jm,1), ' ', ztrid(ji,jm-1,3) 
     806                  ENDIF 
     807 
    805808                  zdiagbis(ji,jm) = ztrid   (ji,jm,2) - ztrid(ji,jm,1) * ztrid   (ji,jm-1,3) / zdiagbis(ji,jm-1) 
    806809                  zindtbis(ji,jm) = zindterm(ji,jm)   - ztrid(ji,jm,1) * zindtbis(ji,jm-1)   / zdiagbis(ji,jm-1) 
     
    838841                write(numout,*)'icethd_zdf_bl99 2: t_s_1d(ji,1), e_s_1d(ji,1) = ',t_s_1d(ji,1), ' ', e_s_1d(ji,1) 
    839842                write(numout,*)'icethd_zdf_bl99 2: zindtbis(ji,nlay_s+1), ztrid(ji,nlay_s+1,3), t_i_1d(ji,1) = ', zindtbis(ji,nlay_s+1), ' ', ztrid(ji,nlay_s+1,3), ' ', t_i_1d(ji,1) 
     843                write(numout,*)'icethd_zdf_bl99 2: zdiagbis(ji,nlay_s+1), h_s_1d(ji) = ',zdiagbis(ji,nlay_s+1), ' ', h_s_1d(ji) 
    840844              END IF 
    841845            END DO 
  • NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/iceupdate.F90

    r12449 r12475  
    154154               &           + wfx_opw(ji,jj) + wfx_dyn(ji,jj) + wfx_res(ji,jj) + wfx_lam(ji,jj) + wfx_pnd(ji,jj) 
    155155 
     156 
     157 
    156158            ! add the snow melt water to snow mass flux to the ocean 
    157159            wfx_snw(ji,jj) = wfx_snw_sni(ji,jj) + wfx_snw_dyn(ji,jj) + wfx_snw_sum(ji,jj) 
     
    161163            emp(ji,jj)    = emp_oce(ji,jj) - wfx_ice(ji,jj) - wfx_snw(ji,jj) - wfx_err_sub(ji,jj)   ! mass flux + F/M mass flux (always ice/ocean mass exchange) 
    162164 
     165 
     166             
    163167 
    164168            ! Salt flux at the ocean surface       
     
    166170            sfx(ji,jj) = sfx_bog(ji,jj) + sfx_bom(ji,jj) + sfx_sum(ji,jj) + sfx_sni(ji,jj) + sfx_opw(ji,jj)   & 
    167171               &       + sfx_res(ji,jj) + sfx_dyn(ji,jj) + sfx_bri(ji,jj) + sfx_sub(ji,jj) + sfx_lam(ji,jj) + sfx_pnd(ji,jj) 
     172 
     173 
    168174             
    169175            ! Mass of snow and ice per unit area    
     
    177183         END DO 
    178184      END DO 
     185 
     186            write(numout,*)'iceupdate: MEAN(wfx_pnd) = ',SUM(wfx_pnd)/SIZE(wfx_pnd) 
     187            write(numout,*)'iceupdate: MEAN(wfx_sum) = ',SUM(wfx_sum)/SIZE(wfx_sum) 
     188            write(numout,*)'iceupdate: MEAN(wfx_bom) = ',SUM(wfx_bom)/SIZE(wfx_bom) 
     189            write(numout,*)'iceupdate: MEAN(wfx_ice) = ',SUM(wfx_ice)/SIZE(wfx_ice) 
     190            write(numout,*)'iceupdate: MEAN(wfx_sni) = ',SUM(wfx_sni)/SIZE(wfx_sni) 
     191            write(numout,*)'iceupdate: MEAN(wfx_snp) = ',SUM(wfx_snp)/SIZE(wfx_snp) 
     192            write(numout,*)'iceupdate: MEAN(fmmflx) = ',SUM(fmmflx)/SIZE(fmmflx) 
     193            write(numout,*)'iceupdate: MEAN(wfx_snw) = ',SUM(wfx_snw)/SIZE(wfx_snw) 
     194            write(numout,*)'iceupdate: MEAN(emp) = ',SUM(emp)/SIZE(emp) 
     195            write(numout,*)'iceupdate: MEAN(emp_oce) = ',SUM(emp_oce)/SIZE(emp_oce) 
     196            write(numout,*)'iceupdate: MEAN(sfx) = ',SUM(sfx)/SIZE(sfx) 
     197            write(numout,*)'iceupdate: MEAN(sfx_sni) = ',SUM(sfx_sni)/SIZE(sfx_sni) 
     198            write(numout,*)'iceupdate: MEAN(sfx_pnd) = ',SUM(sfx_pnd)/SIZE(sfx_pnd) 
     199            write(numout,*)'iceupdate: MEAN(sfx_sum) = ',SUM(sfx_sum)/SIZE(sfx_sum) 
     200            write(numout,*)'iceupdate: MEAN(sfx_bom) = ',SUM(sfx_sum)/SIZE(sfx_bom) 
     201            write(numout,*)'iceupdate: MAXVAL(wfx_pnd) = ',MAXVAL(wfx_pnd) 
     202            write(numout,*)'iceupdate: MINVAL(wfx_pnd) = ',MINVAL(wfx_pnd) 
     203            write(numout,*)'iceupdate: MAXVAL(wfx_bom) = ',MAXVAL(wfx_bom) 
     204            write(numout,*)'iceupdate: MINVAL(wfx_bom) = ',MINVAL(wfx_bom) 
     205            write(numout,*)'iceupdate: MAXVAL(wfx_sni) = ',MAXVAL(wfx_sni) 
     206            write(numout,*)'iceupdate: MINVAL(wfx_sni) = ',MINVAL(wfx_sni) 
     207            write(numout,*)'iceupdate: MAXVAL(sfx_pnd) = ',MAXVAL(sfx_pnd) 
     208            write(numout,*)'iceupdate: MINVAL(sfx_pnd) = ',MINVAL(sfx_pnd) 
     209            write(numout,*)'iceupdate: MAXVAL(sfx_bom) = ',MAXVAL(sfx_bom) 
     210            write(numout,*)'iceupdate: MINVAL(sfx_bom) = ',MINVAL(sfx_bom) 
     211 
    179212 
    180213      ! Storing the transmitted variables 
  • NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/icewri.F90

    r12439 r12475  
    172172      !------------------ 
    173173      ! trends 
    174       IF( iom_use('dmithd') )   CALL iom_put( "dmithd", - wfx_bog - wfx_bom - wfx_sum - wfx_sni - wfx_opw - wfx_lam - wfx_res - wfx_pnd ) ! Sea-ice mass change from thermodynamics 
     174      IF( iom_use('dmithd') )   CALL iom_put( "dmithd", - wfx_bog - wfx_bom - wfx_sum - wfx_sni - wfx_opw - wfx_lam - wfx_res - wfx_pnd - wfx_snp ) ! Sea-ice mass change from thermodynamics 
    175175      IF( iom_use('dmidyn') )   CALL iom_put( "dmidyn", - wfx_dyn + rhoi * diag_trp_vi      )   ! Sea-ice mass change from dynamics(kg/m2/s) 
    176176      IF( iom_use('dmiopw') )   CALL iom_put( "dmiopw", - wfx_opw                           )   ! Sea-ice mass change through growth in open water 
    177177      IF( iom_use('dmibog') )   CALL iom_put( "dmibog", - wfx_bog                           )   ! Sea-ice mass change through basal growth 
    178178      IF( iom_use('dmisni') )   CALL iom_put( "dmisni", - wfx_sni                           )   ! Sea-ice mass change through snow-to-ice conversion 
     179      IF( iom_use('dmisnp') )   CALL iom_put( "dmisnp", - wfx_snp                           )   ! Sea-ice mass change through snow melt into ponds 
    179180      IF( iom_use('dmisum') )   CALL iom_put( "dmisum", - wfx_sum                           )   ! Sea-ice mass change through surface melting 
    180181      IF( iom_use('dmipnd') )   CALL iom_put( "dmipnd", - wfx_pnd                           )   ! Sea-ice mass change through meltponds overflowing to ocean 
  • NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/OCE/IOM/restart.F90

    r10888 r12475  
    281281         neuler = 0 
    282282      ENDIF 
     283 
     284      write(numout,*)'rst_read 2: ub(3,4,1) = ', ub(3,4,1) 
     285 
    283286      ! 
    284287      CALL iom_get( numror, jpdom_autoglo, 'un'     , un, ldxios = lrxios )   ! now    fields 
     
    306309         ! 
    307310      ENDIF 
     311      write(numout,*)'rst_read 3: ub(3,4,1) = ', ub(3,4,1) 
    308312      ! 
    309313   END SUBROUTINE rst_read 
  • NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/OCE/SBC/sbcmod.F90

    r10888 r12475  
    407407      !                                            !        forcing field computation         ! 
    408408      !                                            ! ---------------------------------------- ! 
     409 
     410      write(numout,*)'sbcmod 1: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 
     411 
    409412      ! 
    410413      ll_sas = nn_components == jp_iam_sas               ! component flags 
     
    413416      IF( .NOT.ll_sas )   CALL sbc_ssm ( kt )            ! mean ocean sea surface variables (sst_m, sss_m, ssu_m, ssv_m) 
    414417      IF( ln_wave     )   CALL sbc_wave( kt )            ! surface waves 
     418 
     419      write(numout,*)'sbcmod 2: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 
    415420 
    416421      ! 
     
    436441      !                                            !==  Misc. Options  ==! 
    437442      ! 
     443      write(numout,*)'sbcmod 3: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 
     444 
    438445      SELECT CASE( nn_ice )                                       ! Update heat and freshwater fluxes over sea-ice areas 
    439446      CASE(  1 )   ;         CALL sbc_ice_if   ( kt )             ! Ice-cover climatology ("Ice-if" model) 
  • NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/OCE/SBC/sbcssm.F90

    r10888 r12475  
    5959      !!--------------------------------------------------------------------- 
    6060      ! 
     61      write(numout,*)'sbcssm 1: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 
    6162      !                                        !* surface T-, U-, V- ocean level variables (T, S, depth, velocity) 
    6263      DO jj = 1, jpj 
     
    8384         ! 
    8485         frq_m(:,:) = fraqsr_1lev(:,:) 
     86         write(numout,*)'sbcssm 2: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 
    8587         ! 
    8688      ELSE 
     
    106108            ! 
    107109            frq_m(:,:) = zcoef * fraqsr_1lev(:,:) 
     110            write(numout,*)'sbcssm 3: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 
    108111            !                                             ! ---------------------------------------- ! 
    109112         ELSEIF( MOD( kt - 2 , nn_fsbc ) == 0 ) THEN      !   Initialisation: New mean computation   ! 
     
    116119            e3t_m(:,:) = 0._wp 
    117120            frq_m(:,:) = 0._wp 
     121            write(numout,*)'sbcssm 4: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 
    118122         ENDIF 
    119123         !                                                ! ---------------------------------------- ! 
     
    121125         !                                                ! ---------------------------------------- ! 
    122126         ssu_m(:,:) = ssu_m(:,:) + ub(:,:,1) 
     127         write(numout,*)'sbcssm 5: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 
    123128         ssv_m(:,:) = ssv_m(:,:) + vb(:,:,1) 
    124129         IF( l_useCT )  THEN     ;   sst_m(:,:) = sst_m(:,:) + eos_pt_from_ct( zts(:,:,jp_tem), zts(:,:,jp_sal) ) 
     
    146151            e3t_m(:,:) = e3t_m(:,:) * zcoef     ! mean vertical scale factor [m] 
    147152            frq_m(:,:) = frq_m(:,:) * zcoef     ! mean fraction of solar net radiation absorbed in the 1st T level [-] 
     153            write(numout,*)'sbcssm 6: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 
    148154            ! 
    149155         ENDIF 
     
    167173            ! 
    168174            IF( lwxios ) CALL iom_swap(      cxios_context          ) 
     175            write(numout,*)'sbcssm 7: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 
    169176         ENDIF 
    170177         ! 
     
    180187         CALL iom_put( 'frq_m', frq_m ) 
    181188      ENDIF 
     189 
     190      write(numout,*)'sbcssm 8: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 
     191 
    182192      ! 
    183193   END SUBROUTINE sbc_ssm 
     
    195205      !!---------------------------------------------------------------------- 
    196206      ! 
     207 
     208      write(numout,*)'sbcssm 1: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 
     209 
    197210      IF( nn_fsbc == 1 ) THEN 
    198211         ! 
     
    216229            CALL iom_get( numror, jpdom_autoglo, 'ssh_m'  , ssh_m, ldxios = lrxios  )    !   "         "    height      (T-point) 
    217230            CALL iom_get( numror, jpdom_autoglo, 'e3t_m'  , e3t_m, ldxios = lrxios  )    ! 1st level thickness          (T-point) 
     231 
     232            write(numout,*)'sbcssm 2: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 
     233 
    218234            ! fraction of solar net radiation absorbed in 1st T level 
    219235            IF( iom_varid( numror, 'frq_m', ldstop = .FALSE. ) > 0 ) THEN 
     
    233249               e3t_m(:,:) = zcoef * e3t_m(:,:) 
    234250               frq_m(:,:) = zcoef * frq_m(:,:) 
     251               write(numout,*)'sbcssm 3: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 
    235252            ELSE 
    236253               IF(lwp) WRITE(numout,*) '   mean fields read in the ocean restart file' 
     
    238255         ENDIF 
    239256      ENDIF 
     257 
     258      write(numout,*)'sbcssm 4: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 
     259 
    240260      ! 
    241261      IF( .NOT.l_ssm_mean ) THEN   ! default initialisation. needed by iceistate 
     
    253273         ! 
    254274      ENDIF 
     275 
     276      write(numout,*)'sbcssm 5: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 
     277 
    255278      ! 
    256279      IF( .NOT. ln_traqsr )   fraqsr_1lev(:,:) = 1._wp   ! default definition: qsr 100% in the fisrt level  
     
    267290      ENDIF 
    268291 
     292      write(numout,*)'sbcssm 6: ssu_m(3,4), ub(3,4,1) = ', ssu_m(3,4), ' ', ub(3,4,1) 
     293 
    269294   END SUBROUTINE sbc_ssm_init 
    270295 
  • NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/OCE/step.F90

    r10888 r12475  
    9393      ! 
    9494      IF( ln_timing )   CALL timing_start('stp') 
     95      write(numout,*)'stp 1: ub(3,4,1) = ', ub(3,4,1) 
    9596      ! 
    9697      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     
    110111      ! Update external forcing (tides, open boundaries, and surface boundary condition (including sea-ice) 
    111112      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
     113      write(numout,*)'stp 2: ub(3,4,1) = ', ub(3,4,1) 
    112114      IF( ln_tide    )   CALL sbc_tide( kstp )                   ! update tide potential 
     115      write(numout,*)'stp 3: ub(3,4,1) = ', ub(3,4,1) 
    113116      IF( ln_apr_dyn )   CALL sbc_apr ( kstp )                   ! atmospheric pressure (NB: call before bdy_dta which needs ssh_ib)  
     117      write(numout,*)'stp 4: ub(3,4,1) = ', ub(3,4,1) 
    114118      IF( ln_bdy     )   CALL bdy_dta ( kstp, time_offset=+1 )   ! update dynamic & tracer data at open boundaries 
     119 
     120      write(numout,*)'stp 5: ub(3,4,1) = ', ub(3,4,1) 
     121 
    115122                         CALL sbc     ( kstp )                   ! Sea Boundary Condition (including sea-ice) 
    116123 
Note: See TracChangeset for help on using the changeset viewer.