[Dart-dev] [3258] DART/trunk/models/MITgcm_ocean/trans_pv_sv.f90: initial, barely compiling, structure not complete

[thoar at subversion.ucar.edu]

An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/dart-dev/attachments/20080314/5901b0af/attachment.html
-------------- next part --------------
Added: DART/trunk/models/MITgcm_ocean/trans_pv_sv.f90
===================================================================
--- DART/trunk/models/MITgcm_ocean/trans_pv_sv.f90	                        (rev 0)
+++ DART/trunk/models/MITgcm_ocean/trans_pv_sv.f90	2008-03-14 15:58:36 UTC (rev 3258)
@@ -0,0 +1,144 @@
+! Data Assimilation Research Testbed -- DART
+! Copyright 2004-2007, Data Assimilation Research Section
+! University Corporation for Atmospheric Research
+! Licensed under the GPL -- www.gpl.org/licenses/gpl.html
+ 
+program trans_pv_sv
+
+! <next few lines under version control, do not edit>
+! $URL$
+! $Id$
+! $Revision$
+! $Date$
+
+!----------------------------------------------------------------------
+! purpose: interface between MITgcm_ocean and DART
+!
+! method: Read MITgcm_ocean 'snapshot' files of model state
+!         Get target time from assim_model_state_ic (temp_ic).
+!         Reform fields into a state vector.
+!         Write out state vector in "proprietary" format for DART
+!
+! author: Tim Hoar 3/13/08
+!
+!----------------------------------------------------------------------
+
+use        types_mod, only : r4, r8
+use    utilities_mod, only : get_unit, file_exist, E_ERR, E_WARN, E_MSG, &
+                             initialize_utilities, finalize_utilities, &
+                             error_handler
+use        model_mod, only : MIT_meta_type, read_meta, read_snapshot, &
+                             drop_snapshot, prog_var_to_vector
+use  assim_model_mod, only : assim_model_type, static_init_assim_model, &
+                             init_assim_model, get_model_size, set_model_state_vector, &
+                             write_state_restart, set_model_time, open_restart_read, &
+                             open_restart_write, close_restart, aread_state_restart
+use time_manager_mod, only : time_type, read_time
+
+implicit none
+
+! version controlled file description for error handling, do not edit
+character(len=128), parameter :: &
+   source   = "$URL$", &
+   revision = "$Revision$", &
+   revdate  = "$Date$"
+
+character (len = 128) :: msgstring
+
+! eg. [S,T,U,V,SSH].0000040992.[data,meta]
+!
+! The '.meta' file contains matlab-format information about shapes, etc.
+
+character (len = 128) :: file_base = '0000000672'
+character (len = 128) :: S_filename
+character (len = 128) :: T_filename
+character (len = 128) :: U_filename
+character (len = 128) :: V_filename
+character (len = 128) :: SSH_filename
+
+character (len = 128) :: file_out = 'temp_ud', file_time = 'temp_ic'
+
+! Temporary allocatable storage to read in a native format for cam state
+type(assim_model_type) :: x
+!type(model_type)       :: var
+type(time_type)        :: model_time, adv_to_time
+
+real(r4), allocatable  :: S(:,:,:)
+real(r4), allocatable  :: T(:,:,:)
+real(r4), allocatable  :: U(:,:,:)
+real(r4), allocatable  :: V(:,:,:)
+real(r4), allocatable  :: SSH(:,:)
+
+real(r8), allocatable  :: x_state(:)
+
+integer                :: file_unit, x_size, timestep
+
+type(MIT_meta_type) :: mitmeta
+
+! The strategy here is to use the 'write_state_restart()' routine - 
+! which requires an 'assim_model_type' variable.  
+! Read the individual variables and pack them
+! into one big state vector, add a time type and make nn
+
+call initialize_utilities('trans_pv_sv')
+
+call read_snapshot('XC',SSH,timestep)
+
+! Static init assim model calls static_init_model
+call static_init_assim_model()
+call init_assim_model(x)
+
+write(*,*)'dimensions are ',shape(SSH)
+write(*,*)'model size is ',get_model_size()
+
+stop
+
+! Read the [meta,data] files 
+
+mitmeta = read_meta(file_base,'U')
+write(*,*)'timestep is ',timestep
+
+call read_snapshot(file_base,S,timestep,'S')
+call read_snapshot(file_base,T,timestep,'T')
+call read_snapshot(file_base,U,timestep,'U')
+call read_snapshot(file_base,V,timestep,'V')
+call read_snapshot(file_base,SSH,timestep,'SSH')
+
+! write(8)mitmeta%dimList
+! write(8)S
+! close(8)
+
+! Allocate the local state vector and fill it up.
+
+x_size = size(S) + size(T) + size(U) + size(V) + size(SSH)
+
+if ( x_size /= get_model_size() ) then
+   write(msgstring,*)'data size ',x_size,' /= ',get_model_size(),' model size'
+   call error_handler(E_ERR,'trans_pv_sv',msgstring,source,revision,revdate)
+endif
+
+allocate(x_state(x_size))
+
+! transform fields into state vector for DART
+call prog_var_to_vector(S,T,U,V,SSH,x_state)
+
+! Integration of model was controlled by the restart file,
+! so we use the target time of the restart file (from assim_model_state)
+! as the current model state time.
+ file_unit = open_restart_read(file_time)
+
+! We're done with x_state, so it can be uselessly filled in aread_state_restart,
+! while getting model_time.
+!call aread_state_restart(model_time, x_state, file_unit, adv_to_time)
+!call set_model_time (x, adv_to_time)
+!call close_restart(file_unit)
+
+! Get channel for output,
+! write out state vector in "proprietary" format
+!file_unit = open_restart_write(file_out)
+!call write_state_restart(x, file_unit)
+!call close_restart(file_unit)
+
+call finalize_utilities()
+
+end program trans_pv_sv


Property changes on: DART/trunk/models/MITgcm_ocean/trans_pv_sv.f90
___________________________________________________________________
Name: svn:keywords
   + Date Rev Author URL Id
Name: svn:eol-style
   + native