Changeset 3618 for palm/trunk/UTIL/inifor/src/inifor_grid.f90
- Timestamp:
- Dec 10, 2018 1:25:22 PM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/UTIL/inifor/src/inifor_grid.f90
r3615 r3618 26 26 ! ----------------- 27 27 ! $Id$ 28 ! Prefixed all INIFOR modules with inifor_, removed unused variables 29 ! 30 ! 31 ! 3615 2018-12-10 07:21:03Z raasch 28 32 ! bugfix: abort replaced by inifor_abort 29 33 ! … … 97 101 !------------------------------------------------------------------------------! 98 102 99 MODULE grid100 101 USE control102 USE defs,&103 MODULE inifor_grid 104 105 USE inifor_control 106 USE inifor_defs, & 103 107 ONLY: DATE, EARTH_RADIUS, TO_RADIANS, TO_DEGREES, PI, dp, hp, sp, & 104 108 SNAME, LNAME, PATH, FORCING_STEP, WATER_ID, FILL_ITERATIONS, & 105 109 BETA, P_SL, T_SL, BETA, RD, RV, G, P_REF, RD_PALM, CP_PALM, & 106 110 RHO_L, OMEGA, HECTO 107 USE i o,&111 USE inifor_io, & 108 112 ONLY: get_netcdf_attribute, get_netcdf_dim_vector, & 109 113 get_netcdf_variable, parse_command_line_arguments, & 110 114 validate_config 111 USE netcdf, & 112 ONLY: NF90_MAX_NAME, NF90_MAX_VAR_DIMS 113 USE transform, & 115 USE inifor_transform, & 114 116 ONLY: average_2d, rotate_to_cosmo, find_horizontal_neighbours,& 115 117 compute_horizontal_interp_weights, & … … 119 121 project, centre_velocities, phi2phirot, rla2rlarot, uv2uvrot, & 120 122 phirot2phi, rlarot2rla 121 USE types 122 USE util 123 USE inifor_types 124 USE inifor_util 125 USE netcdf, & 126 ONLY: NF90_MAX_NAME, NF90_MAX_VAR_DIMS 123 127 124 128 IMPLICIT NONE … … 628 632 phirot = project(0.5_dp*ly, y0, EARTH_RADIUS) * TO_DEGREES, & 629 633 rlarot = project(0.5_dp*lx, x0, EARTH_RADIUS) * TO_DEGREES, & 630 polphi = phi_cn * TO_DEGREES, pollam = lambda_cn * TO_DEGREES,&634 polphi = phi_cn * TO_DEGREES, & 631 635 polgam = gam * TO_DEGREES & 632 636 ) * TO_RADIANS … … 660 664 f3 = 2.0_dp * OMEGA * SIN( & 661 665 TO_RADIANS*phirot2phi( phi_centre * TO_DEGREES, lam_centre * TO_DEGREES,& 662 phi_n * TO_DEGREES, lambda_n * TO_DEGREES,&666 phi_n * TO_DEGREES, & 663 667 gam * TO_DEGREES ) & 664 668 ) … … 1825 1829 IF ( number_stretch_level_end >= 1 ) THEN 1826 1830 CALL calculate_stretching_factor( number_stretch_level_end, dz, & 1827 dz_stretch_factor, &1828 1831 dz_stretch_factor_array, & 1829 1832 dz_stretch_level_end, & … … 1903 1906 !> results into an overdetermined system. 1904 1907 !------------------------------------------------------------------------------! 1905 SUBROUTINE calculate_stretching_factor( number_end, dz, dz_stretch_factor,&1908 SUBROUTINE calculate_stretching_factor( number_end, dz, & 1906 1909 dz_stretch_factor_array, & 1907 1910 dz_stretch_level_end, & … … 1911 1914 REAL(dp), DIMENSION(:), INTENT(INOUT) :: dz_stretch_factor_array 1912 1915 REAL(dp), DIMENSION(:), INTENT(IN) :: dz_stretch_level_end, dz_stretch_level_start 1913 REAL(dp) :: dz_stretch_factor1914 1916 1915 1917 INTEGER :: iterations !< number of iterations until stretch_factor_lower/upper_limit is reached … … 4117 4119 END SUBROUTINE get_soil_layer_thickness 4118 4120 4119 END MODULE grid4121 END MODULE inifor_grid
Note: See TracChangeset
for help on using the changeset viewer.