New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 12475 for NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/OCE/SBC/sbcssm.F90 – NEMO

Ignore:
Timestamp:
2020-02-27T12:27:50+01:00 (4 years 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.
File:
1 edited

Legend:

Unmodified
Added
Removed
  • 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 
Note: See TracChangeset for help on using the changeset viewer.