coords_n_const.f90

Go to the documentation of this file.

!/===========================================================================/
! Copyright (c) 2007, The University of Massachusetts Dartmouth 
! Produced at the School of Marine Science & Technology 
! Marine Ecosystem Dynamics Modeling group
! All rights reserved.
!
! FVCOM has been developed by the joint UMASSD-WHOI research team. For 
! details of authorship and attribution of credit please see the FVCOM
! technical manual or contact the MEDM group.
!
! 
! This file is part of FVCOM. For details, see http://fvcom.smast.umassd.edu 
! The full copyright notice is contained in the file COPYRIGHT located in the 
! root directory of the FVCOM code. This original header must be maintained
! in all distributed versions.
!
! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" 
! AND ANY EXPRESS OR  IMPLIED WARRANTIES, INCLUDING,  BUT NOT  LIMITED TO,
! THE IMPLIED WARRANTIES OF MERCHANTABILITY AND  FITNESS FOR A PARTICULAR
! PURPOSE ARE DISCLAIMED.  
!
!/---------------------------------------------------------------------------/
! CVS VERSION INFORMATION
! $Id$
! $Name$
! $Revision$
!/===========================================================================/

SUBROUTINE coords_n_const

  USE mod_input
  USE mod_setup
  IMPLICIT NONE
  INTEGER STATUS, I

!===============================================================================!
! OPEN FILES
! UNLESS OTHERWISE SPECIFED THE ROUTINES CALLED HERE ARE IN mod_input.F
!===============================================================================!
    
  SELECT CASE(startup_type)
     !=================================================
     ! HOTSTART 
  CASE("hotstart")
     !=================================================
     if(dbg_set(dbg_log)) then 
        WRITE(ipt,*)'!               READING FILES FOR HOTSTART                       !'
        WRITE(ipt,*)'!                                                                !'
     end if
               
     ! ALLOCATE SPACE FOR THE LOCAL GRID DATA
     ALLOCATE(y_lcl(0:mt),stat=status)
     IF (status /=0 ) CALL fatal_error("COULD NOT ALLOCATE Y_LCL")
     y_lcl=0.0_sp
     ALLOCATE(x_lcl(0:mt),stat=status)
     IF (status /=0 ) CALL fatal_error("COULD NOT ALLOCATE X_LCL")
     x_lcl=0.0_sp
     ALLOCATE(h_lcl(0:mt),stat=status)
     IF (status /=0 ) CALL fatal_error("COULD NOT ALLOCATE H_LCL")
     h_lcl=0.0_sp

     CALL load_restart_coords(x_lcl,y_lcl)
     CALL coordinate_units(x_lcl,y_lcl)
     CALL setup_center_coords


     DEALLOCATE(x_lcl)
     DEALLOCATE(y_lcl)


     CALL load_restart_depth(h_lcl)
     CALL setup_depth
     DEALLOCATE(h_lcl) ! COULD BE LOADED DIRECTLY - MUST SET MAX/MIN

     CALL load_restart_coriolis(cor)! LOAD DIRECTLY

    !!  ggao for equatoral min (4deg)
