PISCES in OFFLINE and 1D mode

PISCES can be used in offline and 1D. At this point, I don't have a complete understanding of how the combination of offline and 1D differs from the standard PISCES (online & 3D) since many computations are expected to be simplified (i.e. neglected, or read from offline dynamics). Unfortunately, there is no formal documents for offline mode, while the 1D document is available here.

In offline mode, you need the following dynamical (physical) fields according to the "namelist" file:


!-----------------------------------------------------------------------
&namdta_dyn        !   offline dynamics read in files                ("key_offline")
!-----------------------------------------------------------------------
!            !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation !
!            !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  !
   sn_tem  = 'dyna_grid_T' ,    -12             , 'votemper' ,  .true.    , .true. ,   'yearly'  , ''       , ''
   sn_sal  = 'dyna_grid_T' ,    -12            , 'vosaline' ,  .true.    , .true. ,   'yearly'  , ''       , ''
   sn_mld  = 'dyna_grid_T' ,    -12            , 'somixhgt' ,  .true.    , .true. ,   'yearly'  , ''       , ''
   sn_emp  = 'dyna_grid_T' ,    -12            , 'sowaflcd' ,  .true.    , .true. ,   'yearly'  , ''       , ''
   sn_ice  = 'dyna_grid_T' ,    -12            , 'soicecov' ,  .true.    , .true. ,   'yearly'  , ''       , ''
   sn_qsr  = 'dyna_grid_T' ,    -12            , 'soshfldo' ,  .true.    , .true. ,   'yearly'  , ''       , ''
   sn_wnd  = 'dyna_grid_T' ,    -12            , 'sowindsp' ,  .true.    , .true. ,   'yearly'  , ''       , ''
   sn_uwd  = 'dyna_grid_U' ,    -12            , 'vozocrtx' ,  .true.    , .true. ,   'yearly'  , ''       , ''
   sn_vwd  = 'dyna_grid_V' ,    -12            , 'vomecrty' ,  .true.    , .true. ,   'yearly'  , ''       , ''
   sn_wwd  = 'dyna_grid_W' ,    -12            , 'vovecrtz' ,  .true.    , .true. ,   'yearly'  , ''       , ''
   sn_avt  = 'dyna_grid_W' ,    -12            , 'voddmavs' ,  .true.    , .true. ,   'yearly'  , ''       , ''
   sn_ubl  = 'dyna_grid_U' ,    -12            , 'sobblcox' ,  .true.    , .true. ,   'yearly'  , ''       , ''
   sn_vbl  = 'dyna_grid_V' ,    -12            , 'sobblcoy' ,  .true.    , .true. ,   'yearly'  , ''       , ''
   sn_eiw  = 'dyna_grid_W' ,    -12            , 'soleaeiw' ,  .true.    , .true. ,   'yearly'  , ''       , ''
!
   cn_dir      = './'      !  root directory for the location of the dynamical files
   ln_degrad   =  .false.  !  flag for degradation -                requires ("key_degrad")
   ln_dynwzv   =  .true.   !  computation of vertical velocity instead of using the one read in file
   ln_dynbbl   =  .false.   !  bbl coef are in files, so read them - requires ("key_trabbl")
/

Note that "sobblcox", "sobblcoy", and "vovecrtz" are not required if you set "ln_dynbbl" (for the first 2 variables) and "ln_dynwzv" equal ".false.". They will be computed in such case.
 
The offline dynamics are computed in "dtadyn.F90".

To create your own offline dynamics, you can first produce a reference dynamics file "dyna_grid_[TUVW].nc" by running OPA model (e.g. default ORCA_LIM2) for 1 year with the modified "iodef.xml":

      <group id="1y" output_freq="-12"    output_level="10" enabled=".TRUE.">                      <!-- real yearly files -->

        <file id="1y_grid_T" name="auto" description="ocean T grid variables" >

                <field ref="toce"               name="votemper" />

                <field ref="soce"               name="vosaline" />

                <field ref="mldkz5"             name="somixhgt" />

                <field ref="empsmr"             name="sowaflcd" />

                <field ref="ice_cover"  name="soicecov" />

                <field ref="qt"                 name="soshfldo" />

                <field ref="wspd"               name="sowindsp" />

        </file>

        <file id="1y_grid_U" name="auto" description="ocean U grid variables" >

                <field ref="uoce" name="vozocrtx"  />

        </file>

        <file id="1y_grid_V" name="auto" description="ocean V grid variables" >

                <field ref="voce" name="vomecrty" />

        </file>

        <file id="1y_grid_W" name="auto" desciption="ocean W grid variables" >

                <field ref="woce"         name="vovecrtz" />

                <field ref="avs"          name="voddmavs" />

                <field ref="aht2d_eiv"    name="soleaeiw" />

        </file>

      </group>

In Matlab, there is a NetCDF command called "ncwriteschema", which creates a file based on your scheme (the reference dynamics field). Now, assign the values to each variable, and you can use the file as your offline dynamics.


==============

The following netCDF files are used for offline 1D PISCES run:

bathy.orca.nc            data_O2_nomask.nc   dyna_grid_V.nc
data_Alkalini_nomask.nc  data_PO4_nomask.nc  dyna_grid_W.nc
data_DIC_nomask.nc       data_Si_nomask.nc   mesh_mask.nc
data_DOC_nomask.nc       dust.orca.nc        ndeposition.orca.nc
data_Fer_nomask.nc       dyna_grid_T.nc      presatm.nc
data_NO3_nomask.nc       dyna_grid_U.nc      river.orca.nc


ORCA2_OFF_PISCES on placentia (this one works!!!):
key_trabbl key_orca_r2 key_ldfslp key_traldf_c2d key_traldf_eiv key_top key_offline key_pisces key_iomput

OFFLINE_3D (this one does not work):
key_trabbl key_orca_r2 key_ldfslp key_traldf_c2d key_traldf_eiv key_top key_offline key_pisces key_iomput     key_diaeiv key_dynldf_c3d