Changeset 5159 for LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat
- Timestamp:
- Aug 2, 2024, 9:58:25 PM (12 months ago)
- Location:
- LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat
- Files:
-
- 10 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/calfis.f90
r5158 r5159 1 1 ! $Id$ 2 2 3 ! 4 ! 3 4 5 5 SUBROUTINE calfis(lafin, & 6 6 jD_cur, jH_cur, & … … 25 25 pdqfi, & 26 26 pdpsfi) 27 ! 27 28 28 ! Auteur : P. Le Van, F. Hourdin 29 29 ! ......... … … 37 37 USE lmdz_ssum_scopy, ONLY: scopy, ssum 38 38 USE lmdz_comgeom2 39 USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm 40 USE lmdz_paramet 39 41 40 42 IMPLICIT NONE 41 43 !======================================================================= 42 ! 44 43 45 ! 1. rearrangement des tableaux et transformation 44 46 ! variables dynamiques > variables physiques 45 47 ! 2. calcul des termes physiques 46 48 ! 3. retransformation des tendances physiques en tendances dynamiques 47 ! 49 48 50 ! remarques: 49 51 ! ---------- 50 ! 52 51 53 ! - les vents sont donnes dans la physique par leurs composantes 52 54 ! naturelles. … … 64 66 ! ngridmx pour le pole sud 65 67 ! ---> ngridmx=2+(jjm-1)*iim 66 ! 68 67 69 ! Input : 68 70 ! ------- … … 74 76 ! pts surface temperature (K) 75 77 ! callrad clef d'appel au rayonnement 76 ! 78 77 79 ! Output : 78 80 ! -------- … … 81 83 ! pdhfi tendency for the potential temperature 82 84 ! pdtsfi tendency for the surface temperature 83 ! 85 84 86 ! pdtrad radiative tendencies \ both input 85 87 ! pfluxrad radiative fluxes / and output 86 ! 88 87 89 !======================================================================= 88 ! 90 89 91 !----------------------------------------------------------------------- 90 ! 92 91 93 ! 0. Declarations : 92 94 ! ------------------ 93 94 INCLUDE "dimensions.h"95 INCLUDE "paramet.h"96 95 97 96 INTEGER :: ngridmx … … 137 136 REAL :: zplev(ngridmx, llm + 1), zplay(ngridmx, llm) 138 137 REAL :: zphi(ngridmx, llm), zphis(ngridmx) 139 ! 138 140 139 REAL :: zrot(iip1, jjm, llm) ! AdlC May 2014 141 140 REAL :: zufi(ngridmx, llm), zvfi(ngridmx, llm) … … 143 142 REAL :: ztfi(ngridmx, llm), zqfi(ngridmx, llm, nqtot) 144 143 REAL :: zpk(ngridmx, llm) 145 ! 144 146 145 REAL :: pcvgu(ngridmx, llm), pcvgv(ngridmx, llm) 147 146 REAL :: pcvgt(ngridmx, llm), pcvgq(ngridmx, llm, 2) 148 ! 147 149 148 REAL :: zdufi(ngridmx, llm), zdvfi(ngridmx, llm) 150 149 REAL :: zdtfi(ngridmx, llm), zdqfi(ngridmx, llm, nqtot) 151 150 REAL :: zdpsrf(ngridmx) 152 ! 151 153 152 REAL :: zdufic(ngridmx, llm), zdvfic(ngridmx, llm) 154 153 REAL :: zdtfic(ngridmx, llm), zdqfic(ngridmx, llm, nqtot) … … 160 159 REAL :: zsinbis(iim), zcosbis(iim), z1bis(iim) 161 160 REAL :: unskap, pksurcp 162 ! 161 163 162 REAL :: flxwfi(ngridmx, llm) ! Flux de masse verticale sur la grille physiq 164 ! 163 165 164 LOGICAL, SAVE :: firstcal = .TRUE., debut = .TRUE. 166 165 ! REAL rdayvrai 167 166 168 ! 167 169 168 !----------------------------------------------------------------------- 170 ! 169 171 170 ! 1. Initialisations : 172 171 ! -------------------- 173 ! 174 ! 172 173 175 174 IF (firstcal) THEN 176 175 debut = .TRUE. … … 187 186 ENDIF ! of IF (firstcal) 188 187 189 ! 190 ! 188 189 191 190 !----------------------------------------------------------------------- 192 191 ! 40. transformation des variables dynamiques en variables physiques: … … 208 207 209 208 ! 42. pression intercouches et fonction d'Exner: 210 ! 209 211 210 ! ----------------------------------------------------------------- 212 211 ! .... zplev definis aux (llm +1) interfaces des couches .... … … 215 214 216 215 ! ... Exner = cp * ( p(l) / preff ) ** kappa .... 217 ! 216 218 217 unskap = 1. / kappa 219 ! 218 220 219 DO l = 1, llm 221 220 zpk(1, l) = ppk(1, 1, l) … … 241 240 ENDDO 242 241 zplev(ngridmx, llmp1) = pp(1, jjp1, llmp1) 243 ! 242 244 243 ! 245 244 … … 274 273 ! 43.bis traceurs 275 274 ! --------------- 276 ! 275 277 276 itr = 0 278 277 DO iq = 1, nqtot … … 337 336 ! ENDDO 338 337 339 ! 338 340 339 ! 45. champ u: 341 340 ! ------------ … … 452 451 zrfi(ngridmx, l) = 0. 453 452 ENDDO 454 ! 453 455 454 ! On change de grille, dynamique vers physiq, pour le flux de masse verticale 456 455 CALL gr_dyn_fi(llm, iip1, jjp1, ngridmx, flxw, flxwfi) … … 542 541 543 542 CALL gr_fi_dyn(1, ngridmx, iip1, jjp1, zdpsrf, pdpsfi) 544 ! 543 545 544 ! 62. enthalpie potentielle 546 545 ! --------------------- … … 587 586 ! initialisation des tendances 588 587 pdqfi(:, :, :, :) = 0. 589 ! 588 590 589 itr = 0 591 590 DO iq = 1, nqtot -
LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/lmdz_calfis_loc.F90
r5158 r5159 53 53 USE lmdz_ssum_scopy, ONLY: ssum 54 54 USE lmdz_comgeom2 55 55 USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm 56 USE lmdz_paramet 56 57 57 58 !======================================================================= … … 108 109 ! 0. Declarations : 109 110 ! ------------------ 110 111 INCLUDE "dimensions.h"112 INCLUDE "paramet.h"113 111 114 112 INTEGER :: ngridmx -
LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/mod_interface_dyn_phys.F90
r5134 r5159 13 13 SUBROUTINE Init_interface_dyn_phys 14 14 USE lmdz_phys_mpi_data 15 USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm 15 16 IMPLICIT NONE 16 INCLUDE 'dimensions.h' 17 17 18 18 19 INTEGER :: i,j,k -
LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/ce0l.F90
r5136 r5159 48 48 USE lmdz_comgeom2 49 49 50 USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm 51 USE lmdz_paramet 50 52 IMPLICIT NONE 51 53 52 54 !------------------------------------------------------------------------------- 53 55 ! Local variables: 54 INCLUDE "dimensions.h" 55 INCLUDE "paramet.h" 56 57 56 58 57 59 REAL :: masque(iip1, jjp1) !--- CONTINENTAL MASK -
LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/etat0dyn_netcdf.F90
r5136 r5159 45 45 USE lmdz_comgeom2 46 46 47 USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm 48 USE lmdz_paramet 47 49 IMPLICIT NONE 48 50 … … 50 52 PUBLIC :: etat0dyn_netcdf 51 53 52 INCLUDE "dimensions.h" 53 INCLUDE "paramet.h" 54 55 54 56 REAL, SAVE :: deg2rad 55 57 INTEGER, SAVE :: iml_dyn, jml_dyn, llm_dyn, ttm_dyn, fid_dyn -
LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/etat0phys_netcdf.F90
r5158 r5159 52 52 USE comconst_mod, ONLY: pi, dtvr 53 53 USE lmdz_iniprint, ONLY: lunout, prt_level 54 USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm 54 55 USE lmdz_comgeom2 55 56 USE lmdz_clesphys 57 USE lmdz_paramet 56 58 57 59 PRIVATE 58 60 PUBLIC :: etat0phys_netcdf 59 61 60 INCLUDE "dimensions.h" 61 INCLUDE "paramet.h" 62 62 63 INCLUDE "dimsoil.h" 63 64 REAL, SAVE :: deg2rad -
LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/iniphysiq_mod.F90
r5158 r5159 45 45 USE lmdz_tracstoke 46 46 47 USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm 48 USE lmdz_paramet 47 49 IMPLICIT NONE 48 50 … … 52 54 ! ======================================================================= 53 55 54 INCLUDE "dimensions.h" 55 INCLUDE "paramet.h" 56 57 56 58 57 59 REAL, INTENT (IN) :: prad ! radius of the planet (m) -
LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/init_ssrf_m.F90
r5136 r5159 3 3 !******************************************************************************* 4 4 5 USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm 5 6 USE indice_sol_mod, ONLY: is_ter, is_oce, is_oce, is_lic, epsfra 6 7 USE dimphy, ONLY: klon, zmasq … … 14 15 USE lmdz_iniprint, ONLY: lunout, prt_level 15 16 USE lmdz_comgeom2 17 USE lmdz_paramet 16 18 17 19 CHARACTER(LEN = 256), PARAMETER :: icefname = "landiceref.nc", icevar = "landice" 18 20 PRIVATE 19 21 PUBLIC :: start_init_subsurf 20 INCLUDE "dimensions.h"21 INCLUDE "paramet.h"22 22 23 23 CONTAINS -
LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/limit_netcdf.f90
r5136 r5159 84 84 USE lmdz_comgeom2 85 85 86 USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm 87 USE lmdz_paramet 86 88 IMPLICIT NONE 87 89 !------------------------------------------------------------------------------- 88 90 ! Arguments: 89 INCLUDE "dimensions.h" 90 INCLUDE "paramet.h" 91 92 91 93 REAL, DIMENSION(iip1, jjp1), INTENT(INOUT) :: masque ! land mask 92 94 REAL, DIMENSION(iip1, jjp1), INTENT(INOUT) :: phis ! ground geopotential … … 335 337 USE lmdz_comgeom2 336 338 339 USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm 340 USE lmdz_paramet 337 341 IMPLICIT NONE 338 INCLUDE "dimensions.h" 339 INCLUDE "paramet.h" 342 343 340 344 !----------------------------------------------------------------------------- 341 345 ! Arguments: -
LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/test_disvert_m.F90
r5158 r5159 1 1 module test_disvert_m 2 2 3 USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm 3 4 IMPLICIT NONE 4 5 … … 19 20 20 21 ! For llm: 21 INCLUDE "dimensions.h" 22 22 23 23 24 ! Local:
Note: See TracChangeset
for help on using the changeset viewer.