Changeset 8531 for branches/ORCHIDEE_2_2
- Timestamp:
- 2024-04-15T11:23:37+02:00 (7 months ago)
- Location:
- branches/ORCHIDEE_2_2/ORCHIDEE
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/ORCHIDEE_2_2/ORCHIDEE/src_driver/forcing_tools.f90
r7901 r8531 4484 4484 REAL(r_std), ALLOCATABLE, DIMENSION(:,:,:) :: tmp_slab2d 4485 4485 CHARACTER(LEN=80) :: name 4486 CHARACTER(LEN=80) :: coordinates ! temporary variable for attribut coordinates 4486 4487 LOGICAL :: windzero 4487 4488 ! … … 4544 4545 iret = NF90_GET_ATT(force_id(fileindex), varid, 'cell_methods', cellmethod) 4545 4546 IF (iret /= NF90_NOERR) THEN 4546 ! If the attribute is not found then we set a reasonable default : instantaneous and centered.4547 cellmethod=" time: instantaneous"4547 ! If the attribute is not found then we set undefined 4548 cellmethod="undefined" 4548 4549 ENDIF 4550 4551 ! Read also attribute coordinates that might contain information about the name of the time axis. 4552 ! If the forcing file is produced by XIOS, this is the case. 4553 ! Reading of the time axis name can fail if the attribute coordinates contains something else than the time axis name. 4554 iret = NF90_GET_ATT(force_id(fileindex), varid, 'coordinates', coordinates) 4555 IF (iret == NF90_NOERR) THEN 4556 ! The attribute 'coordinates' was found. This is probably containg the name of the time axis. 4557 ! We add this information in the begining of the celllmethod variable. 4558 cellmethod=TRIM(coordinates)//":"//TRIM(cellmethod) 4559 ENDIF 4560 4549 4561 ! 4550 4562 ! … … 4645 4657 ENDDO 4646 4658 ! 4647 ! Go through all the time axes we have to find the right one. 4659 ! Go through all the time axes we have to find the right one based on the information 4660 ! read from attribute cell_methods or coordinates 4648 4661 ! 4649 4662 timeindex=0 … … 4682 4695 ENDDO 4683 4696 ! 4684 ! If there is no "(" then we have to find the centered axis.4685 4697 ELSE 4698 ! There is no "(" so we set the centered axis by defaut 4686 4699 DO im=1,nbtmethods 4687 4700 IF ( INDEX(time_cellmethod(itbase+im), "cent") > 0 ) THEN … … 4691 4704 ENDIF 4692 4705 ! 4693 ! The name of the time axis was found bu no method could be identified4706 ! The name of the time axis was found but no method could be identified 4694 4707 ! 4695 4708 IF ( timeindex < 1 ) THEN … … 4702 4715 ENDIF 4703 4716 ENDDO 4704 ! 4705 ! Should no corresponding time axis name be found, 4706 ! then we use the first centered one. 4707 ! 4708 itax=1 4709 DO WHILE ( timeindex < 1 ) 4710 IF ( INDEX(time_cellmethod(itax), "cent") > 0 ) THEN 4711 timeindex = itax 4712 ELSE 4713 itax = itax + 1 4714 ENDIF 4715 ENDDO 4716 ! 4717 4718 4719 ! 4720 ! If still no corresponding time axis name has be found, then we use the first centered one. 4721 ! This is the case if none of the attributes cell_methods or coordinates were found or if they 4722 ! didn't include the name of the corresponding time axis. 4723 IF ( timeindex < 1 ) THEN 4724 CALL ipslerr(2,'forcing_attributetimeaxe',& 4725 'No information about the name of the time axis was found in the attributes cell_method or coordinantes', & 4726 'The first centered time axis will be used as defaut.', & 4727 'Attribute coordinates is optional but must not contain something else than the name of the time axis.') 4728 itax=1 4729 DO WHILE ( timeindex < 1 ) 4730 IF ( INDEX(time_cellmethod(itax), "cent") > 0 ) THEN 4731 timeindex = itax 4732 ELSE 4733 itax = itax + 1 4734 ENDIF 4735 ENDDO 4736 END IF 4737 4717 4738 END SUBROUTINE forcing_attributetimeaxe 4718 4739 -
branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/sechiba.f90
r8462 r8531 739 739 IF (printlev_loc>=3) WRITE(numout,*) 'Start sechiba_main kjpindex =',kjpindex 740 740 741 CALL xios_orchidee_send_field("nav_lat",lalo(:,1)) 742 CALL xios_orchidee_send_field("nav_lon",lalo(:,2)) 743 741 744 !! 1. Initialize variables at each time step 742 745 CALL sechiba_var_init (kjpindex, rau, pb, temp_air) … … 1055 1058 CALL xios_orchidee_send_field("roughheight",roughheight) 1056 1059 CALL xios_orchidee_send_field("lai",lai) 1060 CALL xios_orchidee_send_field("u",u) 1061 CALL xios_orchidee_send_field("v",v) 1062 CALL xios_orchidee_send_field("zlev",zlev) 1063 1057 1064 histvar(:)=zero 1058 1065 DO ji = 1, kjpindex -
branches/ORCHIDEE_2_2/ORCHIDEE/src_xml/field_def_orchidee.xml
r8468 r8531 33 33 34 34 <field_group id="sechiba" grid_ref="grid_landpoints"> 35 <field id="nav_lon" name="nav_lon" long_name="Longitudes" unit="deg" operation="once" /> 36 <field id="nav_lat" name="nav_lat" long_name="Latitudes" unit="deg" operation="once" /> 35 37 <field id="areas" name="Areas" long_name="Mesh areas" unit="m2" operation="once" freq_offset="0ts" /> 36 38 <field id="LandPoints" name="LandPoints" long_name="Land Points" unit="1" operation="once" freq_offset="0ts"/> … … 65 67 <field id="precip_rain" name="precip_rain" long_name="Rainfall rate" unit="mm/s"/> 66 68 <field id="precip_snow" name="precip_snow" long_name="Snowfall rate" unit="mm/s"/> 69 <field id="u" name="u" long_name="Wind speed north-south" unit="m/s"/> 70 <field id="v" name="v" long_name="Wind speed east-west" unit="m/s"/> 71 <field id="zlev" name="zlev" long_name="Height of the lowest level of the atmospheric model" unit="m"/> 67 72 <field id="netrad" name="netrad" long_name="Net radiation" unit="W/m^2"/> 68 73 <field id="lai" name="lai" long_name="Leaf Area Index per PFT" unit="1" grid_ref="grid_nvm"/> -
branches/ORCHIDEE_2_2/ORCHIDEE/src_xml/file_def_orchidee.xml
r8468 r8531 738 738 </file> 739 739 740 <!-- Sechiba file 4 containing variables read and interpolated from input files-->741 <file id="sechiba4" name=" sechiba_interp_diag" output_level="1" output_freq="1d" enabled=".TRUE.">740 <!-- This file is written in the begining of the execution. The variables are written only once. --> 741 <file id="sechiba4" name="ORCHIDEE_interp_diag" output_level="3" output_freq="1d" enabled=".FALSE."> 742 742 <field_group id="groupremap_sechiba4" group_ref="remap_1d" grid_ref="grid_landpoints_out" > 743 743 … … 1001 1001 1002 1002 </file> 1003 1004 1005 <!-- ORCHIDEE offline forcing file --> 1006 1007 <!-- This is a file that can be used as a forcing file for ORCHIDEE offline simualtions --> 1008 <file id="forcing_by_ORCHIDEE" name="forcing_by_ORCHIDEE" output_level="1" output_freq="3h" enabled=".FALSE."> 1009 <field_group id="groupremap_forcing" group_ref="remap_3h" grid_ref="grid_landpoints_out" > 1010 1011 <!-- level 1 --> 1012 <field field_ref="contfrac" level="1"/> 1013 <field field_ref="nav_lon" level="1"/> 1014 <field field_ref="nav_lat" level="1"/> 1015 <field field_ref="areas" level="1"/> 1016 <field field_ref="temp_air" name="Tair" operation="instant" level="1"/> 1017 <field field_ref="pb" name="PSurf" operation="instant" level="1"/> 1018 <field field_ref="qair" name="Qair" operation="instant" level="1"/> 1019 <field field_ref="u" name="Wind_N" operation="instant" level="1"/> 1020 <field field_ref="v" name="Wind_E" operation="instant" level="1"/> 1021 <field field_ref="zlev" name="Levels" operation="instant" level="1"/> 1022 <field field_ref="precip_rain" name="Rainf" level="1"/> 1023 <field field_ref="precip_snow" name="Snowf" level="1"/> 1024 <field field_ref="lwdown" name="LWdown" level="1"/> 1025 <field field_ref="swdown" name="SWdown" level="1"/> 1026 </field_group> 1027 </file> 1028 1029 1003 1030 </file_definition>
Note: See TracChangeset
for help on using the changeset viewer.