Changeset 2141 for branches/dev_1784_PATM/NEMO/OPA_SRC/step.F90

Timestamp:

2010-09-30T15:22:26+02:00 (14 years ago)

Author:

cbricaud

Message:

add module to read atmospheric pressure and module to compute surface pressure gradient

File:

• branches/dev_1784_PATM/NEMO/OPA_SRC/step.F90 (modified) (5 diffs)

branches/dev_1784_PATM/NEMO/OPA_SRC/step.F90

@@ -28 +28 @@
    USE oce             ! ocean dynamics and tracers variables
    USE dom_oce         ! ocean space and time domain variables 
+   USE sbc_oce         ! surface boundary condition
    USE zdf_oce         ! ocean vertical physics variables
    USE ldftra_oce      ! ocean tracer   - trends
@@ -42 +43 @@
    USE sbcrnf          ! surface boundary condition: runoff variables
    USE sbccpl          ! surface boundary condition: coupled formulation (call send at end of step)
+   USE sbcpatm         ! Atmospheric pressure forcing
    USE cpl_oasis3, ONLY : lk_cpl
 
@@ -70 +72 @@
    USE dynzdf          ! vertical diffusion               (dyn_zdf routine)
    USE dynspg_oce      ! surface pressure gradient        (dyn_spg routine)
+   USE dynapg          ! atmospheric pressure forcing     (dyn_apg routine)
    USE dynspg          ! surface pressure gradient        (dyn_spg routine)
    USE dynnxt          ! time-stepping                    (dyn_nxt routine)
@@ -168 +171 @@
 !      IF ( Agrif_Root() .and. lwp) Write(*,*) '---'
 !      IF (lwp) Write(*,*) 'Grid Number',Agrif_Fixed(),' time step ',kstp
-#endif   
+#endif
       indic = 1                                       ! reset to no error condition
 
@@ -310 +313 @@
                                CALL dyn_zdf( kstp )         ! vertical diffusion
                                indic=0
+        IF( ln_patm_apg  )     CALL dyn_apg( kstp )         ! Atmospheric pressure forcing
                                CALL dyn_spg( kstp, indic )  ! surface pressure gradient
                                CALL dyn_nxt( kstp )         ! lateral velocity at next time step