!     IF(.NOT. EQUATOR_BETA_PLANE)THEN
!      WHERE(COR <  1.e-5_SP .AND. COR > 0.0_SP) COR =  1.e-5_SP
!      WHERE(COR > -1.e-5_SP .AND. COR < 0.0_SP) COR = -1.e-5_SP
!     END IF 

     CALL setup_gravity
     
     CALL load_restart_sponge(cc_sponge) ! LOAD DIRECTLY
     
     stype = stype_restart
     CALL load_restart_sigma(z,z1) ! LOAD DIRECTLY TO ALL_VARS:Z,Z1
     CALL setup_sigma_derivatives

     CALL setup_obctypes

     CALL setup_horizontal_mixing_coefficient

     CALL setup_bottom_roughness

     !=================================================
     ! CRASHSTART 
  CASE("crashrestart")
     !=================================================
     if(dbg_set(dbg_log)) then 
        WRITE(ipt,*)'!              READING FILES FOR CRASHRESTART                    !'
        WRITE(ipt,*)'!                                                                !'
     end if
     ! ALLOCATE SPACE FOR THE LOCAL GRID DATA
     ALLOCATE(y_lcl(0:mt),stat=status)
     IF (status /=0 ) CALL fatal_error("COULD NOT ALLOCATE Y_LCL")
     y_lcl=0.0_sp
     ALLOCATE(x_lcl(0:mt),stat=status)
     IF (status /=0 ) CALL fatal_error("COULD NOT ALLOCATE X_LCL")
     x_lcl=0.0_sp
     ALLOCATE(h_lcl(0:mt),stat=status)
     IF (status /=0 ) CALL fatal_error("COULD NOT ALLOCATE H_LCL")
     h_lcl=0.0_sp

          
     CALL load_restart_coords(x_lcl,y_lcl)
     CALL coordinate_units(x_lcl,y_lcl)
     CALL setup_center_coords


     DEALLOCATE(x_lcl)
     DEALLOCATE(y_lcl)

     
     CALL load_restart_depth(h_lcl)
     CALL setup_depth
     DEALLOCATE(h_lcl) ! COULD BE LOADED DIRECTLY - MUST SET MAX/MIN

     CALL load_restart_coriolis(cor)

     CALL setup_gravity
     
     CALL load_restart_sponge(cc_sponge) ! LOAD DIRECTLY TO ALL_VARS:CC_SPONGE


     stype = stype_restart
     CALL load_restart_sigma(z,z1) ! LOAD DIRECTLY TO ALL_VARS:Z,Z1
     CALL setup_sigma_derivatives

     CALL setup_obctypes

     CALL setup_horizontal_mixing_coefficient

     CALL setup_bottom_roughness

     !=================================================
     ! COLDSTART 
  CASE("coldstart")
     !=================================================
     if(dbg_set(dbg_log)) then 
        WRITE(ipt,*)'!              READING FILES FOR COLDSTART                       !'
        WRITE(ipt,*)'!                                                                !'
     end if
     
     IF(msr) THEN
        ! ALLOCATE SPACE FOR THE GLOBAL GRID DATA
        ALLOCATE(y_gbl(0:mgl),stat=status)
        IF (status /=0 ) CALL fatal_error("COULD NOT ALLOCATE Y_GBL")
        y_gbl=0.0_sp
        ALLOCATE(x_gbl(0:mgl),stat=status)
        IF (status /=0 ) CALL fatal_error("COULD NOT ALLOCATE X_GBL")
        x_gbl=0.0_sp
     END IF

    ! ALLOCATE SPACE FOR THE LOCAL GRID DATA
    ALLOCATE(y_lcl(0:mt),stat=status)
    IF (status /=0 ) CALL fatal_error("COULD NOT ALLOCATE Y_LCL")
    y_lcl=0.0_sp
    ALLOCATE(x_lcl(0:mt),stat=status)
    IF (status /=0 ) CALL fatal_error("COULD NOT ALLOCATE X_LCL")
    x_lcl=0.0_sp

    ALLOCATE(h_lcl(0:mt),stat=status)
    IF (status /=0 ) CALL fatal_error("COULD NOT ALLOCATE H_LCL")
    h_lcl=0.0_sp

    ALLOCATE(c_lcl(0:mt),stat=status)
    IF (status /=0 ) CALL fatal_error("COULD NOT ALLOCATE C_LCL")
    c_lcl=0.0_sp

     
     CALL load_coldstart_coords(x_gbl,y_gbl,x_lcl,y_lcl)
     CALL coordinate_units(x_lcl,y_lcl)
     CALL setup_center_coords


     DEALLOCATE(x_lcl)
     DEALLOCATE(y_lcl)

     CALL load_coldstart_depth(x_gbl,y_gbl,h_lcl)
     CALL setup_depth

     DEALLOCATE(h_lcl)

     CALL load_coldstart_coriolis(x_gbl,y_gbl,c_lcl)
     CALL setup_coriolis

     DEALLOCATE(c_lcl)

     CALL setup_gravity

     ! THESE VARIABLES ARE NOT ALLOCATED YET
     CALL load_coldstart_sponge(x_gbl,y_gbl,nsponge,n_spg,r_spg,c_spg,x_spg,y_spg)
     CALL setup_sponge

 
!     CALL LOAD_COLDSTART_SIGMA
     CALL setup_sigma
     CALL setup_sigma_derivatives

     IF(msr) THEN
        DEALLOCATE(x_gbl)
        DEALLOCATE(y_gbl)
     END IF

     CALL setup_obctypes

     CALL setup_horizontal_mixing_coefficient

     CALL setup_bottom_roughness

  END SELECT
  
  
END SUBROUTINE coords_n_const