New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 2129 for branches – NEMO

Changeset 2129 for branches


Ignore:
Timestamp:
2010-09-29T13:01:23+02:00 (14 years ago)
Author:
rfurner
Message:

pervious merge had mised bits of OBS branch, adding these now

Location:
branches/devukmo2010/NEMO/OPA_SRC
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/devukmo2010/NEMO/OPA_SRC/SBC/geo2ocean.F90

    r1833 r2129  
    2626                                             ! repcmo and repere are keep only for compatibility. 
    2727                                             ! they are only a useless overlay of rot_rep 
     28 
     29   PUBLIC   obs_rot 
    2830 
    2931   REAL(wp), DIMENSION(jpi,jpj) ::   & 
     
    522524   END SUBROUTINE repere 
    523525 
     526 
     527   SUBROUTINE obs_rot ( psinu, pcosu, psinv, pcosv ) 
     528      !!---------------------------------------------------------------------- 
     529      !!                  ***  ROUTINE obs_rot  *** 
     530      !! 
     531      !! ** Purpose :   Copy gsinu, gcosu, gsinv and gsinv 
     532      !!                to input data for rotations of 
     533      !!                current at observation points 
     534      !! 
     535      !! History : 
     536      !!   9.2  !  09-02  (K. Mogensen) 
     537      !!---------------------------------------------------------------------- 
     538      REAL(wp), DIMENSION(jpi,jpj), INTENT( OUT )::   & 
     539         & psinu, pcosu, psinv, pcosv! copy of data 
     540 
     541      !!---------------------------------------------------------------------- 
     542 
     543      ! Initialization of gsin* and gcos* at first call 
     544      ! ----------------------------------------------- 
     545 
     546      IF( lmust_init ) THEN 
     547         IF(lwp) WRITE(numout,*) 
     548         IF(lwp) WRITE(numout,*) ' obs_rot : geographic <--> stretched' 
     549         IF(lwp) WRITE(numout,*) ' ~~~~~~~   coordinate transformation' 
     550 
     551         CALL angle       ! initialization of the transformation 
     552         lmust_init = .FALSE. 
     553 
     554      ENDIF 
     555 
     556      psinu(:,:) = gsinu(:,:) 
     557      pcosu(:,:) = gcosu(:,:) 
     558      psinv(:,:) = gsinv(:,:) 
     559      pcosv(:,:) = gcosv(:,:) 
     560 
     561   END SUBROUTINE obs_rot 
     562 
     563 
    524564  !!====================================================================== 
    525565END MODULE geo2ocean 
  • branches/devukmo2010/NEMO/OPA_SRC/opa.F90

    r2128 r2129  
    5858   USE sshwzv          ! vertical velocity used in asm 
    5959   USE diaptr          ! poleward transports           (dia_ptr_init routine) 
     60   USE diaobs          ! Observation diagnostics       (dia_obs_init routine) 
    6061   USE step            ! OPA time-stepping                  (stp     routine) 
    6162#if defined key_oasis3 
     
    160161         END DO 
    161162      ENDIF 
     163 
     164      IF( lk_diaobs ) CALL dia_obs_wri 
    162165        
    163166      !                            !------------------------! 
     
    294297      !                                     ! diagnostics 
    295298      CALL iom_init                             ! iom_put initialization 
     299 
     300      IF( lk_diaobs ) THEN 
     301         CALL dia_obs_init                  ! Initialize observational data 
     302         CALL dia_obs( nit000 - 1 )         ! Observation operator for restart 
     303      ENDIF       
     304 
    296305      CALL dia_ptr_init                         ! Poleward TRansports initialization 
    297306      CALL trd_mod_init                         ! Mixed-layer/Vorticity/Integral constraints trends 
  • branches/devukmo2010/NEMO/OPA_SRC/step.F90

    r2128 r2129  
    111111   USE diahth          ! thermocline depth                (dia_hth routine) 
    112112   USE diafwb          ! freshwater budget                (dia_fwb routine) 
     113   USE diaobs          ! obs-minus-model (assimilation)   (dia_obs routine) 
    113114   USE flo_oce         ! floats variables 
    114115   USE floats          ! floats computation               (flo_stp routine) 
     
    250251      IF( lk_diahdy  )   CALL dia_hdy( kstp )         ! dynamical heigh diagnostics 
    251252      IF( lk_diafwb  )   CALL dia_fwb( kstp )         ! Fresh water budget diagnostics 
     253      IF( lk_diaobs  )   CALL dia_obs( kstp )         ! obs-minus-model (assimilation) diagnostics 
    252254      IF( ln_diaptr  )   CALL dia_ptr( kstp )         ! Poleward TRansports diagnostics 
    253255      IF( lk_diaar5  )   CALL dia_ar5( kstp )         ! ar5 diag 
Note: See TracChangeset for help on using the changeset viewer.