--- trunk/libf/dyn3d/leapfrog.f90 2008/04/18 14:45:53 10 +++ trunk/libf/dyn3d/leapfrog.f90 2010/03/03 13:23:49 24 @@ -1,41 +1,36 @@ module leapfrog_m - ! This module is clean: no C preprocessor directive, no include line. - IMPLICIT NONE contains - SUBROUTINE leapfrog(ucov, vcov, teta, ps, masse, phis, nq, q, clesphy0, & - time_0) + SUBROUTINE leapfrog(ucov, vcov, teta, ps, masse, phis, q, time_0) ! From dyn3d/leapfrog.F, version 1.6 2005/04/13 08:58:34 - ! Version du 10/01/98, avec coordonnees verticales hybrides, avec - ! nouveaux operat. dissipation * (gradiv2, divgrad2, nxgraro2) + ! Version du 10/01/98, avec coordonnées verticales hybrides, avec + ! nouveaux opérateurs dissipation "*" (gradiv2, divgrad2, nxgraro2) - ! Auteur: P. Le Van /L. Fairhead/F.Hourdin - ! Objet: - ! GCM LMD nouvelle grille + ! Auteurs : P. Le Van, L. Fairhead, F. Hourdin + ! Objet: nouvelle grille - ! ... Dans inigeom, nouveaux calculs pour les elongations cu, cv - ! et possibilite d'appeler une fonction f(y) a derivee tangente - ! hyperbolique a la place de la fonction a derivee sinusoidale. + ! Dans "inigeom", nouveaux calculs pour les élongations cu, cv + ! et possibilité d'appeler une fonction f(y) à dérivée tangente + ! hyperbolique à la place de la fonction à dérivée sinusoïdale. - ! ... Possibilité de choisir le schéma pour l'advection de - ! q, en modifiant iadv dans "traceur.def" (10/02) . + ! Possibilité de choisir le schéma pour l'advection de + ! q, en modifiant iadv dans "traceur.def". - ! Pour Van-Leer + Vapeur d'eau saturee, iadv(1)=4. (F.Codron, 10/99) + ! Pour Van-Leer + vapeur d'eau saturée, iadv(1)=4. ! Pour Van-Leer iadv=10 - use dimens_m, only: iim, jjm, llm, nqmx - use paramet_m, only: ip1jmp1, ip1jm, ijmllm, ijp1llm, jjp1, iip1, iip2 + use dimens_m, only: iim, llm, nqmx + use paramet_m, only: ip1jmp1, ip1jm, ijp1llm, jjp1, iip1 use comconst, only: dtvr, daysec, dtphys use comvert, only: ap, bp use conf_gcm_m, only: day_step, iconser, idissip, iphysiq, iperiod, nday, & offline, periodav - use logic, only: ok_guide, apdiss, apphys, conser, forward, iflag_phys, & - leapf, statcl + use logic, only: ok_guide, iflag_phys use comgeom use serre use temps, only: itaufin, day_ini, dt @@ -48,9 +43,6 @@ use pression_m, only: pression use pressure_var, only: p3d - integer nq - REAL clesphy0(:) - ! Variables dynamiques: REAL vcov(ip1jm, llm), ucov(ip1jmp1, llm) ! vents covariants REAL teta(ip1jmp1, llm) ! temperature potentielle @@ -59,7 +51,7 @@ REAL masse(ip1jmp1, llm) ! masse d'air REAL phis(ip1jmp1) ! geopotentiel au sol - REAL time_0 + REAL, intent(in):: time_0 ! Variables local to the procedure: @@ -95,9 +87,9 @@ REAL tppn(iim), tpps(iim), tpn, tps - INTEGER itau, itaufinp1 + INTEGER itau ! index of the time step of the dynamics, starts at 0 INTEGER iday ! jour julien - REAL time ! Heure de la journee en fraction d'1 jour + REAL time ! time of day, as a fraction of day length REAL SSUM real finvmaold(ip1jmp1, llm) @@ -106,9 +98,8 @@ INTEGER ij, l REAL rdayvrai, rdaym_ini - LOGICAL:: callinigrads = .true. - !+jld variables test conservation energie + ! Variables test conservation energie REAL ecin(ip1jmp1, llm), ecin0(ip1jmp1, llm) ! Tendance de la temp. potentiel d (theta) / d t due a la ! tansformation d'energie cinetique en energie thermique @@ -119,15 +110,13 @@ INTEGER:: ip_ebil_dyn = 0 ! PRINT level for energy conserv. diag. logical:: dissip_conservative = .true. - LOGICAL:: prem = .true. + logical forward, leapf, apphys, conser, apdiss !--------------------------------------------------- print *, "Call sequence information: leapfrog" itaufin = nday * day_step - itaufinp1 = itaufin + 1 - itau = 0 iday = day_ini time = time_0 @@ -136,21 +125,21 @@ iday = iday + 1 ENDIF + dq = 0. ! On initialise la pression et la fonction d'Exner : - dq=0. CALL pression(ip1jmp1, ap, bp, ps, p3d) CALL exner_hyb(ps, p3d, pks, pk, pkf) ! Debut de l'integration temporelle: outer_loop:do - if (ok_guide.and.(itaufin - itau - 1) * dtvr > 21600) then + if (ok_guide.and.(itaufin - itau - 1) * dtvr > 21600.) then call guide(itau, ucov, vcov, teta, q, masse, ps) else IF (prt_level > 9) print *, & - 'Attention : on ne guide pas les 6 dernieres heures.' + 'Attention : on ne guide pas les 6 dernières heures.' endif - CALL SCOPY(ijmllm, vcov, 1, vcovm1, 1) + CALL SCOPY(ip1jm * llm, vcov, 1, vcovm1, 1) CALL SCOPY(ijp1llm, ucov, 1, ucovm1, 1) CALL SCOPY(ijp1llm, teta, 1, tetam1, 1) CALL SCOPY(ijp1llm, masse, 1, massem1, 1) @@ -167,7 +156,6 @@ ! gestion des appels de la physique et des dissipations: apphys = .FALSE. - statcl = .FALSE. conser = .FALSE. apdiss = .FALSE. @@ -196,7 +184,8 @@ ! integrations dynamique et traceurs: CALL integrd(2, vcovm1, ucovm1, tetam1, psm1, massem1, dv, du, & - dteta, dq, dp, vcov, ucov, teta, q, ps, masse, phis, finvmaold) + dteta, dq, dp, vcov, ucov, teta, q, ps, masse, phis, & + finvmaold, leapf) ! calcul des tendances physiques: @@ -218,9 +207,9 @@ teta, q(:, :, 1), q(:, :, 2)) ENDIF - CALL calfis(nq, lafin, rdayvrai, time, ucov, vcov, teta, q, & + CALL calfis(nqmx, lafin, rdayvrai, time, ucov, vcov, teta, q, & masse, ps, pk, phis, phi, du, dv, dteta, dq, w, & - clesphy0, dufi, dvfi, dtetafi, dqfi, dpfi) + dufi, dvfi, dtetafi, dqfi, dpfi) ! ajout des tendances physiques: CALL addfi(nqmx, dtphys, & @@ -305,7 +294,7 @@ ENDIF ENDIF - IF (itau == itaufinp1) exit outer_loop + IF (itau == itaufin + 1) exit outer_loop ! ecriture du fichier histoire moyenne: @@ -319,7 +308,7 @@ ENDIF IF (itau == itaufin) THEN - CALL dynredem1("restart.nc", 0., vcov, ucov, teta, q, masse, ps) + CALL dynredem1("restart.nc", vcov, ucov, teta, q, masse, ps) CLOSE(99) ENDIF @@ -337,7 +326,7 @@ dt = 2. * dtvr END IF ELSE - ! ...... pas leapfrog ..... + ! pas leapfrog leapf = .TRUE. dt = 2. * dtvr END IF