Changeset 3421 for palm/trunk/SOURCE/ocean_mod.f90

Timestamp:

Oct 24, 2018 6:39:32 PM (6 years ago)

Author:

gronemeier

Message:

new surface-data output; renamed output variables (pt to theta, rho_air to rho, rho_ocean to rho_sea_water)

File:

• palm/trunk/SOURCE/ocean_mod.f90 (modified) (9 diffs)

palm/trunk/SOURCE/ocean_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Renamed output variables
+! 
+! 3381 2018-10-19 13:09:06Z raasch
 ! spin-up cooling for ocean surface implemented, see new parameter
 ! surface_cooling_spinup_time
@@ -674 +677 @@
     SELECT CASE ( TRIM( var ) )
 
-       CASE ( 'rho_ocean' )
+       CASE ( 'rho_sea_water' )
           unit = 'kg/m3'
 
@@ -723 +726 @@
           unit = dopr_unit
 
-       CASE ( 'rho_ocean' )
+       CASE ( 'rho_sea_water' )
           dopr_index(var_count) = 64
           dopr_unit             = 'kg/m3'
@@ -797 +800 @@
     SELECT CASE ( TRIM( var ) )
 
-       CASE ( 'rho_ocean', 'sa' )
+       CASE ( 'rho_sea_water', 'rho_sea_water_xy', 'rho_sea_water_xz', &
+              'rho_sea_water_yz', 'sa', 'sa_xy', 'sa_xz', 'sa_yz' )
           grid_x = 'x'
           grid_y = 'y'
@@ -846 +850 @@
        SELECT CASE ( TRIM( variable ) )
 
-          CASE ( 'rho_ocean' )
+          CASE ( 'rho_sea_water' )
              IF ( .NOT. ALLOCATED( rho_ocean_av ) )  THEN
                 ALLOCATE( rho_ocean_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
@@ -867 +871 @@
        SELECT CASE ( TRIM( variable ) )
 
-          CASE ( 'rho_ocean' )
+          CASE ( 'rho_sea_water' )
              IF ( ALLOCATED( rho_ocean_av ) )  THEN
                 DO  i = nxlg, nxrg
@@ -899 +903 @@
        SELECT CASE ( TRIM( variable ) )
 
-          CASE ( 'rho_ocean' )
+          CASE ( 'rho_sea_water' )
              IF ( ALLOCATED( rho_ocean_av ) ) THEN
                 DO  i = nxlg, nxrg
@@ -980 +984 @@
     SELECT CASE ( TRIM( variable ) )
 
-       CASE ( 'rho_ocean_xy', 'rho_ocean_xz', 'rho_ocean_yz' )
+       CASE ( 'rho_sea_water_xy', 'rho_sea_water_xz', 'rho_sea_water_yz' )
           IF ( av == 0 )  THEN
              to_be_resorted => rho_ocean
@@ -1073 +1077 @@
     SELECT CASE ( TRIM( variable ) )
 
-       CASE ( 'rho_ocean' )
+       CASE ( 'rho_sea_water' )
           IF ( av == 0 )  THEN
              to_be_resorted => rho_ocean