Changeset 4990 for trunk/NEMOGCM/NEMO/OPA_SRC/C1D/step_c1d.F90
- Timestamp:
- 2014-12-15T17:42:49+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/C1D/step_c1d.F90
r4313 r4990 72 72 ! Ocean physics update (ua, va, ta, sa used as workspace) 73 73 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 74 CALL bn2( tsb, r n2b )! before Brunt-Vaisala frequency75 CALL bn2( tsn, r n2 )! now Brunt-Vaisala frequency74 CALL bn2( tsb, rab_b, rn2b ) ! before Brunt-Vaisala frequency 75 CALL bn2( tsn, rab_n, rn2 ) ! now Brunt-Vaisala frequency 76 76 ! VERTICAL PHYSICS 77 77 CALL zdf_bfr( kstp ) ! bottom friction … … 115 115 ! Passive Tracer Model 116 116 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 117 CALL trc_stp( kstp )! time-stepping117 CALL trc_stp( kstp ) ! time-stepping 118 118 #endif 119 119 … … 121 121 ! Active tracers (ua, va used as workspace) 122 122 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 123 tsa(:,:,:,:) = 0.e0! set tracer trends to zero123 tsa(:,:,:,:) = 0._wp ! set tracer trends to zero 124 124 125 CALL tra_sbc ( kstp )! surface boundary condition126 IF( ln_traqsr ) CALL tra_qsr ( kstp )! penetrative solar radiation qsr127 IF( ln_tradmp ) CALL tra_dmp ( kstp )! internal damping trends- tracers128 IF( lk_zdfkpp ) CALL tra_kpp ( kstp )! KPP non-local tracer fluxes129 CALL tra_zdf ( kstp )! vertical mixing130 CALL tra_nxt ( kstp ) ! tracer fields at next time step131 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 module125 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 eos( tsn, rhd, rhop, gdept_0(:,:,:) ) ! now potential density for zdfmxl 131 IF( ln_zdfnpc ) CALL tra_npc( kstp ) ! applied non penetrative convective adjustment on (t,s) 132 CALL tra_nxt( kstp ) ! tracer fields at next time step 133 133 134 134 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 135 135 ! Dynamics (ta, sa used as workspace) 136 136 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 137 ua(:,:,:) = 0.e0! set dynamics trends to zero138 va(:,:,:) = 0.e0137 ua(:,:,:) = 0._wp ! set dynamics trends to zero 138 va(:,:,:) = 0._wp 139 139 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 step140 IF( ln_dyndmp ) CALL dyn_dmp ( kstp ) ! internal damping trends- momentum 141 CALL dyn_cor_c1d( kstp ) ! vorticity term including Coriolis 142 CALL dyn_zdf ( kstp ) ! vertical diffusion 143 CALL dyn_nxt_c1d( kstp ) ! lateral velocity at next time step 144 144 145 145 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 146 146 ! Control and restarts 147 147 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 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 file148 CALL stp_ctl( kstp, indic ) 149 IF( kstp == nit000 ) CALL iom_close( numror ) ! close input ocean restart file 150 IF( lrst_oce ) CALL rst_write( kstp ) ! write output ocean restart file 151 151 ! 152 152 END SUBROUTINE stp_c1d
Note: See TracChangeset
for help on using the changeset viewer.