Changeset 4105 for LMDZ6/trunk


Ignore:
Timestamp:
Mar 17, 2022, 11:24:22 PM (2 years ago)
Author:
jyg
Message:

Bug fix in dyn1d/old_lmdz1d.F90: compute the geopotential before it is used

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/dyn1d/old_lmdz1d.F90

    r4104 r4105  
    10501050
    10511051!---------------------------------------------------------------------
    1052 ! Interpolation of forcings in time and onto model levels
    1053 !---------------------------------------------------------------------
    1054 
    1055 #include "old_1D_interp_cases.h"
    1056 
    1057       if (forcing_GCM2SCM) then
    1058         write (*,*) 'forcing_GCM2SCM not yet implemented'
    1059         stop 'in time loop'
    1060       endif ! forcing_GCM2SCM
    1061 
    1062 !---------------------------------------------------------------------
    10631052!  Geopotential :
    10641053!---------------------------------------------------------------------
     
    10691058     &    (play(l)-play(l+1))/(play(l)+play(l+1))
    10701059        enddo
     1060
     1061!---------------------------------------------------------------------
     1062! Interpolation of forcings in time and onto model levels
     1063!---------------------------------------------------------------------
     1064
     1065#include "old_1D_interp_cases.h"
     1066
     1067      if (forcing_GCM2SCM) then
     1068        write (*,*) 'forcing_GCM2SCM not yet implemented'
     1069        stop 'in time loop'
     1070      endif ! forcing_GCM2SCM
     1071
     1072!!!!---------------------------------------------------------------------
     1073!!!!  Geopotential :
     1074!!!!---------------------------------------------------------------------
     1075!!!
     1076!!!        phi(1)=RD*temp(1)*(plev(1)-play(1))/(.5*(plev(1)+play(1)))
     1077!!!        do l = 1, llm-1
     1078!!!          phi(l+1)=phi(l)+RD*(temp(l)+temp(l+1))*                           &
     1079!!!     &    (play(l)-play(l+1))/(play(l)+play(l+1))
     1080!!!        enddo
    10711081
    10721082!---------------------------------------------------------------------
Note: See TracChangeset for help on using the changeset viewer.