source:
branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/NEMO/SAO_SRC/sao_intp.F90
@
5042
| Last change on this file since 5042 was 4852, checked in by andrewryan, 10 years ago | |
|---|---|
| File size: 1.6 KB | |
| Line | |
|---|---|
| 1 | MODULE sao_intp |
| 2 | !!====================================================================== |
| 3 | !! *** MODULE sao_intp *** |
| 4 | !! ** Purpose : Run NEMO observation operator in offline mode |
| 5 | !!====================================================================== |
| 6 | !! NEMO modules |
| 7 | USE in_out_manager |
| 8 | USE diaobs |
| 9 | !! Stand Alone Observation operator modules |
| 10 | USE sao_read |
| 11 | USE sao_data |
| 12 | |
| 13 | IMPLICIT NONE |
| 14 | PRIVATE |
| 15 | |
| 16 | PUBLIC sao_interp |
| 17 | |
| 18 | CONTAINS |
| 19 | |
| 20 | SUBROUTINE sao_interp |
| 21 | !!---------------------------------------------------------------------- |
| 22 | !! *** SUBROUTINE sao_interp *** |
| 23 | !! |
| 24 | !! ** Purpose : To interpolate the model as if it were running online. |
| 25 | !! |
| 26 | !! ** Method : 1. Populate model counterparts |
| 27 | !! 2. Call dia_obs at appropriate time steps |
| 28 | !!---------------------------------------------------------------------- |
| 29 | INTEGER :: istp ! time step index |
| 30 | !! Loop over entire run |
| 31 | istp = nit000 - 1 |
| 32 | nstop = 0 |
| 33 | DO WHILE ( istp <= nitend .AND. nstop == 0 ) |
| 34 | IF (jifile <= n_files + 1) THEN |
| 35 | IF ( MOD(istp, nn_sao_freq) == nit000 ) THEN |
| 36 | !! Read next model counterpart |
| 37 | CALL sao_rea_dri(jifile) |
| 38 | jifile = jifile + 1 |
| 39 | ENDIF |
| 40 | !! Interpolate single time step |
| 41 | CALL dia_obs(istp) |
| 42 | ENDIF |
| 43 | !! Increment model step |
| 44 | istp = istp + 1 |
| 45 | END DO |
| 46 | END SUBROUTINE sao_interp |
| 47 | |
| 48 | END MODULE sao_intp |
Note: See TracBrowser
for help on using the repository browser.
