- Timestamp:
- 2012-04-18T12:42:56+02:00 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_r3337_NOCS10_ICB/NEMOGCM/NEMO/OPA_SRC/ICB/icbrun.F90
r3339 r3359 39 39 PRIVATE 40 40 41 PUBLIC icb_stp ! routine called in xxx.F90 module42 PUBLIC icb_end ! routine called in xxx.F90 module41 PUBLIC icb_stp ! routine called in sbcmod.F90 module 42 PUBLIC icb_end ! routine called in nemogcm.F90 module 43 43 44 44 CONTAINS … … 52 52 !! ** Method : - blah blah 53 53 !!---------------------------------------------------------------------- 54 INTEGER, INTENT( in ) :: kt54 INTEGER, INTENT( in ) :: kt 55 55 ! 56 56 INTEGER :: iyr, imon, iday, ihr, imin, isec 57 LOGICAL :: l err, sample_traj, l_budget,l_verbose57 LOGICAL :: ll_sample_traj, ll_budget, ll_verbose 58 58 !!---------------------------------------------------------------------- 59 59 !! start of timestep housekeeping 60 60 61 ktberg = kt61 nktberg = kt 62 62 63 63 ! read calving data … … 65 65 ! 66 66 CALL fld_read ( kt, 1, sf_icb ) 67 p_calving(:,:) = sf_icb(1)%fnow(:,:,1) ! calving in km^3/year (water equivalent)68 p_calving_hflx(:,:) = 0._wp ! NO heat flux for now67 src_calving(:,:) = sf_icb(1)%fnow(:,:,1) ! calving in km^3/year (water equivalent) 68 src_calving_hflx(:,:) = 0._wp ! NO heat flux for now 69 69 ! 70 70 ENDIF … … 87 87 current_yearday = yearday(imon, iday, ihr, imin, isec) 88 88 89 l _verbose = .FALSE.89 ll_verbose = .FALSE. 90 90 IF( nn_verbose_write .GT. 0 .AND. & 91 MOD(kt-1,nn_verbose_write ) == 0 ) l _verbose = nn_verbose_level >= 092 IF( l _verbose ) WRITE(numicb,9100) iyr, imon, iday, ihr, imin, isec, &91 MOD(kt-1,nn_verbose_write ) == 0 ) ll_verbose = nn_verbose_level >= 0 92 IF( ll_verbose ) WRITE(numicb,9100) iyr, imon, iday, ihr, imin, isec, & 93 93 current_year, current_yearday 94 94 9100 FORMAT('y,m,d=',3i5,' h,m,s=',3i5,' yr,yrdy=',i5,f8.3) … … 112 112 113 113 IF( lk_mpp ) THEN 114 114 CALL mpp_send_bergs () ! Send bergs to other PEs 115 115 ELSE 116 116 CALL lbc_send_bergs() ! Deal with any cyclic boundaries in non-mpp case 117 117 ENDIF 118 118 … … 122 122 !! end of timestep housekeeping 123 123 124 sample_traj = .FALSE.125 IF( nn_sample_rate .GT. 0 .AND. MOD(kt-1,nn_sample_rate) == 0 ) sample_traj = .TRUE.126 IF( sample_traj .AND. &124 ll_sample_traj = .FALSE. 125 IF( nn_sample_rate .GT. 0 .AND. MOD(kt-1,nn_sample_rate) == 0 ) ll_sample_traj = .TRUE. 126 IF( ll_sample_traj .AND. & 127 127 ASSOCIATED(first_berg) ) CALL traj_write ( kt ) ! For each berg, record trajectory 128 128 … … 142 142 143 143 ! Diagnose budgets 144 l _budget = .FALSE.145 IF( nn_verbose_write .GT. 0 .AND. MOD(kt-1,nn_verbose_write) == 0 ) l _budget = ln_bergdia146 CALL icb_budget( l _budget )144 ll_budget = .FALSE. 145 IF( nn_verbose_write .GT. 0 .AND. MOD(kt-1,nn_verbose_write) == 0 ) ll_budget = ln_bergdia 146 CALL icb_budget( ll_budget ) 147 147 148 148 IF( MOD(kt,nn_stock) == 0 ) THEN … … 175 175 DEALLOCATE( ssh_e ) 176 176 177 DEALLOCATE( icbfldpts ) 177 DEALLOCATE( nicbfldpts ) 178 DEALLOCATE( nicbflddest ) 179 DEALLOCATE( nicbfldproc ) 178 180 179 181 IF( lk_mpp ) CALL dealloc_buffers() … … 206 208 DEALLOCATE(berg_grid) 207 209 208 DEALLOCATE( initial_width)209 DEALLOCATE( initial_length)210 DEALLOCATE(first_width) 211 DEALLOCATE(first_length) 210 212 211 213 IF (lwp) WRITE(numout,'(a,i6)') 'icebergs: icb_end complete',narea
Note: See TracChangeset
for help on using the changeset viewer.