Changeset 5837 for branches/2014/dev_r4650_UKMO14.4_OBS_GENERAL_VINTERP/NEMOGCM/NEMO/OPA_SRC/C1D/step_c1d.F90
- Timestamp:
- 2015-10-26T15:59:39+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4650_UKMO14.4_OBS_GENERAL_VINTERP/NEMOGCM/NEMO/OPA_SRC/C1D/step_c1d.F90
- Property svn:keywords set to Id
r4313 r5837 32 32 !!---------------------------------------------------------------------- 33 33 !! NEMO/C1D 3.3 , NEMO Consortium (2010) 34 !! $Id : step_c1d.F90 2382 2010-11-13 13:08:12Z gm$34 !! $Id$ 35 35 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 36 36 !!---------------------------------------------------------------------- … … 72 72 ! Ocean physics update (ua, va, ta, sa used as workspace) 73 73 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 74 CALL bn2( tsb, rn2b ) ! before Brunt-Vaisala frequency 75 CALL bn2( tsn, rn2 ) ! now Brunt-Vaisala frequency 74 CALL eos_rab( tsb, rab_b ) ! before local thermal/haline expension ratio at T-points 75 CALL eos_rab( tsn, rab_n ) ! now local thermal/haline expension ratio at T-points 76 CALL bn2( tsb, rab_b, rn2b ) ! before Brunt-Vaisala frequency 77 CALL bn2( tsn, rab_n, rn2 ) ! now Brunt-Vaisala frequency 76 78 ! VERTICAL PHYSICS 77 79 CALL zdf_bfr( kstp ) ! bottom friction … … 115 117 ! Passive Tracer Model 116 118 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 117 CALL trc_stp( kstp )! time-stepping119 CALL trc_stp( kstp ) ! time-stepping 118 120 #endif 119 121 … … 121 123 ! Active tracers (ua, va used as workspace) 122 124 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 123 tsa(:,:,:,:) = 0.e0! set tracer trends to zero125 tsa(:,:,:,:) = 0._wp ! set tracer trends to zero 124 126 125 CALL tra_sbc ( kstp ) ! surface boundary condition 126 IF( ln_traqsr ) CALL tra_qsr ( kstp ) ! penetrative solar radiation qsr 127 IF( ln_tradmp ) CALL tra_dmp ( kstp ) ! internal damping trends- tracers 128 IF( lk_zdfkpp ) CALL tra_kpp ( kstp ) ! KPP non-local tracer fluxes 129 CALL tra_zdf ( kstp ) ! vertical mixing 130 CALL tra_nxt ( kstp ) ! tracer fields at next time step 131 IF( ln_zdfnpc ) CALL tra_npc ( kstp ) ! applied non penetrative convective adjustment on (t,s) 132 CALL eos( tsb, rhd, rhop, gdept_0(:,:,:) ) ! now (swap=before) in situ density for dynhpg module 127 CALL tra_sbc( kstp ) ! surface boundary condition 128 IF( ln_traqsr ) CALL tra_qsr( kstp ) ! penetrative solar radiation qsr 129 IF( ln_tradmp ) CALL tra_dmp( kstp ) ! internal damping trends- tracers 130 IF( lk_zdfkpp ) CALL tra_kpp( kstp ) ! KPP non-local tracer fluxes 131 CALL tra_zdf( kstp ) ! vertical mixing 132 CALL eos( tsn, rhd, rhop, gdept_0(:,:,:) ) ! now potential density for zdfmxl 133 IF( ln_zdfnpc ) CALL tra_npc( kstp ) ! applied non penetrative convective adjustment on (t,s) 134 CALL tra_nxt( kstp ) ! tracer fields at next time step 135 136 133 137 134 138 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 135 139 ! Dynamics (ta, sa used as workspace) 136 140 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 137 ua(:,:,:) = 0.e0! set dynamics trends to zero138 va(:,:,:) = 0.e0141 ua(:,:,:) = 0._wp ! set dynamics trends to zero 142 va(:,:,:) = 0._wp 139 143 140 IF( ln_dyndmp ) CALL dyn_dmp ( kstp )! internal damping trends- momentum141 CALL dyn_cor_c1d( kstp )! vorticity term including Coriolis142 CALL dyn_zdf ( kstp )! vertical diffusion143 CALL dyn_nxt_c1d( kstp )! lateral velocity at next time step144 IF( ln_dyndmp ) CALL dyn_dmp ( kstp ) ! internal damping trends- momentum 145 CALL dyn_cor_c1d( kstp ) ! vorticity term including Coriolis 146 CALL dyn_zdf ( kstp ) ! vertical diffusion 147 CALL dyn_nxt_c1d( kstp ) ! lateral velocity at next time step 144 148 145 149 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 146 150 ! Control and restarts 147 151 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 148 149 IF( kstp == nit000 ) CALL iom_close( numror )! close input ocean restart file150 IF( lrst_oce ) CALL rst_write ( kstp )! write output ocean restart file152 CALL stp_ctl( kstp, indic ) 153 IF( kstp == nit000 ) CALL iom_close( numror ) ! close input ocean restart file 154 IF( lrst_oce ) CALL rst_write( kstp ) ! write output ocean restart file 151 155 ! 156 #if defined key_iomput 157 IF( kstp == nitend .OR. indic < 0 ) CALL xios_context_finalize() ! needed for XIOS 158 ! 159 #endif 152 160 END SUBROUTINE stp_c1d 153 161
Note: See TracChangeset
for help on using the changeset viewer.